我希望通过使用micein 来运行150次多次插补R.然而,为了节省一些计算时间,我将谎言将并行流中的流程细分(如Stef van Buuren在"针对缺失数据的灵活插补"中所建议的那样).
我的问题是:怎么做？
我可以想象2个选项:
opt.1:
imp1<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
imp2<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
imp...<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
imp150<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
然后通过使用complete和as.mids之后将插补结合在一起
opt.2:
imp1<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
imp2<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
imp...<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
imp150<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
通过添加VAL_1to150否则,在我看来(我可能是错的),如果它们都使用相同的数据集和相同的种子运行,那么您将获得相同结果的150倍.
还有其他选择吗？
谢谢
我正在使用R软件包TraMineR来计算和分析状态序列.我想使用命令seqfplot获得序列频率图.但是,不要使用设置要绘制的最频繁序列的数量
seqfplot(mydata.seq, tlim=1:20)
设置达到的最常见序列的百分比(例如,50%的样本)将是有用的.我试过这个
seqfplot(mydata.seq, trep = 0.5)
但是-不同于seqrep.grp和seqrep-选项trep不被支持的seqfplot命令.我应该创建一个新功能吗？
谢谢.