336913_zh-CN

liliconcom

泰勒级数线性化以NHANES计算方差

你好，

通过结合人口加权和人口阶层与PSU的NHANES数据，是否有办法在JMP中进行调查设计？ 下面的CDC附带了SAS代码以及R＆Stata代码。

https://wwwn.cdc.gov/nchs/nhanes/tutorials/Module4.aspx

码：
SAS

PROC SURVEYMEANS data=one varmethod=taylor nomcar;
  STRATA sdmvstra;
  CLUSTER sdmvpsu;
  WEIGHT WTMEC4YR;
  DOMAIN Select;
   * more statements...;
run;

理想情况下，我可以完成此操作，以便可以运行在JMP中非常常见的制表功能。 
 

脚本编写回复：泰勒级数线性化以NHANES计算方差

我不确定您是否可以仅在JMP / JSL中执行此操作，但是JMP使得运行R代码如此方便，我想知道是否值得在这里重新发明轮子。

我没有对此进行测试，但是在安装R之后，您的代码将如下所示：

 

Names default to here(1);

//Open R
R Init();

//Open your data table here:
dt = Current data table();

//send data to R
R Send(dt);

//run your code
R Submit( "\[
 # install package if needed, really only needs to happen once
 if (!("survey" %in% installed.packages())) install.packages("survey")
 
 #load package
 library(survey)

 #code from the website, you will want to verify column names here
 NHANES_all <- svydesign(data=dt, id=~SDMVPSU, strata=~SDMVSTRA, weights=~WTMEC4YR, nest=TRUE)                    
 NHANES <- subset(NHANES_all, inAnalysis==1)
 m <- svymean(~Depression, NHANES)
 
]\" );

//return data from R
svymean = R Get( m )

//Close R
R Term();

有关将R与JMP集成的更多信息：

https://www.jmp.com/en_au/events/ondemand/mastering-jmp/jmp-and-r-integration.html

https://community.jmp.com/t5/JMPer-Cable/Getting-started-with-the-JMP-to-R-Interface/ba-p/43920

 

以下是他们使用的R包的一些信息：

https://cran.r-project.org/web/packages/survey/survey.pdf