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 進行分析，然後將結果帶回 JMP 以使用表格/圖形產生器。

 

我還沒有對此進行測試，但是安裝 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