336913_zh-TW

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