16得票2回答
使用NeweyWest时如何更新摘要?

我正在使用NeweyWest标准误差来纠正我的lm() / dynlm()输出。例如:fit1<-dynlm(depvar~covariate1+covariate2) coeftest(fit1,vcov=NeweyWest) 系数显示得很好,但不幸的是,我失去了所有回归输出信息,比如...

16得票3回答
使用ggpmisc显示nls模型的方程式

R包ggpmisc可以用来在ggplot2上显示lm模型和poly模型的方程(参见这里)。 我想知道如何使用ggpmisc在ggplot2上显示nls模型方程结果。下面是我的最小工作示例。library(ggpmisc) args <- list(formula = y ~ k * e ...

16得票6回答
R错误提示说:“模型未全部拟合到相同大小的数据集上”。

我已经创建了两个广义线性模型,如下所示:glm1 <-glm(Y ~ X1 + X2 + X3, family=binomial(link=logit)) glm2 <-glm(Y ~ X1 + X2, family=binomial(link=logit)) 然后我使用anov...

16得票2回答
在使用boxcox公式时,lm.fit(x, y, offset = offset, singular.ok, ...)出现错误,有0个非NA值。

我正在尝试用以下代码运行一个boxcox变换:urban1 <- subset(ski,urban <= 4,na.rm=TRUE) ski$gender <- as.numeric((as.character(ski$gender)),na.rm=TRUE) urban1 ...

15得票4回答
如何使用Anova命令(car包)进行Tukey HSD测试

我在处理一个不平衡的设计/样本,并最初学习了 aov()。现在我知道,对于我的ANOVA测试,我需要使用第III型平方和,这需要使用 lm() 进行拟合,而不是使用 aov()。 问题是如何使用 lm() 来获取事后检验(具体来说是 Tukey's HSD)。我所做的所有研究都表明,在 mu...

15得票2回答
数据框替换时发生错误,*tmp* 数据有 x 条数据,但要替换的数据有 y 条。

我是一名R语言的初学者。这里有一个非常简单的代码,我试图保存残差项: # Create variables for child's EA: dat$cldeacdi <- rowMeans(dat[,c('cdcresp', 'cdcinv')],na.rm=T) dat$cldea...

14得票4回答
在R的glmnet中,使用$\lambda=0$的LASSO和OLS会产生不同的结果。

我预期LASSO在没有惩罚因子($\lambda=0$)的情况下会产生与OLS拟合相同(或非常相似)的系数估计值。然而,当我将相同的数据(x,y)放入以下两个模型中: glmnet(x, y , alpha=1, lambda=0) 用于LASSO拟合且不使用惩罚因子 lm(y ~ x) ...

14得票3回答
如何在R中使用公式排除主效应但保留交互作用。

我不需要主效应,因为它与一个更细的因素固定效应共线,所以出现这些 NA 很烦人。在这个例子中:lm(y ~ x * z) 我希望得到x(数值型)和z(因子型)之间的交互作用,但不要得到z的主效应。

14得票3回答
在R中拟合函数

我有一些数据点(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...

14得票1回答
将一个或多个系数设置为特定整数。

我正在使用标准的lm模型,并希望将一个或多个变量的系数设置为特定的整数。 例如,我希望天气和价格变量的系数分别设置为647和15。 我使用带有标准公式的lm函数。 到目前为止,我找到的最接近的解决方法是glm中的offset函数或systemfit中的restrict.rhs函数。 我...