R中混合效果模型的随机效果的标准误差(BLUP)
R中混合效果模型的随机效果的标准误差（BLUP）是指在混合效果模型中，对于随机效果的估计值的标准误差。BLUP代表最佳线性无偏预测，它是一种通过最小化预测值与真实值之间的平方差来估计随机效应的方法。
混合效果模型是一种统计模型，用于分析具有随机效应和固定效应的数据。随机效应是指在数据中存在随机变化的因素，而固定效应是指在数据中存在固定的因素。混合效果模型可以用于处理各种类型的数据，包括实验设计、观测数据和面板数据等。
BLUP在混合效果模型中的应用非常广泛。它可以用于估计随机效应的方差和协方差，从而提供关于随机效应的不确定性的信息。BLUP还可以用于预测随机效应的值，并计算其标准误差，从而评估预测的准确性。
在R中，可以使用lme4包来拟合混合效果模型，并使用ranef函数来计算随机效应的BLUP。通过调用summary函数，可以获取随机效应的标准误差。
以下是腾讯云提供的与混合效果模型相关的产品和产品介绍链接地址：
- 腾讯云机器学习平台（https://cloud.tencent.com/product/tcmlp）：提供了丰富的机器学习工具和算法，可用于构建混合效果模型和进行相关分析。
- 腾讯云数据分析平台（https://cloud.tencent.com/product/dp）：提供了数据处理和分析的工具，可用于预处理和分析混合效果模型所需的数据。
- 腾讯云大数据平台（https://cloud.tencent.com/product/emr）：提供了大数据处理和分析的解决方案，可用于处理大规模数据集和进行复杂的混合效果模型分析。
请注意，以上产品和链接仅为示例，可能并非完全适用于混合效果模型的随机效果的标准误差（BLUP）的应用。在实际应用中，建议根据具体需求选择适合的工具和平台。
相关·内容
在元分析包元中的R中，rma.uni()有一个blup()函数，但没有rma.mv()。我有一个带有随机截距项的rma.mv模型。如果我用：
predict.rma(model,transf=exp)
然后，我从固定效果版主那里得到了对原始数据集中的每个数据点的估计值，并且：
ranef(model,transf=exp)
给我预测我的随机效应的每一个层次。
那么，是否有办法将这两个函数的信息组合成随机和固定效果的组合blup，就像blup()函数所提供的那样？
(我试着对原始数据集中的每个数据点进行固定效果估计和随机效应预测，当我绘制它时，这看起来是正确的。)但这肯定不是那么简单吗？)
浏览 6提问于2021-03-14得票数 1
回答已采纳
1回答
因此，当我认为某一特定特性组存在依赖时，就会使用混合效应回归模型。我附加了Wiki链接，因为它比我解释得更好。()
虽然我认为在很多情况下我们需要考虑混合效果，但支持这种效果的模块不多。R有lme4和Python似乎有类似的模块，但它们都是统计驱动的；它们不使用梯度增强之类的代价函数算法。
在机器学习设置中，如何处理需要考虑混合效果的情况？有没有其他模型可以处理具有混合效应(随机效应)的纵向数据？
(R似乎有一个支持混合效果的包：，但我正在寻找一种Python解决方案。
浏览 1提问于2017-09-08得票数 2
我的混合模型公式在SAS proc mixed语法中的随机效果如下所示：
random intercept color size
/type = vc subject = group solution;
我将其转换为R lmer syntax，如下所示：
((1|group) + (0 + color|group)) + ((1|group) + (0 + size|group))
这是正确的吗？
我可以将sas随机效果公式表示为：
(1|group) + (0 + color|group)) + (0 + size|group) ?
或者它是R中的错误实现
浏览 0提问于2020-11-22得票数 0
我希望将summary()的固定效果部分提取为data.frame。我使用lme4运行以下模型：
SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)
summary(SleepStudy)
我知道我可以使用以下方法提取summary的随机效应部分：
SleepStudy_RE <- as.data.frame(VarCorr(SleepStudy))
对于固定效果是否有类似的代码行，包括估计值、标准误差、自由度和精确p值？
谢谢。
浏览 3提问于2016-05-20得票数 5
回答已采纳
因此，我在R中拟合了两个随机拦截的线性混合模型：
Y = X beta + Z b + e_i,
其中b ~ MVN (0, Sigma)；X和Z分别是固定和随机效应模型矩阵，beta和b是固定效应参数和随机效应BLUP/条件模式。
我想了解一下b的基本协方差矩阵，这在lme4包中似乎不是一件小事。您只能通过VarCorr获得方差，而不是实际的相关矩阵。
根据 (第2页)，您可以计算beta：e_i * lambda * t(lambda)的协方差。所有这些组件都可以从lme4的输出中提取。
我在想这是不是该走的路？或者你还有其他建议吗？
浏览 2提问于2016-06-13得票数 6
回答已采纳
7回答
在lm和glm模型中，我使用函数coef和confint来实现此目标：
m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous
coef(m)
confint(m)
现在，我使用lme4包中的lmer函数将随机效果添加到模型使用的混合效果模型中。但是，函数coef和confint对我不再起作用了！
> mix1 = lmer(resp ~ 0 + var1 + var1:var2 + (1|var3))
# var1, va
浏览 4提问于2012-06-17得票数 33
回答已采纳
我在我的模型中发现了一些异方差，我想用更健壮的标准误差来补偿。我尝试在R中使用merDeriv包中的Huber-White标准误差，但我认为这些只适用于具有二项分布的GLMM。对于负二项分布，有没有办法达到同样的效果？ 型号： library(lme4)
model <- glmer.nb(Jobs ~ 1 + Month + Year + (1|Region), data = df) Huber-White鲁棒标准误差： library(merDeriv)
bread.glmerMod(model) 错误： Error in vcov.lmerMod(object, full = f
浏览 76提问于2020-06-28得票数 3
回答已采纳
1回答
我正在尝试做一个在r或stata中具有随机效果的栅栏模型。我已经查看了glmmADMB包，但在R中下载它时遇到了问题，并且我在Cran中找不到任何关于该包的文档。这个套餐还能用吗？有没有人用它成功地估计了一个具有随机效应的栅栏模型？
或者，有没有一种方法可以在stata中估计这个值？在stata中，有没有一种方法可以用任何类型的计数数据来估计随机效应？
任何建议都将不胜感激。
詹妮弗
浏览 5提问于2013-05-25得票数 2
我正在处理具有随机效果的多变量数据。
我的假设是：D对A1和A2有影响，其中A1和A2是二进制数据，D是连续变量。
我还有一个随机效应，R，它是一个因子变量。
所以我的模型应该是这样的：A1andA2~D, random=1=~1|R
我尝试使用mvabund包中的函数manyglm，但它不能处理随机效果。或者我可以使用lme4，但它不能处理多变量数据。
我可以将我的多变量数据转换为4级因子变量，但我没有找到任何方法来使用不是二进制而是因子数据作为响应变量。我还可以将连续的D转换为因子变量。
对于在这种情况下使用什么，你有什么建议吗？
浏览 0提问于2018-06-30得票数 1
2回答
在统计中，如果随机变量的所有随机变量都有相同的有限方差，那么随机变量序列是同方差的。这也称为方差的同质性(维基百科)。
这个假设是否仅适用于线性回归模型？
如果没有，是否有任何相关的方法来检查非线性模型？
浏览 0提问于2021-04-12得票数 2
回答已采纳
我正在使用Proc glimmix对五个时间点上的连续结果运行随机截取模型。表类型III表示性别不显着，而表lsmeans中95%的男性和女性CI没有重叠。我想知道如何解释这种差异？当我观察男性和女性随着时间的推移的结果时，我也预计男性的结果会显着更高，就像第二个表中的顺位所说的那样。非常感谢！ ?
浏览 48提问于2021-03-30得票数 0
回答已采纳
我想要估计一个固定效应模型，同时使用面板校正的标准误差以及普拉斯-温斯顿(AR1)变换，以解决面板异方差，同时空间相关性和自相关性。
我有时间序列横截面数据，想要执行回归分析。我能够分别估计固定效应模型、面板校正标准误差和Prais-winsten估计。我能够将面板校正后的标准误差包含在固定效果模型中。但我想一次把它们都拿出来。
# Basic ols model
ols1 <- lm(y ~ x1 + x2, data = data)
summary(ols1)
# Fixed effects model
library('plm')
plm1 <- plm(y
浏览 9提问于2019-05-08得票数 1
1回答
我有一个数据集，其中包含有关天气、空气污染和健康结果的信息。我想回归温度(T)和温度滞后(T1)与心源性死亡(CVD)的关系。我之前在R中使用过glm模型，使用的脚本如下： #for mean daily temperature and temperature lags separately.
modelT<-glm(cvd~T, data=datapoisson, family=poisson(link="log"), na=na.omit) 我得到了效果估计值和标准误差值，我用它们来换算成风险比率。 现在我想使用动态线性模型或分布式线性模型来检查预测结果和滞后预测结
浏览 22提问于2020-07-22得票数 1
我试着尽可能多地向你们展示数据的结构和我所产生的结果。
数据的结构如下：
GroupID Person Factor2 Factor1 Rating
<int> <int> <fctr> <fctr> <int>
1 2 109 2 0 1
2 2 109 2 1 -2
3 2 104 1 0 4
4 2 236 1 1
浏览 0提问于2017-11-06得票数 1
回答已采纳
我正在尝试将NA传递给R函数，例如，使用仅使用固定效果(即没有随机效果)的lme4混合模型进行预测：
import rpy2.rinterface as ri
from rpy2.robjects.packages import importr
rstats = importr('stats')
rstats.predict( mymodel, re_form=ri.NA_Logical )
但是，由于某些原因，re_form=ri.NA_Logical无法将NA传递给re.form (我还尝试使用别名REform、ReForm等)。有什么想法吗？
此R函数：
浏览 0提问于2019-02-24得票数 1
我有一个数据集，数据被保留审查，我想应用多级混合效应tobit回归，但我只找到信息，如何在Stata。在R中可以这样做吗？
我找到了'VGAM‘和'CensREG’包，但我不知道如何添加固定的和随机的效果。
另外，我的数据是日志正态分布的，有没有办法将这个添加到模型中？
谢谢!
浏览 2提问于2022-09-29得票数 0
回答已采纳
我想知道如何写出两个相互作用因素的随机效应。例如，我有6个种，在48个地块中种植，并在两个块中复制。共有48个组合(所有可能的1，2，3，5和6种组合)的6个物种。
通过设计，组合实际上是在多样性内嵌套的:组合是每个多样性水平的正确复制，因此多样性的效果应该根据多样性levels.Since中组合之间的随机变化来检验，多样性或物种丰富度是一个连续变量，我不能认为组合是嵌套在多样性内的。
我想把块和组合之间的相互作用作为随机效应。我不知道如何在R.中混合模型中正确地写
浏览 6提问于2022-03-03得票数 0
1回答
如何修剪R中的树？
、、、
我正在使用R中的rpart进行分类。树模型的训练方法是：
> tree <- rpart(activity ~ . , data=trainData)
> pData1 <- predict(tree, testData, type="class")
此树模型的精度为：
> sum(testData$activity==pData1)/length(pData1)
[1] 0.8094276
我读了一个通过交叉验证修剪树的：
> ptree <- prune(tree,cp=tree$cptable[which.min(tree$cpt
浏览 1提问于2013-03-10得票数 13
回答已采纳
1回答
我使用的plm函数(来自plm包)使用的是固定效果。
我需要保证在使用函数plm时，我使用的是异方差健壮的标准误差。
lag1.1 <- plm(E.primary ~ lv18_bank_c + Expense + lag(Expense, 1), data=panel_exog1, index=c("Country", "Year"), model="within")
如何获得该模型的异方差鲁棒标准误差？
浏览 6提问于2020-08-26得票数 1
1回答
我想请一位同事复制一个第一差分线性面板数据模型，该模型是我用R中的plm软件包(或其他软件包)用Stata估计的。
在Stata中，xtreg没有第一个差异选项，因此我运行：
reg D.(y x), nocons cluster(ID)
在R中，我正在做：
plm(formula = y ~ -1 + x, data = data, model = "fd", index = c("ID","Period"))
系数匹配，但R中的标准误差大于Stata中的标准误差。我查看了plm帮助和pdf文档，但我肯定遗漏了什么。
浏览 0提问于2013-09-26得票数 10
回答已采纳
2回答
我有一个lme4模型，我运行了一个分层逻辑回归模型，并且我正在使用effects包绘制效果。我想用平均值的标准误差作为误差条来创建一个效果图。我可以得到点估计，95%的置信区间，和标准错误的数据。然而，标准误差似乎与置信限参数不一致，参见下面的一个普通glm示例。
library(effects)
library(dplyr)
mtcars <- mtcars %>%
mutate(vs = factor(vs))
glm1 <- glm(am ~ vs, mtcars, family = "binomial")
(glm1_eff <- Effect(
浏览 1提问于2016-04-17得票数 2
回答已采纳
1回答
最近，我进行了一项研究，将参与者分为两种情况。他们给不同的物品打分。我想测试两种情况下他们的收视率是否不同。由于这些项目有很大的不同，而且我对每个参与者进行了许多度量，所以我决定为项目(编码为"item")和参与者(编码为"id")添加至少随机的拦截。
我将以下模型与lme4包相匹配：
lmer(y ~ condition + (1|id) + (1|item), dfLong)
但是，我将以下简单的线性模型与通过平均每个参与者的所有响应来创建的指数度量相匹配：
lm(itemMean ~ condition, dfWide)
结果是和我在LMM中得到的完全相同
浏览 7提问于2022-03-20得票数 4
1回答
我有一个固定的效果模型，只有很少的观测值，因此我想要引导以获得更准确的标准误差。同时，我假设SE是集群的，因此我也想纠正集群，即执行集群引导。我找到了一个用于lm模型的函数(vcovBS)，但是找不到任何用于plm模型的函数。有没有人知道一个类似的函数来为固定效果模型获得集群自举SE？
浏览 12提问于2021-03-31得票数 0
2回答
lm和lme结果的标准误差
、、、、
我有以下线性模型
library(nlme)
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)
fm2.lm <- lm(distance ~ age + Sex,data = Orthodont)
如何获得距离随年龄和性别的标准误差？
浏览 2提问于2014-03-10得票数 1
回答已采纳
1回答
随机效应Cox模型
、、、
我正在尝试绘制pm2.5暴露与高血压发病率之间的剂量-反应关系。我拟合了一个随机效应cox模型，其中为研究中心添加了一个随机效应。然后使用plotHR函数绘制剂量-反应关系图。但我遇到了一个错误。下面是我使用的示例R代码。
library(survival)
library(coxme)
library(splines)
library(Greg)
data("eortc")
eortc$age<-rnorm(2323,40,10)
efit1 <- coxph(Surv(y, uncens) ~ ns(age) + trt , eortc)
plotHR(e
浏览 13提问于2018-09-05得票数 1
回答已采纳
我对R很陌生，我正试图计算线性模型的95%置信区间，用引导法重新采样响应变量，然后在原来的解释变量上对这999个自举响应变量进行回归，从而形成线性模型的剩余标准误差。
首先，我不确定是否应该计算原始线性模型的95% CI和残差标准差(没有引导数据)，因为这是没有意义的--对于这个线性模型，R平方值是100%精确的，计算它的CI是没有意义的。
是对的吗？
重要的是，我不知道如何为我通过引导创建的999个线性模型的R平方值和残差标准误差值计算CI。
浏览 3提问于2013-09-25得票数 1
这可能是一个显而易见的问题，我还没有找到(R newbie)的答案，但是当使用lmer函数生成混合效果模型，然后使用以下命令显示结果时：
screenreg(list(model4re), single.row = TRUE)
我们得到了以恒星形式表示的beta估计值、标准误差和显著性水平的列表。
使用什么测试来确定这些p值来标记恒星(重要的是，我认识到关于如何使用这些模型适当地确定显着影响存在一些争议)，以及我们如何提取用于这些恒星的p值？
浏览 1提问于2017-05-12得票数 1
1回答
生存包在R中的鲁棒性标准误差
、、、、
我试图从R.中的生存包中获得一个严重的标准错误，在这样做的时候，我尝试用clogit选项复制Stata vce(robust)命令报告的标准错误。
我的公式R是
conditional_logit <- clogit(dependent_variable ~ independent_variable + some_controls + strata(year), method= "exact", data = data_frame)
向函数添加robust = TRUE参数将失败，错误如下：
Error in residuals.coxph(fit2, type = &#
浏览 4提问于2016-08-10得票数 4
回答已采纳
1回答
通过sklearn来发现ML世界，我正在我的数据集上测试一个大的模型面板。这是为了学习的目的，也是为了工作，所以我希望最终的模型尽可能准确，而我可以在我的理解ML。
我已经将我的数据集(16k行)划分为80%的训练和20%的测试，我正在测试至少KNN、Logistic、DecisionTree、RandomForest、NaivesBayes和SVC (如果我的计算机能处理的话)，当我知道如何处理时，用Bagging和Boosting测试。
我的培训示例有4种方式:2组特性(95和11)，标准的(带有StandardScaler的)与否。
我的结果是二进制的，我使用了一个自定义的计分器“改进”
浏览 0提问于2019-05-03得票数 3
回答已采纳
即使状态变为"UP"，而待定状态仍将继续。为什么我不能在循环中使用变量作为条件？ export confup=false
export blup=false
timeout 30 bash <<EOF || false
while [[ $( curl http://0.0.0.0/ | jq -r '.status') != "UP" ]];
do
if [[ $( curl http://0.0.0.0/ | jq -r '.components."bridge.lock".status') =
浏览 42提问于2020-09-01得票数 0
回答已采纳
proc glm使添加固定效果变得很容易，而无需为类变量的每个可能值创建虚拟变量。
proc reg能够计算健壮(白色)标准误差，但它要求您创建单独的虚拟变量。
有没有办法把这些功能结合起来？我希望能够添加一些类变量，并在我的输出中接收白色标准错误。例如：
使用proc glm，我可以进行这种回归。这将给出正确的结果，无论类变量中包含多少级别，但它不会计算健壮的标准误差。
proc glm data=ds1;
class class1 class2 class3;
weight n;
model y = c class1 class2 class3 / solution;
run;
浏览 3提问于2014-05-09得票数 1
1回答
我试图用一个包含随机效应的模型来预测一个二元结果，这个模型使用了调查数据。我已经在下面描述了抽样设计，所以请随意评论我的调查权重方法。我的主要问题是如何在调查加权模型中包含随机效应。以下是到目前为止的代码：
# Libraries
library(survey)
# Make dataframe object where d is the working dataframe
dfobj <- svydesign((id = ~cluster+household, strata = ~interaction(region, urban) weights = ~chweight, stra
浏览 5提问于2012-08-17得票数 4
3回答
我需要从standard error的输出中提取lmer的方差分量。
library(lme4)
model <- lmer(Reaction ~ Days + (1|Subject), sleepstudy)
以下是差异部分的估计数：
s2 <- VarCorr(model)$Subject[1]
这是，而不是，方差的标准误差。我要标准错误。我怎么能得到它？
编辑：
也许我无法让您理解我所说的“方差分量的标准错误”是什么意思。所以我在编辑我的帖子。
在第12章，随机因素的实验，由道格拉斯C.蒙哥马利(DouglasC.Montgomery)著的“”一书中，在本章末尾，示例12-2
浏览 2提问于2015-07-29得票数 3
回答已采纳
我试图在SAS中进行一个简单的单向随机效应方差分析。我想知道总体方差是否与零显著不同。
在上，它们声明使用PROC混合如下：
proc mixed data = in.hsb12 covtest noclprint;
class school;
model mathach = / solution;
random intercept / subject = school;
run;
考虑到我以前使用PROC混合的经验，这对我来说是有意义的。
然而，在莫里·洛根的“使用R的生物统计设计和分析”一书中，他说，对于单向的方差，固定的和随机的影响是没有区别的，并且在R中进行了“标准”
浏览 3提问于2014-10-05得票数 2
回答已采纳
1回答
我正在使用广义线性混合模型(glmer，lme4-package)使用R分析我的二项式数据集。我想使用Tukey's post-hoc测试(glht，multcomp package)对某个固定效果("Sound")进行成对比较。
它大部分都工作得很好，但我的一个固定效果变量("SoundC")根本没有变化(96倍于"1“，0倍于"0")，似乎Tukey的测试无法处理这一点。所有与这个"SoundC“的成对比较都会给出1.000的p值，而有些显然是有意义的。
作为验证，我将96个“1”中的一个更改为"0“，之后我又
浏览 0提问于2015-06-24得票数 8
回答已采纳
我有一个具有二元结果的回归模型。我用glmnet对模型进行了拟合，得到了选定的变量及其系数。
由于glmnet不计算变量重要性，我想将确切的输出(选定的变量及其系数)提供给glm以获取信息(标准误差等)。
我搜索了r个文档，似乎可以使用glm中的"method“选项来指定用户定义的函数。但是我没有做到这一点，有人能帮我吗？
浏览 2提问于2012-10-17得票数 21
回答已采纳
4回答
R中的大固定效应二项回归
、、、、
我需要在一个相对较大的数据框架上运行一个逻辑回归，其中有480.000个条目，有3个固定的效果变量。固定效应变量A为3233级，B级为2326级，C级为811级。总之，我有6370个固定效果。数据是横截面的。如果我无法使用普通的glm函数运行这个回归，因为回归矩阵对我的内存来说太大了(我得到了消息"Error: cannot allocate vector of size 22.9 Gb")。我在我的Macbook Air (OSX10.9.58GBRAM)上寻找其他方法来运行这个回归。我还可以访问一个拥有16 to的服务器。
我曾尝试以几种不同的方式解决这个问题，但迄今都没有
浏览 10提问于2015-02-20得票数 18
我们能否从R中的混合模型计算BLUP的标准误差(SE)？我正在使用一个名为ASReml的自定义包，它可以根据预测值计算SE。然而，我不确定在哪里可以找到BLUP值的SE。
到目前为止，我有这样的东西...
df<-data.frame(inbred=c("x1","x2","x3","x4","x5","x1","x2","x3","x4","x5","x1","x2","x3",&
浏览 6提问于2020-02-27得票数 2
我正在学习线性回归，我运行了一个线性回归的阶跃函数，并检验了最终方程的残差和拟合图。残差看起来是同步性的，但它不是随机分布在线的上面和下面。我不知道该怎么解释。残差与杠杆图并不表示我有一个高厨师的距离点。我的R平方值非常糟糕.以下是总结输出：
剩余标准误差: 58588自由度上的2.337
倍数R-平方: 0.3692，调整R-平方: 0.3671
F-统计量: 174.1对197和58588 DF，p-值：< 2.2e-16
但我不太明白残差和贴图意味着什么。我在开发一个更好的模型的下一步将会是什么。
📷
浏览 0提问于2019-02-17得票数 3
1回答
国家-年对上的聚类鲁棒标准误差
、、、、
我想在R中复制Stata do.file (面板模型)，但不幸的是，我最终得到了错误的标准误差估计。这些数据是专有的，所以我不能在这里发布。使用的Stata代码如下所示：
xtreg Y X, vce(cluster countrycodeid) fe nonest dfadj
使用fe表示固定效果，使用nonest表示面板不嵌套在集群中，使用dfadj表示发生了某种DF调整-到目前为止不可能找出是哪种排序。
我的R-Code看起来像这样，并使我最终得到正确的系数值：
model <- plm(Y~X+as.factor(year),data=panel,model="wit
浏览 2提问于2015-11-30得票数 0
1回答
我有几个使用plm和pooling的回归模型。我的数据是一个集合的横截面/时间序列数据，与债券发行的数据。这些数据包括对大约2000个债券发行情况的观察，约有25个债券描述变量.
我想要计算一个或所有回归模型的稳健标准误差，以便将它添加到我的天文望远镜可视化中。我的回归如下：
#Regression
primaryreg4 <- plm(issueyield ~ issuer + exchange + yearmonth + maturity.cat + size.cat + coupontype,
data = data,
浏览 5提问于2019-11-20得票数 1
1回答
R中黄土的标准误差
、、、
我试图找到一个参考资料，解释如何计算局部多项式回归的标准误差？具体地说，在R中，可以使用loess函数来获取模型对象，然后使用预测函数来检索标准误差。有没有对实际发生的事情的引用？在残差中可能存在序列相关的情况下，必须使用Newey-West类型方法来调整它，有没有一种方法可以使用三明治包来实现这一点，就像您使用lm对常规OLS所做的那样？
我尝试查看源代码，但标准误差计算调用了一个C函数。
浏览 1提问于2012-10-03得票数 5
回答已采纳
我正在使用Spark ml库来做一些生存分析。是文档。在训练AFT生存模型后，我不能像在R中那样直接获得p值。可用于该模型的是系数向量、比例、截距、预测和输入数据(Y和特征)。如何使用这些变量来计算p值？欢迎任何建议。谢谢
浏览 3提问于2016-05-28得票数 0
当我使用R运行lm()回归时，我从summary()得到了“残差标准误差”。为什么只有一个残差标准误差值，而不是每个观测值的残差标准误差列表？
summary()中显示的此值的含义是什么？摘要()中显示的“残差标准误差”是每个观察值的残差标准误差列表的平均值吗？谢谢。
Residual standard error: 0.8498 on 44848 degrees of freedom
(7940 observations deleted due to missingness)
Multiple R-squared: 0.4377, Adjusted R-squared: 0.4
浏览 2提问于2019-07-31得票数 2
1回答
对于具有固定效果和随机效果的模型，有没有与R的plm包等效的Python或PySpark？
浏览 14提问于2020-05-27得票数 0
我想要计算预测值和标准误差，但我不能简单地使用predict()，因为我使用了15个倍数推算的数据集(生成的Amelia包)。我在每个数据集上运行回归模型。然后，使用Amelia函数mi.meld()将结果合并到一组模型系数和标准误差中，该函数使用鲁宾规则。
示例数据和代码：
dd<-list()
for (i in 1:15){
dd[[i]] <- data.frame(
Age=runif(50,20,90),
Cat=factor(sample(0:4, 50, replace=T)),
Outcome = sample(0:1, 50, replace=T)
浏览 3提问于2017-07-13得票数 3
1回答
我在R中运行一个混合模型，但是我在理解我应该为我拥有的数据运行的模型类型时遇到了一些困难。
让我们将因变量称为计算机化实验中早期按下按钮的次数。一个实验是由多个实验组成的。在每个试验中，参与者必须按下一个按钮来对屏幕上出现的目标做出反应。然而，他们可能过早地按下按钮，这是作为结果变量来衡量的。因此，例如，参与者A可能在整个试验中总共有3次提前按钮按下，而参与者B可能有15次。
在使用R中的lm命令的简单线性回归模型中，我认为这个结果是一个连续的数值变量。还有..。这是参与者在实验中打分的数字。然而，我不是试图运行线性回归，我是试图运行具有随机效果的混合模型。我对R中混合模型的理解是，模型采用的
浏览 0提问于2020-06-19得票数 1
如何在C++中获得给定均值和标准误差的随机数？有什么内置的功能可以直接使用吗？
我知道一些例子，它们用均值和标准差来得到一个随机数。但是，我不知道复制的次数，只有平均值和标准误差，所以我无法计算标准偏差。有什么想法吗？
我感谢你的时间和帮助。谢谢。:-)
浏览 1提问于2016-04-20得票数 0
回答已采纳
2回答
我有一个关于使用lmer在R中正确设置模型的问题。 这是一个重复测量的实验。 每个受试者(总共20人)在每种刺激条件(阳极、阴极和假刺激)下完成4个不同的任务。因变量为反应时间(rt) 我使用了这个模型，但我不确定它是否正确。我更关心随机效果是否被正确分配 model<- lmer(rt ~ task * stimulation + (task * stimulation|subject), data=dat) 任何帮助都将不胜感激。 谢谢
浏览 21提问于2019-03-12得票数 2
我正在尝试使用Python statsModels线性混合效果模型来拟合具有两个随机截获的模型，例如两组。我不知道如何初始化模型，这样才能做到这一点。
下面是示例。我有如下数据(取自)：
subject gender scenario attitude frequency
F1 F 1 pol 213.3
F1 F 1 inf 204.5
F1 F 2 pol 285.1
F1 F 2 inf 259.7
F1 F 3 pol 203.9
F1 F 3 inf 286.9
F1 F 4 pol 250.8
F1
浏览 1提问于2016-08-26得票数 5
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
