0

私は2つの小さなデータセットを持っています:

infected.data.r.p <- structure(list(MLH = c(0.520408163265306, 0.436170212765957, 
0.344086021505376, 0.423076923076923, 0.406976744186047), ColGrowthCL_6 = c(5.923728814, 
0.283950617, 0.377358491, 1.728070175, 0.2)), .Names = c("MLH", 
"ColGrowthCL_6"), row.names = c("12", "22", "28", "30", "34"), class = "data.frame")

と

uninfected.sampling <- structure(list(MLH = c(0.524271844660194, 0.457446808510638, 
0.354838709677419, 0.398058252427184, 0.436893203883495), ColGrowthCL_6 = c(4.401639344, 
4.827586207, 6.387096774, 6.320754717, 4.225490196)), .Names = c("MLH", 
"ColGrowthCL_6"), row.names = c("218", "18", "21", "212", "99"
), class = "data.frame")

R で anova() 構文を使用してこれら 2 つのモデルを比較しようとすると (以下を参照)、p 値を生成できません。問題を引き起こしているのが 2 つのデータ セットの性質であるとは確信していません (ただし、2 つのデータ セットの構造の正確な違いにも興味があります)。 . ありがとうございました！

モデル比較構文:

infected.model<-glm(formula=as.formula(ColGrowthCL_6~MLH), family=poisson, infected.data.r.p)
uninfected.model<-glm(formula=as.formula(ColGrowthCL_6~MLH), family=poisson, uninfected.sampling)    

compare<-anova(infected.model,uninfected.model,test="Chisq")
print(compare)
summary(compare)
4

1 に答える 1

3

同じデータセットからのモデルしか比較できないと思います。そのため、2 つの (ネストされた) モデルをカイ 2 乗で比較する場合、同じデータセットからのものである必要があります。おそらく、p 値が計算されていないのはそのためです。

于 2013-08-24T23:05:50.097 に答える