如何在观星计划中报告标准化模型?
我制作了以下简单的回归模型,并使用Stargazer输出绘制标准化与非标准化回归模型系数,标准误差和P值的表。 library(lm.beta) mod <- lm(mpg ~ cyl + …
r晶格面板功能拟合,如何获得每个面板的模型摘要
如何将摘要的打印输出以适合每个晶格面板的模型? 例如,使用LME4的SleepStudy Dataste: library(lme4) data(sleepstudy) 我可以将线性模型拟合到所…
为什么我会在predive.lm中遇到错误。
Closed. This question is not reproducible or was caused by typos. It is not currently accepting answers. 这个问题是由类错字或无法再现的问题…
如何从不同模型的方差分析输出上运行AIC
这是我使用ANOVA的一小部分数据。 dput(c) structure(list(Mouse.ID = c("DO.1307", "DO.1309", "DO.1311", "DO.1312", "DO.1316", "DO.1318", "DO.06…
在多列条件下运行ANOVA
这是我的数据框子集的样子。 a <- dput(head(mrna.pcs)) structure(list(Mouse.ID = c("DO.0661", "DO.0669", "DO.0670", "DO.0673", "DO.0674", "DO.…
将 R 线性模型重写为 Python
帮助将 R 线性模型重写为 Python。 R 代码: x <- rnorm(10) y <- 1+x+rnorm(10) model <- lm(y~x) res = summary(model)$r.squared print(res) Pytho…
数据帧的引导函数 - 将函数作为 R 中的参数传递
我正在尝试为我的作业创建一个引导功能。要求如下: 计算以下内容的引导程序标准错误: - meand()和 - 中位数()和 - 顶部四分位数和 - max() 价…
如何与个体级变量和个体固定效应绘制相互作用的影响?
我有一个数据集,其中包含单独的id(因子)、时间t(因子)、因变量y(连续)和自变量x(连续),可以在时间 t xt 测量,也可以在单独的级别 xi 上设…
我正在尝试使用 Predict.lm() 基于回归来预测新值,但结果显示新值的数量错误
我正在使用一个带有几个天气站的温度数据的数据框架,以给定的时间间隔。 Here are the 10 first observations data hora a001 a046 241 2021-03-20 0…
R nloptr 重建线性回归问题:优化过早停止
我尝试使用 nloptr 重建多元线性回归,但优化停止得太早。这是为什么呢? (我知道 nloptr 适用于非线性模型,但我仍然期望正确的解决方案,例如 opti…
如何处理 lm 中以多级因子为控制的矢量大小过大的情况
我正在尝试拟合一个包含大约 900,000 个观测值和两个解释变量的线性模型。然而,我还需要包含一个控制变量,它是一个多级因子变量(11,135 个级别)。…