预(误差比例减小)的计算意义

Question

我试图计算两个线性模型的PRE的重要性。我所有的变量都是连续的(所以我不能使用DAMisc包)。我试着用下面的R代码手工完成它，但是在测试它时，我一直得到NaN。我做错什么了？

library(psych)
modelC<-lm(reason.4~reason.16,data=iqitems)
modelA<-lm(reason.4~reason.16+reason.17,data=iqitems)

SSEC<-sum(modelC$residuals^2) #2019
SSEA<-sum(modelA$residuals^2) #1977
PRE=(SSEC-SSEA)/SSEC #0.02
pA<-2
pC<-1
n<-min(summary(modelA)$df)+1 #4
deg_freedom<-min(summary(modelA)$df) #3
Fvalue <- (PRE/(pA-pC))/((1-PRE)*(n-pA)) #0.0105
pf(Fvalue,pC-pA,deg_freedom,lower.tail=FALSE) #NaN

Answer 1

pf函数的第二个参数是f-统计量分子的自由度，第三个参数是f-统计量分母的自由度。

第二个论点的自由度是-1，现在是pC - pA = 1 - 2 = -1。这就是给你带来错误的原因，因为自由度是非负的，我对此没有经验，我认为你需要：

#second argument being pA - pC = 1
> pf(0.0105, 1, 3,lower.tail=FALSE)
[1] 0.9248493