使用plyr拟合lm模型列表
要使用plyr包拟合lm模型列表，可以按照以下步骤操作：
1. 安装和加载plyr包
首先，确保你已经安装了plyr包。如果没有安装，可以使用以下命令进行安装：
install.packages("plyr")然后加载plyr包：
library(plyr)2. 准备数据
假设你有一个数据框（data frame），其中每一列代表一个因变量，每一行代表一个观测值。例如：
data <- data.frame(
y1 = c(1, 2, 3, 4, 5),
y2 = c(5, 4, 3, 2, 1),
x = c(1, 2, 3, 4, 5)
)3. 使用plyr拟合lm模型列表
使用dlply函数可以方便地对数据框的每一列进行操作，并返回一个列表，其中每个元素是一个lm模型。
models <- dlply(data, .(variable), function(df) {
lm(y ~ x, data = df)
})在这个例子中：
dlply函数会对data数据框的每一列进行操作。.(variable)指定了分组的变量名。这里我们假设每一列的名称是y1,y2, 等等。- 函数内部的
lm(y ~ x, data = df)会为每一列拟合一个线性回归模型。
4. 查看模型结果
你可以遍历模型列表并查看每个模型的摘要：
lapply(models, summary)这将输出每个模型的详细摘要信息。
完整示例代码
# 安装和加载plyr包
install.packages("plyr")
library(plyr)
# 准备数据
data <- data.frame(
y1 = c(1, 2, 3, 4, 5),
y2 = c(5, 4, 3, 2, 1),
x = c(1, 2, 3, 4, 5)
)
# 使用plyr拟合lm模型列表
models <- dlply(data, .(variable), function(df) {
lm(y ~ x, data = df)
})
# 查看模型结果
lapply(models, summary)通过这种方式，你可以轻松地对数据框中的每一列进行线性回归分析，并得到一个包含所有模型的列表。
相关·内容
2回答
我正试图结合使用lm()函数和plyr包来拟合数据中的每个片段的模型，因为我的数据是按键分组的。#Sample data
test_data <- data.frame(key = c("a","a","a","a","a&qu
浏览 5提问于2016-03-17得票数 0
回答已采纳
1回答
一旦拆分，我基本上想在子集上运行线性回归模型。下面是我正在使用的代码，但它绝对不正确。此外，如果我能将模型结果输出到新的数据框中，这样每个站点都有一行模型参数估计，那就太好了--这只是一个愿望，目前并不是必须的。谢谢你的帮助！les_events <- split(les, les$link5)
coef <- coef(
浏览 4提问于2013-04-16得票数 0
我用lm和plyr函数按因素级别拟合了一些回归，但在尝试访问summary.lm类的内容时我发现了一些问题。sum = dlply (mtcars, .(cyl), lm, formula = mpg ~ wt) %>%这给出了按cyl变量级别进行的每个回归的总结。在这个例子中，如果我想在fstatistic对象中获得summary.lm元素，我唯一的机会是手动编写类似sum$'4'$fstat
浏览 1提问于2017-09-15得票数 1
回答已采纳
2回答
如何修改dlplyr帮助页中的以下代码，以便连续地将没有截取的线性模型拟合到data.frame中的多个预测器。#From the dlplyr help pagelm(rbi ~ year, data = mutate(df, year = year - min(id), linmod)library(plyr)
grouping<-factor(sample(c('A', 'B', 'C&#x
浏览 6提问于2015-02-24得票数 0
回答已采纳
1回答
我在一个变量(a)上拟合了几个模型作为mod1。Output mod1包含5个模型的回归拟合列表。我在模型列表中做了几个操作。func <-function(z){ fit2 <- lm( y~x + I(z^2)) fit4 <- lm( y~ns(x
浏览 24提问于2019-06-02得票数 0
回答已采纳
<- c(rep("a", 4), rep("b", 4), rep("c", 4))我想为x和y共享相同ID的子集创建一个单独的lm下面的代码完成了这项工作：b.lm <- lm(x~y, data=subset(df, ID=="b&q
浏览 0提问于2011-09-14得票数 6
回答已采纳
1回答
R公式
、、
对于lm()来说，事情很好，但是对于survfit()来说，它们失败了。我正在根据一些数据拟合一系列公式。因此，我用作为变量传递的公式调用建模函数。稍后，我想处理拟合对象的公式。详细实例library("plyr")
lm(mpg ~ wt, data =formula
##
浏览 0提问于2016-10-16得票数 1
回答已采纳
2回答
我喜欢在我的数据集上运行method=“like”，如果出错，我喜欢尝试lm方法。下面是我的代码来做这件事，但它仍然失败，有什么想法吗？theme_solarized(light = FALSE) + geom_smooth(method=ifelse(class(try(loess(Duration~Date, df)))=="try-error", "lm
浏览 3提问于2012-12-15得票数 2
回答已采纳
1回答
3L, 4L, 5L, 6L, 69L, 70L, 71L, 72L, 73L, 74L), class = "data.frame")我可以接近，但我不认为这是正确的，因为它使用的是Tissue type (just one here for e.g.)#Perform Anova
anova(lm
浏览 1提问于2014-05-30得票数 0
回答已采纳
在R中，我试图动态地改变我的线性模型中的变量。我已经保存了一个变量的字符向量，我想在我的lm中使用它作为调节变量。： standardModel <- lm(data=studyData, yVar ~ xVar)
summary.aov(standardModel) 我想从zList中动态地包含一个与此模型一起使用的调节变量列表因此： zList <- c("numVar1","numVar2","facV
浏览 14提问于2020-06-20得票数 0
回答已采纳
我正在尝试使用tapply()函数按几个类别运行模型，但效果不是很好。我的数据有20个诊所，我想按每个诊所运行模型。这是我的模型：rq(logA~ dose+ chtcm + cage +raceth + sex,tau=.9)
我的数据作为一个变量(值为1-20)。有没有人知道如何像在其他统计软件包中一样，在R中运行这个模型？
浏览 0提问于2013-01-17得票数 0
3回答
我的任务是使用70年代、80年代、90年代和00年代这4个不同时期的值，将每个国家的股票收益(因变量)与美国股票收益(自变量)进行回归。数据集(.csv)可从下载谁能提供我必须采取的步骤，以编写一个函数，将达到预期的结果？一些R代码的片断也会受到欢迎。我也在考虑
浏览 0提问于2013-01-29得票数 1
回答已采纳
1回答
我想用lm来拟合R中的线性模型，得到系数估计和总模型拟合的p-值+p-值，所以基本上是summary.lm的输出。有没有办法‘欺骗’lm，给
浏览 0提问于2016-02-18得票数 6
回答已采纳
1回答
用R？
、、、
如何使用R计算每个系数在N=1000最大似然比模型拟合上的平均值和标准差？, sd=e.sd) data.frame(X1,X2,Z) coef(lm(Z~X1+X2,data=samples))raply(1000,statfun(simfun()))
浏览 4提问于2013-10-01得票数 0
1回答
我刚刚开始学习R，并试图展示汽车的每加仑数和马力之间的关系。我试图得到一条最适合的线来匹配数据，这样我就可以用它来预测其他值。 #read data from csv file df <- head(data, 350)
x <- df$horsepower x <- strtoi(x)
#create scatterplot
浏览 13提问于2021-04-07得票数 0
我在尝试弄清楚如何使用plyr而不是循环传递给predict.lm的新数据来获得预测的逻辑。有人能帮上忙吗？示例：require(plyr)
# settin
浏览 0提问于2013-09-05得票数 1
许多函数的输出都是以列表结构形式出现的--例如lm()。因此，您可以使用"$“或使用方括号进行索引来检索输出的单独部分。我的问题是如何以列表形式创建输出，而不公开显示它是一个列表。如您所知，当列表打印在屏幕上时，它通常具有指示的子列表名称或索引，如下所示：L1
[1] 1 2 3
[1] "a""b" "c&q
浏览 0提问于2012-11-02得票数 0
回答已采纳
我有一个回归模型的列表，所有的回归模型都有相同的词数(也就是相同数量的预测变量)。从实质上讲，它们都有不同的模型术语，这是正确的。一些指示性数据d1 <- data.frame(y = rnorm(100), x2= runif(100), x4 = runif(100))
拟合<e
浏览 1提问于2013-08-23得票数 2
回答已采纳
1回答
使用相同的因变量，使用相同的模型矩阵，得到相同的精确拟合。让我们使用维基百科文章中的定义来确定决定系数(如果您愿意的话解释的方差百分比)：[1(lm2_)^2) / (sum(fitted(lm2_)^2) + sum(residuals(lm2_)^2))因此，很明显
浏览 0提问于2013-06-03得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
