小编DrX*_*DrX的帖子

我试图从ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/daily/ghcnd-stations.txt上读取气候站信息.但是,由于第一行未完全填充(缺少最后两个cols)并且第5列包含空格,因此我无法完成阅读:

fread('ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/daily/ghcnd-stations.txt',sep=)

它返回错误消息:

 Expected sep (' ') but new line, EOF (or other non printing character) ends 
 field 5 when detecting types from point 0: AGE00135039  35.7297    0.6500   
 50.0    ORAN-HOPITAL MILITAIRE     

如何在阅读此txt文件时正确应用fread？谢谢!

我有一个小数据，其中有 3 组（A、B、C），每组有 5 名参与者。所有这些参与者在 7 项不同的考试中每项都被测量 6 次，因此每个参与者总共得到 6*7=42 分。建立了简单的线性混合模型mylmm<-lmer(score ~1+group+exam+group*exam+(1|participant), data = mydata)。anova(mylmm)我可以使用多重比较函数获得分组、考试及其交互的方差分析结果和事后成对比较。

但数据量很小（只有5个参与者），残差mylmm也不正常，所以威力不足。我知道使用稳健的混合模型robustlmm和使用残差引导混合模型lmeresampler。但是，我无法使用这些方法进行方差分析和多重比较。有人可以帮我解决以下问题吗？真的很感激。

1. 是否有方法和可用的 R 包来执行线性混合模型的引导方差分析（和事后比较）？
2. 是否还需要计算引导程序或非参数方差分析的功效？如果可以的话，功率怎么计算呢？
3. 我能够使用simranova 方法来计算测试组、考试和模型对象交互的功效lme。另外，simr还可以用于查找事后成对比较的功效还是emmeans应该使用？谢谢