我正在使用NeweyWest标准误差来纠正我的lm() / dynlm()输出。例如:fit1<-dynlm(depvar~covariate1+covariate2) coeftest(fit1,vcov=NeweyWest) 系数显示得很好,但不幸的是,我失去了所有回归输出信息,比如...
我已经创建了两个广义线性模型,如下所示:glm1 <-glm(Y ~ X1 + X2 + X3, family=binomial(link=logit)) glm2 <-glm(Y ~ X1 + X2, family=binomial(link=logit)) 然后我使用anov...
我正在尝试用以下代码运行一个boxcox变换:urban1 <- subset(ski,urban <= 4,na.rm=TRUE) ski$gender <- as.numeric((as.character(ski$gender)),na.rm=TRUE) urban1 ...
我是一名R语言的初学者。这里有一个非常简单的代码,我试图保存残差项: # Create variables for child's EA: dat$cldeacdi <- rowMeans(dat[,c('cdcresp', 'cdcinv')],na.rm=T) dat$cldea...
我预期LASSO在没有惩罚因子($\lambda=0$)的情况下会产生与OLS拟合相同(或非常相似)的系数估计值。然而,当我将相同的数据(x,y)放入以下两个模型中: glmnet(x, y , alpha=1, lambda=0) 用于LASSO拟合且不使用惩罚因子 lm(y ~ x) ...
我不需要主效应,因为它与一个更细的因素固定效应共线,所以出现这些 NA 很烦人。在这个例子中:lm(y ~ x * z) 我希望得到x(数值型)和z(因子型)之间的交互作用,但不要得到z的主效应。
我有一些数据点(x和y),它们之间似乎存在对数关系。> mydata x y 1 0 123 2 2 116 3 4 113 4 15 100 5 48 87 6 75 84 7 122 77 > qplot(x, y, data=mydata...
我正在使用标准的lm模型,并希望将一个或多个变量的系数设置为特定的整数。 例如,我希望天气和价格变量的系数分别设置为647和15。 我使用带有标准公式的lm函数。 到目前为止,我找到的最接近的解决方法是glm中的offset函数或systemfit中的restrict.rhs函数。 我...