如何使用huxtable()在tibble中打印回归模型的列表-列
huxtable()是一个R语言中的包，用于在tibble中打印回归模型的列表-列。它提供了一种简单的方法来创建和格式化表格，以展示回归模型的结果。
使用huxtable()在tibble中打印回归模型的列表-列，可以按照以下步骤进行：
- 首先，确保已经安装了huxtable包。可以使用以下命令安装：
install.packages("huxtable")- 导入huxtable包：
library(huxtable)- 创建一个包含回归模型结果的tibble。假设我们有一个包含回归模型结果的数据框df，其中包含了模型的系数、标准误差、t值和p值等信息。
df <- tibble(
Coefficients = c(1.23, 2.34, 3.45),
Std.Error = c(0.12, 0.23, 0.34),
t.value = c(4.56, 5.67, 6.78),
p.value = c(0.001, 0.01, 0.05)
)- 使用huxtable()函数创建一个huxtable对象，并将数据框df作为参数传递给它：
ht <- huxtable(df)- 可以使用huxtable包提供的一些函数来格式化表格，例如设置表头、设置列名等。以下是一些常用的函数：
- set_all_borders()：设置所有边框
- set_bold()：设置粗体文本
- set_font_size()：设置字体大小
- set_width()：设置列宽度
例如，可以使用以下代码设置表头和列名：
ht <- set_header(ht, "Regression Model Results")
ht <- set_colnames(ht, c("Coefficients", "Std.Error", "t.value", "p.value"))- 最后，使用print()函数打印huxtable对象，以在控制台中显示表格：
print(ht)这样就可以在tibble中打印回归模型的列表-列了。
对于huxtable()在tibble中打印回归模型的列表-列的应用场景，它可以用于展示回归模型的结果，包括系数、标准误差、t值和p值等信息。这对于数据分析、统计建模和机器学习等领域非常有用。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云计算服务：https://cloud.tencent.com/product/cvm
- 腾讯云数据库：https://cloud.tencent.com/product/cdb
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网：https://cloud.tencent.com/product/iot
- 腾讯云移动开发：https://cloud.tencent.com/product/mobdev
- 腾讯云存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链：https://cloud.tencent.com/product/baas
- 腾讯云元宇宙：https://cloud.tencent.com/product/mu
相关·内容
我正在尝试生成一个Rmarkdown文档，该文档打印我目前存储在tibble的list-column中的一系列回归模型。回归模型如下所示。 #generate data
covar1<-rnorm(100)
covar2<-rnorm(100)
depvar1<-rnorm(100)
depvar2<-rnorm(100)
#generate models
model1<-lm(depvar1~covar1)
model2<-lm(depvar1~covar1+covar2)
model3<-lm(depvar2~covar1)
model4<
浏览 6提问于2019-09-12得票数 0
回答已采纳
2回答
我通常希望对数据的两个不同子组和多个模型(例如，双变量模型和带有控件的模型)运行相同的回归。在我看来，嵌套、映射和整理的purrr/tidyr/broom工作流应该能很好地工作。我了解如何在嵌套工作流中创建子组，但我不知道如何运行多个模型并输出每个模型的有序回归结果列表。
例如：
mtcars %>%
nest(data=-c(vs)) %>%
mutate(
fit = map(data,~lm(mpg ~ cyl, data = .x)),
fit1 = map(data,~lm(mpg ~ cyl + gear + wt, data = .x)),
浏览 4提问于2020-08-03得票数 0
回答已采纳
1回答
R中数据帧的回归
、、
嘿，我有以下测试数据：
test = data.frame(Date = c(as.Date("2010-10-10"), as.Date("2010-10-10"), as.Date("2010-12-10"), as.Date("2010-12-10")), Rate = c(0.1, 0.15, 0.16, 0.2), FCF = c(100,200,150,250))
现在我想按日期将数据分组，在每个组中进行线性回归FCF~Rate，然后在每个组中执行这些回归，以获得每个日期和比率的回归值。我有以下代码：
output =
浏览 3提问于2021-05-05得票数 2
回答已采纳
我在R Markdown中工作，建立了几个回归模型，并用huxreg制作了一个汇总表。这些模型有很多特性，因此汇总表非常大。当我检查代码中的输出时，它看起来完全没有问题。但是，当我将代码“编织”成pdf时，表格很大，它既不适合pdf页面的宽度，也不适合pdf页面的高度。它也不会延长到下一页。 代码： {r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, error=FALSE, warning=FALSE, message=FALSE)
setwd('/Users/ladarudnitckaia/Data')
浏览 35提问于2020-04-09得票数 0
我有一个包含3个回归模型的列表，称为logregs。我的数据有一个名为type的列，它只有整数1、2和3，用于确定应该使用logregs中的哪个回归模型，还有一个名为adstock的列，它是回归模型中使用的唯一自变量。
我正在尝试做一些类似这样的事情：
dataframe %>% mutate(probability = predict(logregs[[type]], type = "prediction", newdata = adstock) )
示例数据帧：
structure(list(type = c(3L, 3L, 3L, 3L, 3L, 3L), adst
浏览 9提问于2020-05-22得票数 0
回答已采纳
1回答
潮间带双重筑巢
、、
使用Wickhams介绍的在R中提供数据科学，我正在尝试创建一个双重嵌套列表。
library(gapminder)
library(purrr)
library(tidyr)
gapminder
nest_data <- gapminder %>% group_by(continent) %>% nest(.key = by_continent)
我如何进一步嵌套这些国家，使nest_data包含by_continent和一个新级别的嵌套by_contry (最终包括Tibbleby_year)？
此外，在为gapminder数据创建此数据结构之后--您将如何从中为每个国
浏览 3提问于2016-08-30得票数 7
回答已采纳
1回答
我正在寻找一种在R中用扫帚和dplyr进行“多步”回归的方法，我使用“多步”作为回归分析的占位符，在回归分析中，您可以将以前回归模型的最终回归模型元素(如拟合或残差)进行集成。这种“多步”回归的一个例子是工具变量(IV)回归的2 2SLS方法。
我的(分组)数据如下所示：
df <- data.frame(
id = sort(rep(seq(1, 20, 1), 5)),
group = rep(seq(1, 4, 1), 25),
y = runif(100),
x = runif(100),
z1 = runif(100),
z2 = runif(100)
浏览 4提问于2021-12-09得票数 2
回答已采纳
我有一个由93个变量的308个观察值组成的R数据帧(即，组织成93列的308行数据)。一列是位置，有19个因素，另一列是观察日期。其余91列(3到93)包含各种化学品的浓度测量。 我想针对采样日期(x变量)对浓度测量值(y变量)的每一列应用最小二乘回归(lm函数)，并按位置对回归进行分组。 我尝试了以下代码，它适用于一列浓度测量，给出了19个列表中每个位置的lm函数结果-但是我如何调整代码，使其在所有91列测量上运行？ df.2 <- split(df, f=df$location)
model_lm <- lapply(df.2, function(x)
浏览 21提问于2021-09-12得票数 0
我有一个学生属性和考试成绩的数据，我试图拟合每个年级的线性模型(从1到12)。我正在使用扫帚包，以有效地创建一个模型，为每个年级。下面是一个简化的示例数据集和我正在使用的代码。
#start df creation
grade <- rep(1:12, each = 40)
attendance_rate <- round(runif(480, min=25, max=100), 1)
test_growth <- round(runif(480, min = -12, max = 38))
binary_flag <- round(runif(480, min =
浏览 2提问于2020-06-09得票数 1
回答已采纳
“图像中的数据结构”
RStudio中的问题
我有33个数据帧的列表：
长度为33的growth_in list33列表
18005 list122 x 16 tibble 176行和16列
18013 list122 x 16 tibble 176行和16列
18019 list122 x 16 tibble 176行和16列
18025 list122 x 16 tibble 176行和16列
.依此类推，直到33个数据帧(总共)
我可以导出多个CSV文件，但我需要一个具有不同表名的excel文件(即一个包含33张工作表的excel文件)。在使用excel_export命令后面对下面提到
浏览 2提问于2019-10-22得票数 0
1回答
我在一个变量(a)上拟合了几个模型作为mod1。Output mod1包含5个模型的回归拟合列表。我在模型列表中做了几个操作。现在我想像fit1、fit2、...fit5等一样，在一个回归模型输出中取消列出mod1。 func <-function(z){
fit1 <- lm( y~ x + z )
fit2 <- lm( y~x + I(z^2))
fit3 <- lm( y~poly(x,3) + z)
fit4 <- lm( y~ns(x, 3) + z)
fit5 <- lm( y~ns(x, 9) + z)
return
浏览 24提问于2019-06-02得票数 0
回答已采纳
我正在用一组一致的自变量对结果变量列表进行多变量回归。对于单变量回归，我遵循在嵌套数据框架上使用来自gtsummary的tl_uvregression，但我试图将其推广到嵌套数据框架上使用tbl_regression的多变量回归，当我试图对表进行unnest时，我得到了“输入必须是向量列表”的错误。下面是我尝试过的--我想我错过了一些很小但很关键的一步，但我不知道它是什么。我想要的输出是一个多变量回归输出表，其中每个模型作为列，所有协变量作为行(类似于在这些多变量模型的列表上执行tbl_merge，这些模型在tbl_regression中单独运行)。
library(tidyverse)
li
浏览 5提问于2021-02-20得票数 1
回答已采纳
3回答
R循环残差回归
、、、
data(mtcars)
head(mtcars)
Y = c("mpg")
X = c("cyl", "disp", "hp")
for(V in 1:length(X)){
MODEL=lm(Y~X[V],data=mtcars)
mtcars$paste(X[V], "residual")=MODEL$resid
}
我有数据mtcar，并希望估计许多回归模型，并将残差作为新变量存储在数据mtcar中。我在下面提供了一个例子。基本上，我希望分别回归所有三个变量预测器上的'mpg‘变量，因此生
浏览 4提问于2021-10-09得票数 1
正如所解释的，我们可以使用来自broom包的tidy函数来跨组运行回归。例如，下面列出了iris dataset的演示代码，但是如果我们打算同时遍历多列并使用不同的因变量(Sepal.Length、Sepal.Width、Petal.Length)和这个group_by操作运行回归，我如何将(S)应用函数集成到这种情况中并获得这些回归模型(3*3=9)的结果呢？
library(dplyr);library(broom)
res1=iris%>%
group_by(Species)%>%
do(res=lm(Sepal.Length~Petal.Width,data=.))
tid
浏览 0提问于2019-07-26得票数 1
回答已采纳
2回答
如何替换矢量整理器
、、、、
我正在寻找另一个函数，它可以在broom::tidy()函数被删除后替代它。下面是扫帚包装警告的内容：
整洁的原子矢量清理器被废弃，并将从即将发布的扫帚中移除。
以下是对功能的描述：
tidy()生成一个tibble()，其中每一行都包含有关模型的一个重要组件的信息。对于回归模型，这通常对应于回归系数。如果您想检查模型或创建自定义可视化，这可能很有用。
谢谢你，
约翰
浏览 4提问于2020-08-26得票数 2
我试图按组运行回归，然后将回归模型对象传递给另一个函数。
library("lmtest")
library("broom")
library("tidyr")
library("dplyr")
library("purrr")
fitted_models <- mtcars %>%
group_by(gear) %>%
do(fit = lm(mpg ~ ., data = .),
test = coeftest(fit)) # from lmtest
浏览 0提问于2018-03-01得票数 3
回答已采纳
3回答
按组添加模型预测
、、、、
我在我的数据集中按组估计回归模型，然后我希望为所有组添加正确的拟合值。
我正在尝试以下几点：
library(dplyr)
library(modelr)
df <- tribble(
~year, ~country, ~value,
2001, "France", 55,
2002, "France", 53,
2003, "France", 31,
2004, "France", 10,
2005, "France", 30,
2006, "France&#
浏览 1提问于2018-04-19得票数 2
回答已采纳
2回答
我有一个数据，让我们称它为因变量，各种自变量(指标)和过滤变量。我的目标是通过过滤变量中的不同类别来运行回归。例如，如果我想为code == "all"运行回归，我只需取我的数据，过滤代码，并运行一个回归：
sample_tib %>%
filter(code == "all") %>%
glm(love ~ ., data = ., family = "gaussian")
但是，我面临着几个问题：
在我的示例中，--我的glm()将接受所有列，不包括code。回归的理想输入是love ~ ind1 + ind2 +
浏览 3提问于2021-06-17得票数 2
回答已采纳
2回答
我正在运行一个国家固定效果模型，进一步包括一个虚拟变量。我正在使用Stargazer软件包，但似乎不知道如何同时报告置信区间和确切的p值。如果我运行我的模型：
stargazer(dummy_CPP, title = "xx",align = TRUE,
no.space = TRUE, ci = TRUE,
report = ('vc*p'), single.row = TRUE)
只有p值在我的输出中。此外，由于我有30+变量，所以即使我使用single.row = TRUE，整个表格也不能放在一页中，因为p值是在
浏览 1提问于2018-04-03得票数 2
1回答
我正在计算45个不同id值的因变量的线性回归的模型输出。如何使用整洁(dplyr，等)完成这个任务的代码？
我有一个包含三个变量data = c(id、distance、actPct)的数据集，使得id == 1:45；-10 <= distance <= 10；0 <= actsPct <= 1。
我需要在model0n的每个值上运行一个回归( id )，这样model0n就有了一个新的tibble/df。我已经完成了一次回归：
model01 <- data %>%
filter(id == 1) %>%
filter(distance <
浏览 3提问于2019-01-14得票数 0
回答已采纳
1回答
我列出了14个因变量，其中我正在运行相同的回归模型(相同的模型类型，相同的自变量)。我使用gather作为单个outcome变量获取所有因变量，在该变量上运行tbl_uvregression或tbl_regression，然后使用gtsummary包中的tbl_stack来组织输出。我试图弄清楚如何使用每个模型的outcome变量的值来命名每个表。我知道我可以将一个名称列表传递给tbl_stack(group_header)，但是我注意到这很容易出错，因为我必须小心outcome变量中的值是如何排列的，然后确保我按相同的顺序输入它们，而且我已经犯了很多错误，所以我担心这种方法。有没有一种方法可
浏览 2提问于2021-03-07得票数 2
回答已采纳
感谢来自新手的模型摘要包的创建者。
我调用datasummary_df是为了在RMarkdown表中打印描述性统计信息，并且希望在某些列中不显示显着数字，在其他列中不需要两个重要数字。Datasummary_df只有一个fmt参数，但是在调用模型上的模型摘要函数时，这个用于这个参数。我不知道如何在调用只打印dataframe (datasummary_df)而不是使用模型的函数时应用这个答案。
我也希望允许列标题包装，但我不确定如何做到这一点。
谁能给我指明正确的方向？非常感谢。
浏览 1提问于2022-03-10得票数 0
我浏览了官方的，除了每个回归表的一个主标题(caption或set_caption)之外，我没有看到任何关于如何设置其他标题的内容。如果任何人知道如何在表格顶部的主标题下设置副标题，我将不胜感激。
这是一个回归表的玩具示例。我可以很容易地把标题放在最上面。如何添加字幕？
data(mtcars)
model1 <- lm(mpg~cyl, data = mtcars)
model2 <- lm(mpg~cyl + hp, data = mtcars)
huxtable::huxreg(model1, model2) %>% set_caption("TABLE A
浏览 1提问于2019-11-08得票数 1
5回答
有没有一种方法可以以更易读的方式将小data.frames打印到控制台？
例如，是否可以输出到控制台：
library(MASS)
iris[1:5, ]
Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1 5.1 3.5 1.4 0.2 setosa
2 4.9 3.0 1.4 0.2 setosa
3 4.7 3.2 1.
浏览 3提问于2015-09-03得票数 27
回答已采纳
3回答
整齐的多元回归
、、、
db = tibble(a = rnorm(100), b = rnorm(100), c = rnorm(100))
如果我想要一个整齐的多元线性回归，我可以去
lm(data = db, 0 + a ~ b + c) %>% tidy()
但是如果我想要多个单变量回归，我会去
lm(data = db, a ~ 0 + b) %>% tidy() %>%
add_row(lm(data = db, a ~ 0 + c) %>% tidy())
现在，考虑到许多回归者列，我想避免将每一个回归者编码为一个新的add_row，我应该如何使代码更加合成呢？
这里有一个部
浏览 13提问于2022-03-27得票数 1
回答已采纳
有人能就何时使用map() (所有map_.()函数)以及何时使用summarise_at()/mutate_at()提出建议吗？
例如，如果我们对向量列做了一些修改，那么我们不需要考虑map()？如果我们有一个df /有一个列，其中包含一个列表，那么我们需要使用map()
map()函数总是需要与nest()函数一起使用吗？任何人都可以建议一些有关这方面的学习视频。另外，如何将列表放在df中，同时建模多个列表，然后将模型结果存储在另一列中？
非常感谢!
浏览 2提问于2018-05-04得票数 7
回答已采纳
1回答
我有一个数据集，我正在用线性回归(只是通过geom_smooth())来绘制。我还列出了其他线路的坡度和截距。我想知道测试这些其他行是否落入回归的CI中的最好方法是什么。真正的数据集将有许多行要测试，所以我希望有一些方法可以集体测试。 library(tidyverse)
# sample data
set.seed(0)
df = tibble(
x = runif(20,min=5,max=10),
y = x + rnorm(20,mean=2,sd=3)
)
# df of lines
lines.df <- data.frame(line = c("
浏览 10提问于2020-04-29得票数 2
回答已采纳
2回答
我希望我的函数打印漂亮的表，无论它是从基R调用还是从RStudio笔记本(.Rmd文件)调用。该函数应该知道从哪里调用它，并相应地调整表。我希望该函数易于使用，并且不希望用户指定调用函数的位置。
我可以使用huxtable实现其中的一些功能，但是用户仍然需要稍微修改代码。(我认为这类似于kable。)
以下是函数定义：
library(huxtable)
func = function() {
table = hux(head(iris))
# Color of table border: white on screen / base R, black in HTML
浏览 11提问于2022-02-09得票数 0
1回答
我想调整jtools回归表中的行填充。我尝试了set_row_height，因为我读到底层结构是huxtable。但我不确定，但没起作用。
---
title: "Untitled"
author: "Name"
date: "10 5 2021"
output: html_document
---
## R Markdown
```{r, warning=FALSE, message=FALSE}
mymodel <- lm(mpg ~ .，data=mtcars)
图书馆(Tidyverse)
图书馆(Jtools)
图书馆(Hux
浏览 6提问于2021-05-10得票数 1
回答已采纳
2回答
是否有一种方法来计算个人的回归系数，而不是只计算一个群体的回归系数。计算一个很小样本的回归系数？
背景
我的目标是测试不同的方法来聚集/分割一群人。在所讨论的领域中，用于分段的最佳变量将是最能表达个人在某些属性上的重要性的变量。
然而，我没有对这一重要性的直接衡量。我要做的是，根据属性等级的shapley值回归和最终的应答选择，对整个样本进行间接测量。
很明显，我问了几个人，他们会如何评价一个品牌的属性，以及他们最终会选择哪个品牌。
现在，我想知道是否也有可能在个人层面上获得这种间接重要性衡量，这意味着：
从整体回归模型中提取单个系数估计量(我认为这是不可能的)
计算每个人的回归(不知道如何做
浏览 0提问于2020-06-08得票数 1
2回答
我正在用R/Sweave写一本会议摘要小册子。我已经制作了印刷程序手册，其中只包含id，作者，标题。
现在，我想修改它以包括摘要(不是打印)。但摘要是冗长的。我的想法是使用带有抽象信息的单元格，并将其显示在包含作者信息的行下面--扩展到页面的整个宽度。
因此，每一个其他行都只有一列跨越整个表的宽度。是否有方法将multicolmn{x}添加到其他行？
如果无法找到解决方案，那么建议如何以一种很好的方式打印完整的摘要是值得欢迎的。(“只需使用景观”或“调整列宽”)
而且，它不必是PDF格式的。我可以切换到markdown/html，让它看起来更接近真正的会议日程安排，上面有完整的摘要。同样，
浏览 2提问于2017-05-01得票数 0
1回答
我的数据是由ID指定的大型参与者集合的x值(shift.magnitude)和y值(response.sign)的集合。我使用以下命令分别对每个参与者执行逻辑回归：
do(fitmod = tidy(glm(response.sign~shift.magnitude,family=binomial(link="logit"), data = .))) %>%
unnest(fitmod) %>%
pivot_wider(id_cols = ID, names_from = term, values_from = c(estimate))
我想使用glm
浏览 2提问于2022-05-10得票数 0
通过在不同的表中给出的颜色来编码表(用huxtable生成)是可能的吗？
head_mtcars <- head(mtcars)
tbl_A <- as_hux(head_mtcars[1:2])
tbl_B <- as_hux(head_mtcars[3:4])
( 1)首先，我想为每一列分别用颜色空间对tbl_A第1和第2列中的元素进行着色
tbl_A %>%
huxtable::add_colnames() %>%
map_background_color("for each column of tbl_A", by_colors
浏览 3提问于2019-11-25得票数 0
回答已采纳
1回答
所以我有这样的数据：
fruit cost count
banana 23 1
orange 13 13
grape 10 32
banana 64 42
orange 23 24
grape 10 2
banana 112 12
orange 23 42
grape 64 1
基本上，我想为每个水果建立一个简单的线性回归模型。
因此，香蕉的回归模型将只使用：
banana_df
fruit cost count
banana 23
浏览 3提问于2020-11-19得票数 0
回答已采纳
我见过一个列表应用(lapply)的示例，它很好地获取数据对象的列表，并返回回归输出的列表，我们可以将该列表传递给Stargazer以获得格式良好的输出。Using stargazer with a list of lm objects created by lapply-ing over a split data.frame library(MASS)
library(stargazer)
data(Boston)
by.river <- split(Boston, Boston$chas)
class(by.river)
fit <- lapply(by.river, f
浏览 22提问于2017-02-08得票数 4
回答已采纳
如何将huxtable保存为LaTeX (.tex)格式的文件？我找到了将LaTeX表打印到控制台的this code： library(huxtable)
ht <- huxtable(
a = 1:3,
b = letters[1:3]
)
print_latex(ht) huxtable包有没有将表保存到磁盘的原生方法，或者我们需要一个变通方法？
浏览 16提问于2021-07-07得票数 0
回答已采纳
1回答
不使用整齐提取限制回归系数
、、、、
我使用restriktor包来执行受限的回归；但是，与此同时，我使用dplyr按组执行受限的回归。为了提取这些系数并将它们格式化成一个很好的面板格式，我使用了整洁和扫帚，但是这个整洁包装不能在restriktor上工作，所以我不知道如何提取这些系数：
library(restriktor)
library(dplyr)
reg =
mtcars %>%
group_by(cyl) %>%
do(model = restriktor(lm(mpg ~ wt + hp, data =.), constraints = ' wt < -4 '))
我
浏览 4提问于2022-02-19得票数 1
回答已采纳
我希望将(mutate)多列附加到数据框架中，其中这些列存储在矩阵中。是否有一种使用tidyverse中的函数来完成此操作的方法？(但请注意，通过诉诸base::函数是可能的。)同样地，我要问的是使用tidyverse中的函数来实现这一点的最自然(或惯用的)方法是什么。
例如，假设我们估计了一个分位数回归：
library(dplyr)
tibble(x = runif(100)) %>%
mutate(y = rnorm(n())) ->
EstimationData
library(quantreg)
taus <- (1:9)/10
rq_fit <
浏览 0提问于2018-08-20得票数 2
回答已采纳
2回答
我知道有几种方法可以比较回归模型。创建模型(从线性到多个)和比较R2、调整后的R2等的一种方法：
Mod1: y=b0+b1
Mod2: y=b0+b1+b2
Mod3: y=b0+b1+b2+b3 (etc)
我知道有些包可以进行逐步回归，但是我试图用purrr来分析这个。我可以创建几个简单线性模型 ()，现在我想知道如何创建回归模型将特定的IV添加到方程中
可复制码
data(mtcars)
library(tidyverse)
library(purrr)
library(broom)
iv_vars <- c("cyl", "disp", "
浏览 3提问于2017-11-28得票数 5
回答已采纳
我正在尝试进行回归分析，然后将结果保存在csv文件中。
两个问题：
为什么在单个单元格中以csv格式打印回归分析结果集/列表(regr.predict来自sklearn import linear_model)？我试图将每个预测都打印在另一行下面的一个单独的行中(在我的例子中，我有418个结果，所以应该是418行在同一列中)。以prediction1=regr.predict(xtest)作为f:字段名=‘PassengerId’，‘幸存者’= csv.DictWriter(f，字段名) writer.writeheader()作为prediction1: writer.writerow({‘
浏览 4提问于2021-01-06得票数 0
编辑:我转而使用mtcar数据集，因为它更好地代表了我真实数据的结构。忽略下面的虹膜例子：
我希望运行几个具有多个因变量和自变量的回归模型。
说我在处理虹膜数据集。
我的因变量将是Petal.Length和Petal.Width。我的自变量是Sepal.Length和Sepal.Width。
我想为数据集中的每个物种分别运行以下模型(setosa、versicolor、virginica)：
lm(Petal.Length ~ Sepal.Length)
lm(Petal.Length ~ Sepal.Width)
lm(Petal.Width ~ Sepal.Length)
lm(Petal
浏览 5提问于2022-03-14得票数 0
回答已采纳
我正在尝试运行一个简单的曲线图列表，来评估R中赋值的正确线性回归模型。基本上，在11列中，第一列是y，不同的列x= 3:11是θ，或者y所依赖的变量。第2列的分类值基本上为1-9 (9个类别)。
我正在尝试创建图形矩阵。我可以一个接一个地做图形，这应该适用于赋值，但我试图使用一个循环来完成它。代码如下：
第一列是y，从属和第3-11列是变量。如果我省略了根据第二列的着色，我可以用一个基图来完成它。如果我尝试上面的代码，我得不到任何结果
par(mfrow(3,3))
for(i in 3:11){
ggplot(hd,aes(x=hd[1,],y=hd[,i],color=hd[,2]))
浏览 16提问于2019-11-12得票数 0
我正在尝试编写一个函数，它将运行多个回归，然后将输出存储在一个向量中。我想要的是函数从我将提供的列表中选择因变量，然后在相同的右侧变量上运行回归。不知道该怎么做。如有任何提示将不胜感激。
my_data <- data.frame(x1=(1:10) + rnorm(10, 3, 1.5), x2=25/3 + rnorm(10, 0, 1),
dep.var1=seq(5, 28, 2.5), dep.var2=seq(100, -20, -12.5),
dep.var3=seq(1, 25,
浏览 5提问于2021-09-18得票数 1
回答已采纳
1回答
我有一个要打印的huxtable，它在第一页的末尾。由于其大小和宽度，huxtable将打印在第二页上。有没有什么可行的方法可以拆分表格的一部分，并将表格的前半部分投影到第一页，将下半部分投影到第二页？我尝试了longtable和其他tex方法。此外，我不希望使用kable包，因为它降低了表的美观性。一个样例Huxtable演示作为一个解决方案将对我有很大的帮助，并将非常感谢。谢谢。
浏览 11提问于2020-10-30得票数 1
1回答
我试图在R中的一个数据帧上运行多重回归，但在如何迭代每列中的每个字符串时，我有很多令人困扰的想法。 我有以下数据帧： Category Zone Season P1 P2 P3 Value
-------------------------------------------------------------
Blue D1 Winter 1 4 4 55
Blue D1 Winter 3
浏览 10提问于2021-09-02得票数 1
我有一个数据集，类似于下面生成的tb：
x = seq(0, 10, by=.4)
y1 = x + rnorm(x)
y2 = .25 + 1.25*x + rnorm(x)
tb = bind_rows(tibble(x=x, y=y1, z=F), tibble(x=x, y=y2, z=T))
注意，z是一个指示变量。我在R中计算了线性回归：
summary(lm(y ~ x + z, data=tb))
结果是
Call:
lm(formula = y ~ x + z, data = tb)
Residuals:
Min 1Q Median 3Q
浏览 0提问于2021-06-17得票数 1
我需要运行600多个回归，每个回归都在不同的MECE数据组上运行(组的值为{1,2，...,623})。从每个回归，我需要存储所有自变量的系数估计。我通过回溯(见下文)能够做到这一点；但是，我发现这个过程很慢，我相信还有一个更好的方法：
# loop prep
formula <- "dv ~ iv_1 + iv_2 + iv_3 | fe"
ols_stored_coef <- matrix(0, 623, 3)
ols_stored_coef <- as.data.frame(ols_stored_coef )
# loop
for(i in 1:623
浏览 6提问于2022-11-16得票数 2
回答已采纳
列表列数据结构是否是R语言中的一个概念？
什么是列表列？当将dataframe视为列列表时，列表列是否是dataframe中的列？
我从那里得到了困惑
列表列在数据帧的定义中是隐式的:数据帧是等长向量的命名列表。列表是向量，所以使用列表作为数据框架的列是合法的。但是，基本R不便于创建列表列，data.frame()将列表视为列列表：
data.frame(x = list(1:3，3:5) #> x.1.3 x.3.5 #> 1 1 3 #> 2 2 4 #> 3 3 5
您可以防止data.frame()使用I()进行此操作，但是结果打印得不太好：
data.fra
浏览 2提问于2020-06-03得票数 0
我有一个使用blmer创建的回归模型的列表，但是我无法提取p-值，因为它们不在anova(model)中，或者使用summary(model)$coefficients，因为我需要为多个系数提取p-值。只有当它是一个lmer模型时，我才可以提取p值列。从blmer的这些回归模型中是否有单独的函数或平均值来计算p值？
这里有一个例子，除了我有一个模型列表：
m1 = blme::blmer(Y ~ sex + age + (1|id/Group), data=df)
summary(m1)$coefficients
anova(m1)
我的输出不显示p值列，只显示t值，我知道lmer模型所显
浏览 1提问于2019-09-23得票数 1
回答已采纳
2回答
我需要对同一数据集上的许多线性模型进行估计，并将回归结果都放在一个表中。对于一个可重复的例子，下面是一个使用mtcars的简化
formula_1 = "mpg ~ disp"
formula_2 = "mpg ~ log(disp)"
formula_3 = "mpg ~ disp + hp"
目前，我的做法是：
创建一个包含所有公式的列表。
使用purrr:map()来估计所有的lm模型。
使用stargazer::生成输出表。
library(tidyverse)
library(stargazer)
formula
浏览 0提问于2019-07-12得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
