固定效应联合显着性F检验统计量的提取
固定效应联合显着性F检验统计量是一种用于统计分析中的假设检验方法，用于评估固定效应模型中的自变量对因变量的影响是否显著。它是通过计算F统计量来进行判断的。
在固定效应模型中，自变量被认为是固定的，即其取值是确定的，而不是随机的。固定效应联合显著性F检验统计量的提取过程如下：
- 首先，建立固定效应模型。固定效应模型是一种多元线性回归模型，其中自变量是固定的，因变量是随机的。该模型可以表示为：Y = Xβ + ε，其中Y表示因变量，X表示自变量矩阵，β表示自变量的系数，ε表示误差项。
- 然后，进行方差分析。方差分析是一种用于比较组间差异的统计方法，它将总体方差分解为组间方差和组内方差。在固定效应模型中，方差分析用于评估自变量对因变量的整体影响是否显著。
- 接下来，计算F统计量。F统计量是组间方差与组内方差的比值，用于判断组间差异是否显著。计算公式为：F = (组间平方和 / 组间自由度) / (组内平方和 / 组内自由度)。
- 最后，进行假设检验。根据F统计量的计算结果，与给定的显著性水平进行比较，来判断自变量对因变量的影响是否显著。如果F统计量大于临界值，则拒绝原假设，认为自变量对因变量的影响是显著的。
固定效应联合显著性F检验统计量的提取在实际应用中具有广泛的应用场景，例如社会科学研究、经济学分析、医学研究等领域。在云计算领域中，该统计方法可以用于评估不同自变量对云计算性能、资源利用率等指标的影响。
腾讯云提供了一系列与云计算相关的产品，例如云服务器、云数据库、云存储等，可以帮助用户构建和管理云计算环境。具体产品介绍和链接地址如下：
- 云服务器（ECS）：提供弹性计算能力，支持多种操作系统和应用场景。了解更多：https://cloud.tencent.com/product/cvm
- 云数据库（CDB）：提供高可用、可扩展的数据库服务，支持多种数据库引擎。了解更多：https://cloud.tencent.com/product/cdb
- 云存储（COS）：提供安全可靠的对象存储服务，适用于存储和管理各种类型的数据。了解更多：https://cloud.tencent.com/product/cos
腾讯云的这些产品可以帮助用户构建稳定、高效的云计算环境，并提供丰富的功能和服务来满足不同的需求。
相关·内容
我正在运行以下代码：
library(lme4)
library(nlme)
nest.reg2 <- glmer(SS ~ (bd|cond), family = "binomial",
data = combined2)
coef(nest.reg2)
summary(nest.reg2)
它产生以下输出：
系数
$cond
bd (Intercept)
LL -1.014698 1.286768
no -3.053920 4.486349
SS -5.300883 8.011879
摘
浏览 1提问于2016-12-30得票数 1
1回答
三组假设检验
、、、、
根据这些数据，A、B和C组的平均销售额在统计上是否相同？
我在AB，BC，CA上做了t.test。对于CA，p值>0.05，所以我得出结论，对于CA，我们不能拒绝零假设，而且平均值可能是相同的。
H1-替代假说-- 36-45组和46-50组之间的均值差异并不等于0。
我的问题是-我做得对吗?还是有其他方法来检验三组的假设？
浏览 15提问于2022-12-03得票数 0
我是统计学新手。我正在试验一个样本Kolmogorov-Smirnov测试。我能在D max之前找到。但我很困惑于假设检验的进展。
为了确定这个假设，我是否应该--
表中的临界值
如果测试统计量D大于从表中获得的临界值，则拒绝
或
KS状态值P值。
哪一个更好?读那个p值更好。在中，他们告诉"kstest决定通过将p值p与显着性水平Alpha进行比较，而不是通过将测试统计量ksstat与临界值cv进行比较，从而拒绝空假设。由于kstest是近似的，因此将kstest与cv进行比较，有时会得出与p与Alpha不同的结论。“
但却找不到任何关于相同的等式。
在上面的参考文献中，我怀疑他们
浏览 3提问于2015-10-08得票数 0
我不明白为什么lmerTest中的处理步骤(模型)的结果是异常的。
m0 <- lmer(seed ~ connection*age + (1|unit), data = test)
step(m0)
注意：“连接”和“年龄”都设置为as.factor()
Random effects:
Chi.sq Chi.DF elim.num p.value
unit 0.25 1 1 0.6194
Fixed effects:
Analysis of Variance Table
Response: y
浏览 5提问于2014-11-23得票数 1
回答已采纳
我已经手工计算了测试数据，但想确认我的答案是正确的。我不知道从哪里开始，在给向量赋值之后。
在大样本情况下，Wald检验统计量、分数检验统计量和似然比检验统计量近似均具有chisquare df=1分布。Y~二项(1，π)，P(Y=1)=pi，P(Y=0)=1-pi.观察到，i=1 to m与Y相同。
对于m= 10和数据(0，1，0，0，1，0，0，1，0)，使用这些统计量检验空假设pi=0.3。
免责声明:虽然这是与作业有关的，但R部分不是作业的一部分.我很好奇我是否能在R中运行这个问题来确认我的答案。
浏览 2提问于2016-10-02得票数 1
回答已采纳
我正在尝试拟合几个混合线性模型。
summary(model3 <- lme4::lmer(LPP2POz ~ COND + (COND|ID), data = dataLPP2POz))
其输出是
误差:观测次数(=75)项的随机效应数(=75)；随机效应参数和剩余方差(或尺度参数)可能无法识别。
任何人都能弄清楚为什么这个模型会逆转这个错误吗？
以下是数据集：
> dput(head(dataLPP2POz))
structure(list(ID = structure(c(1L, 1L, 1L, 2L, 2L, 2L), .Label = c("01"
浏览 2提问于2021-10-27得票数 2
回答已采纳
我正在尝试对面板数据OLS回归(在R中)的固定效应(个人特定的虚拟变量)的联合显着性进行F检验，但是我还没有找到一种方法来实现大量固定效应。理想情况下，我会使用plm包中的函数，但是我还没有找到任何专门执行此测试的函数。
这是Stata在使用xtreg, fe命令时自动执行的操作。在Stata中，结果如下所示：
------------------------------------------------------------------------------
F test that all u_i=0: F(49, 498) = 12.00 Prob > F = 0.
浏览 7提问于2011-05-30得票数 3
回答已采纳
我正在处理面板数据集，并试图运行具有固定效果的logit回归。
我发现来自lme4包和bife包的glmer模型适合这种工作。
然而，当我对每个模型进行回归时，我没有相同的结果(估计、标准误差等)。
下面是带有截距的glmer模型的代码和结果：
glmer_1 <- glmer(CVC_dummy~at_log + (1|year), data=own, family=binomial(link="logit"))
summary(glmer_1)
Estimate Std. Error zvalue Pr(>|z|)
(I
浏览 1提问于2019-02-20得票数 3
我正在用图书馆重复测量R中的变异数：
library(ordinal)
library(car)
library(RVAideMemoire)
我有两组:月份和距离，因变量是CO2：
距离月CO2
0米7月234
我为CO2建立了一个clmm模型，解释了距离、月份以及月份和距离之间的交互作用：
model_CO2 = clmm(CO2.f ~ month + distance + month:distance + (1|nest),
data = field_data,
threshold = "equidistant")
浏览 2提问于2020-04-24得票数 1
回答已采纳
建立了具有固定效应和随机截距效应的广义加性混合效应模型(即范畴变量)。在运行模型之后，我能够使用ranef(m1$lme)$x[[1]]提取每个类别的随机截取。但是，当我试图使用se.ranef(m1$lme)提取随机效应的标准错误时，该函数不起作用。其他使用se.ranef(m1)和se.ranef(m1$gam)的尝试也不起作用。我不知道这是不是因为这些函数只适用于类lmer的模型？
有人能帮我从“游戏”类中提取出随机拦截的标准错误吗？我想使用随机截取和标准错误来绘制我的gamm模型的最佳线性无偏预测器。
我的初始模型是形式：gamm(y ~ s(z), random = list(x =
浏览 2提问于2018-01-19得票数 7
1回答
我有一个数据集plink.raw，我正在测试maker CN00020133 (有3级0、1、2)是否与phenotype5相关联。我想使用GLM或fisher提取测试来比较0 vs 1和2 vs1：
table(plink.raw$phenotype5,plink.raw$CN00020133)
1 0 2
0 3559 0 7
1 14806 54 123
使用GLM进行测试，我可以看到0和1的p值是0.912894。
plink.raw$CN00020133 <- f
浏览 2提问于2019-12-13得票数 1
2回答
嗨，伙计们，我想提取图像中最突出的部分。为此，我希望在java中使用显著检测器算法，并希望在位图图像上实现此算法。让我告诉你一些事情，我不能使用openCV来做这件事，所以请帮我用简单的算法将它们存档。这就是我想要的东西..
我先试过sobel滤波器，但没有得到这样的结果，顺便说一下，只讨论显着性部分，其他结果不应该被考虑。提前感谢
浏览 1提问于2013-05-08得票数 0
回答已采纳
嘿，有人知道如何遍历RTAQ包中的MedRV函数吗？我有过
##find days in data
ddx.f = endpoints(full, on="days");
days.full = format(index(full)[ddx.f], "%Y-%m-%d");
for (day in days.full) {
x = full[day]
}
提取每一天，但是不知道每天通过MedRV函数向循环中添加什么。有什么想法吗？还有人知道在R中有一个显着性测试来测试跳跃吗？
提前谢谢。
浏览 0提问于2012-06-16得票数 0
回答已采纳
我想知道是否有人知道从R.中glmmadmb的输出中提取色散参数估计值的方法，我使用负二项式模型，并且希望对几个不同的物种使用这个代码，而不必手动提取代码的其余部分的值。
这是我输出的一个例子：
> summary(mod1)
Call:
glmmadmb(formula = species ~ (1 | year) + (1 | site), data = cs,
family = "nbinom2", link = "log")
AIC: 8131.7
浏览 1提问于2013-10-02得票数 0
回答已采纳
dataset dt由16万6千行组成，列为: log_price:依赖的数字回归器sku: 381级的独立分类回归器；具有15级transaction_type的独立分类回归器；2级的独立分类回归器；购买者:具有1001级regressor_01、.、regressor_04:四个独立的数字回归体。
我们只考虑沿数据集显示超过200行的skus。
model_1回归花了几分钟时间，得到了很好的结果：
model_1<- lmer (formula = log_price ~
0 +
浏览 5提问于2020-11-23得票数 1
回答已采纳
1回答
我正试着用等价物来画我的模型。
c1<-lmer(sqrt(co22)~t.air+t.soil+tdr+ph+TC+TN+wtd+Fe+0+(1|treatment/site), data=data)
然而，当我试图提取方程时，我总是得到一个错误：
extract_eq(c1,wrap = TRUE, ital_vars = FALSE)
Error in `[<-`(`*tmp*`, v[1], v[2], value = greek_vcov) :
subscript out of bounds
但是，如果我将模型转换为一个没有随机效应的简单线性模型，效果会很好：
c
浏览 2提问于2021-06-14得票数 1
回答已采纳
3回答
我想控制一个包含超过100个级别的因素变量，而不将该控件的结果输出到汇总表中。注意，我还对复制Stata命令的速度感兴趣，而不仅仅是对输出进行表面上的更改。
在Stata中，我可以像这样使用“吸收”：
use http://www.stata-press.com/data/r14/abdata.dta, clear
. xtreg n w k i.year, fe
Fixed-effects (within) regression Number of obs = 1,031
Group variable: id
浏览 0提问于2018-04-03得票数 2
回答已采纳
1回答
我希望运行一个混合效应模型在R的基础上，我如何运行统计在SPSS与重复的测量方差。是我在SPSS中建立重复测量方差的方法。我将如何在R中将其转换为lme4？
Key：EBT100.是任务的名称，基因型是我的IV，我的学科内因素是天(5级)和Cue (9级)。是我的DV。
在R中，这是我试图运行的代码:在R中，这是我的代码：
lmeModel <- lmer(Att ~ Genotype*Day*Cue + (1|Subject)
我的基因型效应在R和SPSS之间是相同的(p~0.12)，但我的所有交互作用都是不同的(基因型x日，基因型x Cue，基因型x型x日x Cue)。
R (lme
浏览 0提问于2018-04-15得票数 0
1回答
我有一个小任务要做，不幸的是我不熟悉这个统计领域.事实上，我确实需要计算(我不是在寻找现成的解决方案)，但是我不知道它们是否正确，也不知道我的思维方式，因此如果你看一看并指出我的错误，我将非常感激。
这里是假数据，显示狗和猫的增长率(完全是虚构的)：
time <- c(1,2,3,4,5,6,7,8,9,10,1,2,3,4,5,6,7,8,9,10)
a <- rep('dog', 10)
b <- rep('cat', 10)
animal <- c(a,b)
val <- c(2.00,8.00,17.00,21.00,29
浏览 1提问于2018-08-07得票数 1
回答已采纳
我的任务是将35万个文档分为“已签名”和“未签名”两类。使用开源工具搜索看起来像人类签名的东西最快的方法是什么？为了使问题更加复杂，我需要假设每个文档的长度和签名位置都是唯一的。有谁有什么想法吗？
浏览 3提问于2013-10-08得票数 0
1回答
因此，我正在尝试预测某事物(MS_TOT)随时间的变化。我有我的时间变量(ExamStage: 1，3，6和12个月)和我的药物使用变量(acstatus，三个水平的因素)。我认为药物的使用会影响这种变化作为交互作用的时间，我把它放在我的lmer模型中(忽略其他变量)：
>model6<-lmer(MS_TOT~acstatus+ExamStage+acstatus*ExamStage+AIS.1+Level.1+F+ (ExamStage|id), E4)
下面是输出：
> summary(model6)
Linear mixed model fit by R
浏览 3提问于2016-09-22得票数 0
我有一个从数据仓库中提取的非常大的数据集。要将数据集下载到我想要运行lme4的机器上，需要很长时间。我想知道是否可以将数据处理成协方差矩阵，下载该数据(小得多)，并将其用作lme4的数据输入。我已经使用SAS为多个回归模型做了类似的事情，我希望我可以为lme4创建这种类型的输入。
谢谢。
浏览 9提问于2013-02-06得票数 0
我试图在熊猫中实现一个随机效果模型，但我的回归系数与Stata的输出不匹配。我使用的是航空公司路线和机票价格的面板。下面是我的Python代码：
import pandas as pd
import pandas.stats.plm as plm
airline = pd.read_csv("C:...\Airline.csv")
airline['constant'] = 1.0
airline = airline.set_index(['route', 'time'])
airlinePanel = airline.to_p
浏览 2提问于2015-06-27得票数 3
SAS编码：-我对两组(独立但来自同一人群)的差异进行测试。差异量和t-stat匹配的符号(即两组之间的数学差异为负数，tstat为负值)。或者，如果两组之间的数学差异为正，则tstat为正)。
然而，当我进行wilcoxon秩和检验时，我的z-分数的符号与组差的符号(-/+)不匹配。(也就是说，两组之间的数学差异是负的，而z-得分是正的)。如果两组之间的数学差异为正，则z得分为负数)。
我尝试过对数据集进行常规排序和降序排序。
这是我的密码：
*proc sort data = fundawin3t;
by vb_nvb_TTest;
run;
**Wilcoxon rank su
浏览 1提问于2018-03-29得票数 0
回答已采纳
1回答
我对斯塔塔的NBA数据进行了GLS随机效应回归，我被告知这是错误的，因为我没有使用混合效应模型。也许每件事都是这样的，但我对这一解释感到十分困惑。下面是我使用Stata所做的工作：
因变量是球队获胜的变量，自变量是不同类型的对手投篮数据。
使用xtset来解释不同NBA赛季的面板数据
针对固定和随机效应的鲁棒xtreg
Hausman检验有较高的chi^2，所以我坚持GLS随机效应回归。
如果这是不正确的，请让我知道，以便我可以修复这个。当我的因变量是绝对变量时，我一直认为混合效应模型适用于logistic回归。我被告知在您的分析中没有考虑到这个团队分组结构，但是组变量是团队。
浏览 0提问于2017-01-18得票数 3
1回答
我正在寻找一个代码或应用程序，它可以在考虑上下文和运动的情况下从视频中提取出突出的对象，
或
一种仅用于运动显着性地图检测(运动对比度)的算法，因此我可以将其与我所拥有的context_aware显着目标检测器融合。
实际上，我已经测试过context_aware显着性地图检测器，但是它在某些帧中检测到背景的某一部分作为突出对象，我希望在这个检测中涉及到运动和时间，这样我就可以尽可能地提取出确切的突出对象。
有谁可以帮我？
浏览 2提问于2015-08-31得票数 0
回答已采纳
1回答
我要模拟来自多层Logistic分布的数据，每次1000次，每次估计参数，并计算估计值的平均值。但是，在glmer()函数中，结果似乎不能像lm()函数那样提取，比如lm(y~x)$coef。如何从glmer()函数中提取结果？
以下是R代码：
#Simulating data from multilevel logistic distribution
library(mvtnorm)
library(lme4)
set.seed(1234)
## J = number of groups
## n = group size
##
浏览 10提问于2015-07-09得票数 1
回答已采纳
有人知道如何在使用(1 \…)时提取随机项的参数估计值吗？语法在glmer模型(包括se，t比和p值)？我只能访问平均方差和性病。与摘要函数的偏差。
一些背景:我使用了队列和周期随机项(都是因式分解的)，其中周期=每个调查年，队列=8个出生队列。我的模型空模型如下所示: glmer(pid ~ age + age2 + (1欧元队列)+(1个周期))
浏览 3提问于2021-12-06得票数 0
2回答
我试图用plm软件包估计R中固定效应的模型。我的数据如下所示，它是坚定的、城市的、年份的、季度的。其中每一项我都按公司和市级按年季度观察销售和收入。我的回归是收入和销售额。这是销售收入，但希望控制公司，并具体城市的不可观测性。我的实际数据集中有1000+公司。
fid = c(1,1,1,1,
2,2,2,2,
3,3,3,3,3,3,3,3,
4,4,4,4,5,5,5,5,
5,5,5,5)
cityid = c(101,101,101,101,
102,102,102,102,102,102,102,102,103,103,103,103,
浏览 3提问于2015-12-08得票数 1
回答已采纳
我是新的编码以及张贴在论坛上，但我会尽我最大的努力解释问题，并给出足够的背景，以便你能够帮助我完成它。我已经对类似的错误进行了大量的搜索，但是产生它的所有代码和它背后的数据格式都是非常不同的。
我正在处理由各种增长类别组成的生物数据，但我感兴趣的只是长度(代码中的SCL)和年龄(代码中的年龄)。我有许多长度和年龄估计的每一个人通过时间，我是拟合一个线性的nlme模型的少年年龄和Von伯特曲线的成熟年龄。我的少年模型工作得很好，我提取了h(直线的斜率)和t (x截距)。我现在需要使用这些参数以及T(成熟时已知的年龄)来适应成熟阶段。成熟的模型将估计K(这是我唯一的未知)。我已经包含了一个人(ID
浏览 0提问于2016-03-24得票数 1
1回答
我想比较以下模型的估计值，但我不知道是否可以用plm软件包来估计模型(5)：
1/池-> Yit
fit.pooling <- plm(form, data = pdata, model = "pooling")
2在(ind/time/双效) -> Yit - Yi./ Yit - Yt./易-易。- Yt.
fit.fe <- plm(form, data = pdata, model = "within", effect="ind"/"time"/"twoways")
3(ind/时间
浏览 1提问于2018-02-23得票数 0
回答已采纳
3回答
我正在详细介绍我为背景所做的实验--我很清楚lmer的方法，只是不清楚如何提取我需要的一些值/手工计算它们，因此我在SO上发布了这个，而不是简历。我希望这是正确的地方张贴！
。
我的实验有一个平分图设计，有层次:块/图/子图。
有六个街区。每块有两个样地，每个块有两个子样地。处理1有两个层次(A和B)，并在地块一级应用:在每个区块有一个小区接受处理，1个A级和1个接受处理1级B。
处理2适用于子小区级，也有两个级别(C和D)：每个小区有一个子小区接收处理，两个A级和一个子小区接收处理2级B。
实验进行了3年。我感兴趣的是，这两种治疗方法的每一种组合如何影响我的因变量(DV)。
因此，我有4种治
浏览 2提问于2013-08-24得票数 2
回答已采纳
1回答
multi.sanctions.bust.full.ag <- glmer(allbuster ~ lageutradeshare100 + lnlaggdpp + lagtradeopenP + colonial + lagtradesharePT + lnlaggdpt + duration + lndist + nobust + nobustsq + nobustcb + (1 | partnercode) + (1 | case
浏览 5提问于2020-07-31得票数 2
回答已采纳
2回答
我正在用C语言开发一个SLAM算法，并且实现了快速的角点查找方法，给出了图像中的一些强关键点。下一步是获得具有亚像素精度的关键点中心，因此我在每个关键点周围提取一个3x3补丁，并对二维二次型进行最小二乘拟合：
[]()
其中f(x，y)是每个像素的角显着性度量，类似于在上提出的快速评分，但修改后也提供了一个非角像素的显着性度量。
最小二乘：
[]()
[]()
使用
估计的参数。
现在，我可以通过取梯度等于零来计算拟合二次曲线的峰值位置，达到我最初的目标。
这一问题出现在某些拐角情况下，即局部峰值更接近窗口边缘，结果是一个低残差的拟合，但
浏览 9提问于2020-04-14得票数 2
回答已采纳
我试图从混合效应模型(logistic回归)中预测固定效应。我的模型是：
MOD.MIX.1 <- glmer(PATCH_TYPE~PC1+PC2+PC3+JUL.DAY+(1|Study_area)+
(1|ID),family=binomial,data=FOR.MODEL)
我试图用这种方式预测模型：
newdata <- with(MOD.MIX.1, expand.grid(PC1=unique(PC1),
PC2=mean(FOR.MODEL$PC2),
浏览 0提问于2014-11-03得票数 0
回答已采纳
我想从模型公式中提取信息，特别是我想要去除随机效应，以便从混合模型(lme4-表示法)中得到“固定效果部分”。
要做到这一点，我搜索公式中的最后一个+，然后才能找到偏执的(。直到+之前的一切都必须是公式的“固定”部分。对于具有固定效应预测器/变量的模型来说，这很好。
然而，对于空模型(仅在固定效果中拦截)，可能没有+，例如，如果公式是Reaction ~ (Days | Subject)的话。在本例中，我检查是否存在+-sign。但这并不适用于多个随机部件的模型。在下面的示例中，grepl()for f2应该返回FALSE，但返回TRUE，因为在随机部分中，+是为第二个开始的偏执找到的。
我的
浏览 0提问于2019-04-30得票数 2
回答已采纳
当我在plm包中做两种方式的模型时，我的理解是它应该对组和时间有固定的影响，但是当我手动查看固定的效果时，它只对组变量产生固定的效果。例如，在plm包中使用罐装数据：
> data("Produc", package = "plm")
> zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
data = Produc, index = c("state","year"),model='pooling')
浏览 0提问于2019-08-13得票数 0
回答已采纳
3回答
我使用stats() 3.5.2运行manova：
参与者1:20性别在被试因子中的组作为被试因子中的焦虑作为独立测量标准BAC作为受抚养度量<code>H 210</code><code>F 211<//code>
数据集如下：
treat4 = data.frame (
participant = rep(1:20,3),
gender = factor (rep(c(rep("male", 10), rep ("female", 10)),3)),
group = factor (c(rep
浏览 5提问于2020-03-04得票数 0
我试图复制Stata的代码，它估计了一个线性回归模型。
问题是，有两个固定的效应变量(Domaine d‘étude E.F. UniversitéE.F.)。
我现在拥有的是：
import statsmodels.formula.api as smf
results = smf.ols('discriminant ~ diff_eval_formfr + presse + trav_sup + recrut_seul + proced_auditions + taux_insertion_30mois + taux_stable_30 + taux_plei_30 + sal_m
浏览 9提问于2022-06-29得票数 0
SendDelegateMessage:委托在等待10秒后无法返回。主运行循环模式: kCFRunLoopDefaultMode，如果在整个时间间隔内没有使用触摸屏(这会延长等待时间)，请提交一个bug。
- (CGFloat)tableView:(UITableView *)tableView heightForRowAtIndexPath:(NSIndexPath *)indexPath {
/*
//NSString *label = [self.aNote length] == 0 ? kDefaultNoteLabel : self.aNote;
NSSt
浏览 3提问于2009-09-09得票数 0
1回答
我学会了如何使用R来进行F-检验，因为回归模型缺乏拟合，其中$H_0$：“回归模型中并不缺少拟合”。
其中，df_1是SSLF (缺少拟合平方和)的自由度，df_2是SSPE (纯误差平方和)的自由度。
在R中，F检验(例如对于有两个预测因子的模型)可以用以下方法计算。
anova(lm(y~x1+x2), lm(y~factor(x1)*factor(x2)))
示例输出：
Model 1: y ~ x1 + x2
Model 2: y ~ factor(x1) * factor(x2)
Res.Df RSS Df Sum of Sq F Pr(>F)
1
浏览 1提问于2017-07-28得票数 2
回答已采纳
1回答
我试图运行glmmLasso来使用命令估计混合模型：
glm1_final <- glmmLasso(Activity~Novelty + Valence + ROI, rnd =
list(Subject=~1), data = KNov, lambda=lambda[opt],switch.NR=F,final.re=TRUE)
这段代码基本上是从demo("glmmLasso-soccer")中提取的，但是我的变量被替换了。活动是衡量大脑活动的指标，新奇和价是编码刺激类型的分类变量，用于激发反应，而ROI是一个范畴变量，编码大脑的三个区域，我们对这个活动进行了取样
浏览 4提问于2016-06-19得票数 1
回答已采纳
2回答
我有一些平衡的面板数据，并希望将趋势变量包含在我的回归中。然而，我在7年的时间内有60个地区，我不确定如何包括趋势变量。年份变量与预期的一样，对于2005-2011年是重复的。我在思考以下几点：
gen t = .
replace t = 1 if year==2005
replace t = 2 if year==2006
直到2011年，它给了我从1到7的t变量，用于数据中的180个不同的面板。
我的问题是:是否可以像我上面描述的那样包含趋势变量，或者我应该直接将year变量抛入回归？
浏览 0提问于2013-07-27得票数 3
1回答
我需要评估脂肪之间的纵向变化，2次与线性混合效应模型。
我有一些会从第一次访问到第二次访问，因为他们是高血压状态，糖尿病，体重指数，腰围，smoking_status等，以及其他变量不会改变从访问1到访问2，因为他们是性别或种族。
以下变量是虚拟的(高血压状态，糖尿病状况，smoking_status，性别)，而接下来的变量是连续的(bmi，腰围，年龄)。
我最初使用nlme包的方法是：
lme(fat~ diabetes_status + hypertension_status + bmi + waist + smoker + gender + ethnicity + Visit, rand
浏览 2提问于2020-08-12得票数 0
回答已采纳
1回答
我正试着回归一个具有团体时间固定效果的模型，还有很多假人。
egen id_t = concat(id year), format(%15.0f)
areg y u2j* j2j* d1* d2* x1 x2, absorb(id_t) vce(r)
d1和d2是虚拟的，每一个都有数百个可能的值。u2j*是一个变量u2j与时间假人的交互：
forvalues t=1980/2000 {
gen y_`t' = (year==`t')
gen u2jXy`t' = y_`t'*u2j
(...)
我在做这件事时遇到了一个内存错误。我所有的虚
浏览 2提问于2015-05-01得票数 0
2回答
我试图从一个包含许多因素的线性模型中得到拟合值，我想用R包lfe中的lfe函数来估计它。除非我误解了函数返回的fitted.values的含义，否则这些值看起来与手动构建它们时得到的输出不匹配。下面是从包文档中改编的一个示例：
library(lfe)
set.seed(42)
nn = 10
n1 = 3
x <- rnorm(nn)
f1 <- sample(n1, length(x), replace=TRUE)
y <- 2.13*x + cos(f1) + rnorm(length(x), sd=0.5)
est <- felm(y ~ x | f1)
es
浏览 3提问于2017-11-14得票数 3
回答已采纳
1回答
R中具有时不变变量的PLM
、、、、
我正在尝试分析一个面板数据，其中包括美国每个州在45年内收集的观察数据。我有两个随时间变化的预测变量(A，B)和一个不变的变量(C)。我特别感兴趣的是知道C对因变量Y的影响，同时控制A和B，以及不同状态和时间的差异。
这是我的模型，使用R中的plm包。
random <- plm(Y~log1p(A)+B+C, index=c("state","year"),model="random",data=data)
我的推理是，对于时间不变的变量，我应该使用随机而不是固定的效果模型。我的问题是:我的模型和思维是否正确？
提前感谢您的帮助。
浏览 1提问于2014-06-28得票数 0
1回答
可能是最基本的问题，但我不确定该怎么做。如果一个人想要估计一个多级模型，但只有固定的效果，而不是随机的效果估计，你会怎么做呢？例如，以下是我的模型：
model4 <- lmer(m ~ x + ( | id) , data = data_stacked2, REML = TRUE)
谢谢。
浏览 6提问于2020-09-12得票数 1
我想将不同的drl文件加载到kiebase中，如果规则没有附加显著属性，默认执行顺序是否与它们在kmodule中指定的顺序对齐？我的意思是，Drools在幕后，只是从多个drl文件自动创建一个大的drl文件，然后执行。我的理解是，如果规则文件中的规则没有显着性，那么当它们编写时，它们将被执行，除非对现有事实进行任何修改/插入，这可能会重新触发规则的顺序。
谢谢。
浏览 4提问于2021-04-21得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
