Displaying 1 result from an estimated 1 matches for "my_df2".
Did you mean:
my_df
2008 Oct 19
2
definition of "dffits"
...<- residuals(lm.out)
ey <- fitted(lm.out)
py <- ey + res1/(1-lev1)
df1 <- dffits(lm.out, infl = lm.influence(lm.out))
df1 <- dffits(lm.out)
print("df1: dffits")
print(df1)
my_df1 <- (ey-py)/(sig1*sqrt(lev1))
print("my_df1")
print(my_df1)
my_df2 <- -lev1*(ey-py)/(sig1*sqrt(lev1))
print("my_df2")
print(my_df2)
}
[1] "df1: dffits"
1 2 3 4 5
-1.3333333 0.4082483 0.6000000 -1.0475699 0.2672612
[1] "my_df1"
1 2 3 4...