我目前正在尝试使用bootstrapping在R中进行双向ANOVA。在我写的引导函数中,我遇到了一个试图调用引导系数的问题。我目前处理的模型正在研究两个变量:EtOH处理(3个水平)和特定轴突(4个水平)。我当前查看系数的函数是
summary.lm(lm(distance~treatment*axon)
但我得到下表:
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.35893 0.05247 25.901 < 2e-16 ***
treatmentChronic -0.60041 0.08109 -7.404 1.46e-12 ***
treatmentControl 0.13440 0.07002 1.919 0.0559 .
axonAxon 2 0.06613 0.07420 0.891 0.3735
axonAxon 3 0.19200 0.07420 2.588 0.0102 *
axonAxon 4 0.18240 0.07420 2.458 0.0146 *
treatmentChronic:axonAxon 2 -0.03650 0.11468 -0.318 0.7505
treatmentControl:axonAxon 2 0.01220 0.09903 0.123 0.9020
treatmentChronic:axonAxon 3 -0.08533 0.11468 -0.744 0.4574
treatmentControl:axonAxon 3 -0.02033 0.09903 -0.205 0.8375
treatmentChronic:axonAxon 4 0.02501 0.11468 0.218 0.8275
treatmentControl:axonAxon 4 0.05260 0.09903 0.531 0.5957
如您所见,Axon 1和急性治疗及其相互作用缺失。无论如何要解决这个问题,以便可以调用和引导值吗?