R中的方差分析用于提取有效代码

代码：

结果:

pred_model = anova(m1, m2, test="Chisq")
pred_model

---

我不熟悉R，有人能建议如何获取模型的重要代码吗？函数通常会返回一个包含一列p值的矩阵。然而，在这种情况下，

anova

的结果作为数据帧返回，该数据帧作为“anova”对象被进一步“class（）-ed”（以便它可以有自己的

print

-method）。运行

？lm

中的第一个示例，然后：

Analysis of Variance Table

Model 1: male_birth ~ male_death + female_birth + female_death
Model 2: male_birth ~ male_death + female_birth
Res.Df    RSS Df  Sum of Sq Pr(>Chi)
 1     48 3.4883                       
 2     49 3.4951 -1 -0.0068189   0.7594

---

summary

函数通常会返回一个包含一列p值的矩阵。然而，在这种情况下，

anova

的结果作为数据帧返回，该数据帧作为“anova”对象被进一步“class（）-ed”（以便它可以有自己的

print

-method）。运行

？lm

中的第一个示例，然后：

Analysis of Variance Table

Model 1: male_birth ~ male_death + female_birth + female_death
Model 2: male_birth ~ male_death + female_birth
Res.Df    RSS Df  Sum of Sq Pr(>Chi)
 1     48 3.4883                       
 2     49 3.4951 -1 -0.0068189   0.7594

---

什么意思？您想获取什么？@user20650

##Signif。代码：0'***'0.001'***'0.01'*'0.05.'0.1''1.

p值。如果您查看

str（pred_模型）

您可以看到可以从

pred_模式

提取的对象的名称，因此对于pvalue，您可以使用

pred_模型$`Pr（>Chi）

。星星不是重新调谐的，而是印刷的。您可以捕获方差分析输出，或者搞乱产生这些结果的打印方法，编写自己的wee函数将pvalue转换为signfance星号，或者最好省略星号。事实上，从查看

printCoefmat

，您可以使用以下方法生成星号：如上面的注释所示提取pvalue，并将其分配给名为

pv

：

pv=pred_model$`Pr（>Chi）

，然后生成星号

symnum（pv，corr=FALSE，na=FALSE，cutpoints=c（0,0.001,0.01,0.05,0.1,1），symbols=c（“***”，“***”，“***”，“*”，”），“）

你的意思是什么？你想取什么？@user20650

##Signif.code:0'***'0.001'***'0.01'***'0.05.'0.1''1.

p值。如果你看一下

str（pred\u模型）

你可以看到可以从

pred\u模式

中提取的东西的名称，因此对于pvalue，你可以使用

pred\u模式$`Pr（>Chi）

。星星不是重新调谐的，而是打印出来的。你可以捕获方差分析输出，或者乱搞产生这些结果的打印方法，编写你自己的wee函数将pvalue转换为signfance星星，或者最好省略星星。事实上，通过查看

printCoefmat

，你可以生成星星：提取pvalue如上面的注释所示，分配给名为

pv

的对象：

pv=pred_model$`Pr（>Chi）

，然后生成星号

symnum（pv，corr=FALSE，na=FALSE，cutpoints=c（0，0.001，0.01，0.05，0.1，1），symbols=c（“***”，“***”，“***”，“*”，“，”）

anova(lm.D9,lm.0)$'Pr(>F)'[2]
[1] 0.2490232