我试图以最“准确”的方式理解和使用PISA数据。 最后，我需要的是国家方法(但我想控制个人水平的变量，所以我有国家净效应，因为我想分析教育系统的特征)。我目前正在使用ranef从lmer多级模型中获取它们(使用(1|CNT/SCHOOLID)和(1 |CNT)嵌套)。我现在明白了，我应该使用参议院权重，而不是学生的最终权重(到目前为止，这不是问题，我可以将其包括在lmer中)。如果我使用m
浏览 25提问于2020-07-15得票数 0
我正在试图弄清楚如何从一个简单的lmer模型中组织随机效果的ranef.mer列表，该模型只有一个随机截获和一个变量(性别)。fit.b <- lmer(Math ~ 1 + Sex + (1+Sex|SchoolID), data=pisa_com, REML=FALSE)
我已经使用qqmath绘制了随机效果图，但是我需要能够通过集群编号(在本例中是学校)来标记每个随机效果，或者组织ranef.mer输出。
浏览 1提问于2013-11-07得票数 1
我正在使用几个大型数据库(例如PISA和NAEP)，这些数据库使用具有重复权重和多个似是而非的值的复杂调查设计。我可以使用调查包来解决前者。然而，有没有一个R包/函数来分析后者？
浏览 1提问于2012-06-04得票数 2
回答已采纳
我希望能够输入我将在lmer函数的“weight”选项中使用的变量名。因此，我可以更改数据集，循环“权重”并提取正确的变量。 我想在for循环中提取正确的权重列。因此，对于y，等式为: lmer(y~x+(1|study)，weights = weight.var)和y1: lmer(y1~x+(1|study)，weights = weight.var1)，因此我将权重变量命名为(weight.opt)，然后希望在for循环内的公式
浏览 37提问于2020-10-22得票数 1
回答已采纳
2回答
有人知道如何在Lumley的调查包中使用BRR权重来估计方差，如果您的数据集已经有BRR权重的话？
我正在处理PISA数据，它们已经在数据集中包含了80个BRR副本。
浏览 4提问于2012-10-16得票数 6
回答已采纳
我正在尝试使用lmer()为混合效果模型确定正确的格式。下面的三个模型都试图实现相同的东西-一个简单的固定效果模型-因此肯定有两个是错误的。
在模型中，shoppers是商业中心的消费者支出数据。这反映在gravity模型中，该矩阵将较近的地区区域的权重高于较远的地区区域，以模拟距离的影响。还给出了居民区的Population大小。我认为r3应该是正确的格式，但它比r2产生更高的AIC分数。> r1 <- lmer
浏览 0提问于2012-09-05得票数 1
回答已采纳
在有权重的数据上执行线性混合模型(LMM)时，我遇到了困难(我的意思是不同组的权重不同)。此外，在使用glmulti()时似乎没有保留权重。fm2 <- lmer(angle ~ recipe + temperature + (1|replicate), cake, weights=wght,REML= FALSE)第二个问题发生在使用加权数据进行平
浏览 5提问于2014-10-31得票数 1
1回答
我正在尝试处理一个包含510,000行和636列的数据集。我使用dask dataframe方法将其加载到一个dataframe中，但是不能显示条目。当我试图得到形状时，它会导致延迟。有什么方法可以让我不使用大数据技术来分析整个数据集呢？.zip') == False:
r = requests.get('https://s3.amazonaws.com/udacity-hosted-downloads/ud507&#x
浏览 3提问于2022-07-25得票数 0
我正在使用R Studio中的'intsvy'包来分析PISA。我实际上是在从其他人那里复制代码，所以很想知道为什么下面的命令不能工作(当它在其他编码人员中工作时)：class(DEVCON8a$PVxMATH) for the 5 different plausible values
#[1] "
浏览 1提问于2015-07-10得票数 0
我应该如何在第一个@和空格之间提取文本？我已经在下面提到了我的代码。它提取了所有的单词，但我只想要文本第一的@。@pisa, what's up?理想结果italyregex <- "@([A-Za-z]+[A-Za-z0-9])"
输出@italy @spain
浏览 3提问于2020-04-11得票数 1
回答已采纳
1回答
我正在尝试提取与lmer一起运行的3级分层回归的3级随机效果。问题是，我的数据集被归因于(似是而非的值)，我试图只提取L3随机效果到MItools池中。我可以在一个普通的lmer模型上提取它们，只需要一个PV。m0 <- lmer(PV4MATH ~ 1 + (1 | SCHOOLID) + (1 | CNT), data = pisa, REML = F) ranef(m0)$CNT 但是，当我尝试从一个模型中准确地提取出那些我在所有PV上使用mitool拟合的
浏览 19提问于2020-07-12得票数 0
1回答
我的问题是，我想测试不同组合的固定和随机效果的混合模型，并记录模型的评价度量，以确定哪一个是最佳组合的特征。例如，对于dataset 果树，它具有以下特性：在代码中，与R中的lmer()类似：mod6 = lmer(soldStatus ~ smell+taste + (1|tree), Fru
浏览 1提问于2021-08-27得票数 0
我已经在2013版本中运行了R中的lmerTest和lmer：> data1.frame <- read.delim("colorness.txt", fileEncoding当我检查str(data1.frame)中的数据时，没有任何错误。Error", "df", :
length of 'dimnames' [
浏览 15提问于2014-02-28得票数 0
2回答
我正在使用2010-2014年的5年PUMS数据。我尝试使用复制权重来查找标准误差。标准误差公式在文档中，但我在将其转换为Microsoft Excel或R中的公式时遇到了困难，这两个程序我正在使用。我可以在Excel中对两个变量进行交叉分析，以获得每个复制权重和PWGTP的总和，但我想还有一种更简单的方法。
这里有没有人处理过PUMS数据<em
浏览 3提问于2016-10-18得票数 0
2回答
使用这个虚拟数据集： infected =dummy.data)int.nlme <- interactionMeans(nlme.fit)但这并不是：lmer.model<- as.formula(v
浏览 3提问于2015-07-27得票数 1
回答已采纳
1回答
R中采样纬向的纵向分析
、、、、
另一方面，lme4::lmer占重复测量的比重，而不是抽样权重。任何关于R中的包/函数以及如何使用/编写它们的提示都将不胜感激。我给出了一些具有二项和多项结果的数据集：set.seed(1)
id=rep(1:5,2),
timelmer和glmer来使用变量来计算权重作为预测值。，所以我在
浏览 3提问于2021-07-11得票数 1
1回答
我通常是在广义最小二乘估计框架内工作，这是Wooldridge (2013)的介绍性(2013)所称的，随机效应和固定效应模型，由个人和时间维度索引的纵向数据。我一直在使用plm()中可行的GLS估计，从软件包中，估计随机效应模型--一些统计文献所称的混合模型。plm()函数采用一个index参数，其中我指示单个索引和时间索引。我发现使用来自lmer()包的来建立这样的模型是可能的，但是我对术语和可能性框架的不同感到有点困惑，我想知道是否正确地指定了模
浏览 0提问于2018-03-15得票数 0
回答已采纳
1回答
我在r中使用lmer函数构建了一个模型：这是正确指定的，我得到了我预期的结果。我现在正尝试在spss中复制这些结果。如何在SPSS中复制IV1 + IV1:IV2？
浏览 1提问于2020-06-06得票数 0
我有一个包含测量数据的文件。这里，OnPrescrp是指示受试者是否正在服用处方药的二进制变量，BSW1 - BSW500是引导权重，SURV_WGT是每个受试者的测量权重。大约有20000个条目。例如，有多少17岁到24岁的男性在服用处方药。对于每种类型的计算，我需要一个计数N和95%的CI。我对调查方法一点也不熟悉。相反，在构建最终计数N和置信区间时，我必须考虑测量权重和自举权重。然后在STAT
浏览 5提问于2017-08-12得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
