问题描述
我有需要分析的连续（非计数）响应数据。此数据（请参阅给定可重复数据中此数据的摘录）约有 58% 的响应变量 (ShannonDiv) 为零。如果是计数数据，我将使用零膨胀/障碍模型。是否有可以处理此类数据的等效方法。我曾尝试阅读，但恐怕我没有成功。
df<-structure(list(SiteID = structure(c(27L,28L,29L,3L,30L,31L,32L,33L,34L,35L,9L,6L,7L,8L,4L),.Label = c("1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24","25","26","27","28","29","30","31","32","33","34","35"),class = "factor"),WindSpeed = c(1.7,0.3,2.1,0.6,0.8,1.7,2.5,3.4,4.8,1.2,0.4,2,1.6,2.6,0.5,1.4,2.3,4.1,1.3,0.7,3.8,0.1,0.2,1.5,NA,3.7,2.6),Humidity = c(19.7,10.2,10.1,21.1,8.3,18.3,11.2,8.7,16.3,12.7,23.6,20.4,17.5,40.4,15.7,14.5,42.5,11.3,32.5,12.5,16.7,16.6,11.7,31,54.6,48.3,61.8,49.2,42.1,63.3,69.8,82.6,65.2,51.7,69.2,13.6),Temp = c(25.1,34.7,31.2,24.6,35,31.9,25,26.3,28.2,24.1,34.5,27,30.3,25.7,31.3,33,24.3,29.9,30.9,22,34,33.6,33.2,23.8,28.7,29,27.1,21.7,25.3,30.4,23.1,26.4),ShannonDiv = c(0,0.693147180559945,0.934769897858279)),class = c("data.table","data.frame"),row.names = c(NA,-41L),.internal.selfref = <pointer: 0x00000221ef9d1ef0>)
#如果它是零充气模型，我会使用
sh <- mixed_model(ShannonDiv ~ Windspeed+Humidity+Temp,data = df,random = ~ 1|SiteID,family = zi.poisson(),zi_fixed = ~Windspeed+Humidity+Temp)
summary(sh)
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）