问题描述
我有一个看起来像这样的df：
set.seed(42)
ID <- sample(1:30,100,rep=T)
Trait <- sample(0:1,rep=T)
Year <- sample(1992:1999,rep=T)
df <- cbind(ID,Trait,Year)
df <- as.data.frame(df)
ID是单个生物，特征是表型的存在/不存在，Year是进行观察的年份。
我想对个体之间的特征是否随机进行建模，像这样
library(MCMCglmm)
m <- MCMCglmm(Trait ~ ID,random = ~ Year,data = df,family = "categorical")
现在，想改组Trait列并运行x个置换，以检查我观察到的均值和CI是否超出随机预期的范围。 我可以使用for循环来做到这一点，但我宁愿使用tidyverse解决方案。 我读过lapply是bette（？）的替代品，但我一直在努力寻找足够具体的演练以供遵循。
我很感谢在这里提供的任何建议。
干杯！
杰米
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）