-
0 votesanswersviews
r:来自glmnet和插入符号的系数对于相同的lambda是不同的
我已经阅读了一些有关这方面的问答,但我仍然不确定我理解为什么基于相同样本和相同超参数的glmnet和插入符号模型的系数略有不同 . 非常感谢解释! 我正在使用插入符号来训练脊回归: library(ISLR) Hitters = na.omit(Hitters) x = model.matrix(Salary ~ ., Hitters)[, -1] #Dropping the intercept ... -
0 votesanswersviews
空间多边形数据帧的相关系数
我有一个空间多边形数据框,我对我的变量的相关系数矩阵感兴趣 . 命令 >cor(df) 返回以下错误: >Error in cor(MergedData) : supply both 'x' and 'y' or a matrix-like 'x' 如果我运行以下命令,我可以获得成对系数 >cor.test(df$var1, df$var2) 但是,由于我有15个变量,我需... -
0 votesanswersviews
如何在应用两个或更多修改系数的时段上拆分单个didgit?
例如,我知道我们在 previous month 中出售了 5000 件 . 我知道那个月分成几周,每个月我们 sell on 20% more than in previous month . 这是多年来的某种 global trend . 所以我可以通过使用这个方案每周知道确切的销售量:月分为4周,我们在那个月卖了5000件东西 . 第一周我们卖100%,第二个是105,第三个是110,最... -
0 votesanswersviews
在泊松GLM R中改变Y截距
Background :我有以下数据,我运行 glm 函数: location = c("DH", "Bos", "Beth") count = c(166, 57, 38) #make into df df = data.frame(location, count) #poisson summary(glm(count ~ loc... -
0 votesanswersviews
回归循环和存储系数
我将(1)多次循环回归某个标准; (2)从每个回归中存储一定的系数 . 这是一个例子: clear sysuse auto.dta local x = 2000 while `x' < 5000 { xi: regress price mpg length gear_ratio i.foreign if weight < `x' est sto model_`x... -
0 votesanswersviews
计算Stata中的t检验
我目前正在尝试对变量进行t检验,并确定它是否与1在统计上有显着差异 . 以下是我使用的代码: ttest dm1=1 它正在吐出这个输出: 我不希望我的零假设是mean = 1,我希望它是dm1 = 1 . 当我在ttest上进行常规计算({Beta(dm1)-1} / SE(Beta(dm1)))时,我得到新的t统计量应该在-48.89左右 . 如果这不是正确的方法,那么确定系数在统计上是否与... -
1 votesanswersviews
R多元回归循环和提取系数
我必须对同一自变量矩阵的许多因变量向量进行多元线性回归 . 例如,我想创建3个模型,以便: lm( d ~ a + b + c ) lm( e ~ a + b + c ) lm( f ~ a + b + c ) 来自以下矩阵(a,b,c是自变量,d,e,f是因变量) [,1] [,2] [,3] [,4] [,5] [,6] [1,] a1 ... -
0 votesanswersviews
如何在matlab中手动编写相关系数?
以下是采用两个相等大小的矢量X和Y的函数,并且假设返回包含用于图像对应的单个相关系数的矢量 . 如果给出两个相等大小的向量,该函数应该与matlab中的内置corr(X,Y)函数类似地工作 . 现在我的代码正在生成一个包含多个双数向量的向量,而不是包含单个数字的向量 . 我该如何解决? function result = myCorr(X, Y) meanX = mean(X); meanY =... -
-1 votesanswersviews
获得以下相关系数是否合乎逻辑且有意义?
在给定的表中,我们有从独立变量到因变量的一对多关系?相关系数是否有意义?如果所有值都是12,我们可以推断出什么? -
0 votesanswersviews
当回归R中存在一个预测变量时,提取预测变量名称
所以我目前正在运行逻辑回归,并试图提取各个预测变量的摘要信息,而不包括如下截距: zscores1<-summary(step1)$coefficients[-1,"z value"] > zscores1 i3 i2 i1 0.5011802 2.7834229 2.0239975 Step1是这里感兴趣的模型,用于提供上下文... -
2 votesanswersviews
最小二乘回归系数非线性函数的标准误差和置信区间
我在R中运行OLS回归,从中得到几个系数 . 这是代码的一部分: Attacks <- Treat.Terr.Dataset$Attacks[2:30] Attackslag <- Treat.Terr.Dataset$Attacks[1:29] TreatmentEffect <- Treat.Terr.Dataset$TreatmentEffect[2:30] Treatme... -
0 votesanswersviews
提取R插入符号内的系数
library(caret) data(iris) train_control <- trainControl(method="repeatedcv", number=10, repeats=10) model <- train(Sepal.Length~Sepal.Width+Petal.Length+Petal.Width, data=iris, trCont... -
0 votesanswersviews
使用前向回归和/或LARS从插入符号中的R的列车函数中检索系数
我在R工作,并使用几种方法探索使用插入符号进行变量选择和加权 . 在这里,我正在探索使用前向逐步和最小角度回归(LARS),使用每个的调整参数 . 在下面的代码中,我任意选择了一个因变量(y)和一个预测变量子集(x's),并使用70%的数据子集通过训练算法运行它们 . 为此,我正在应用重复的10倍交叉验证 . 我正在努力的是定位命令以识别从列车功能导出的最终模型参数(例如,截距,β权重) . 当我... -
0 votesanswersviews
在lmer中选择一个拦截
我正在使用lmer模型(来自lmerTest)来了解大小是否与基因表达显着相关,如果是这样,哪些特定基因与大小相关(也将'女性'和'笼'视为随机效应): lmer(Expression ~ size*genes + (1|female) + (1|cage), data = df) 在摘要输出中,我的一个基因被用作截距(因为它在字母表中最高,'ctsk') . 在阅读之后,建议我选择最高(或最低... -
0 votesanswersviews
线性回归不返回所有系数
我正在运行所有预测变量的线性回归(我有384个预测变量),但只从汇总中得到373个系数 . 我想知道为什么R不返回所有系数,我怎么能得到所有384个系数? full_lm <- lm(Y ~ ., data=dat[,2:385]) #384 predictors coef_lm <- as.matrix(summary(full_lm)$coefficients[,4]) #only...