Hmisc:如何将reformM的结果作为参数包含在aregImpute中？
Hmisc是一个在R语言中广泛使用的包，它提供了许多数据处理和统计分析的函数和工具。reformM和aregImpute是Hmisc包中的两个函数。
reformM函数用于将模型的系数矩阵转换为一个向量，以便可以将其作为参数传递给其他函数。aregImpute函数是Hmisc包中用于多重插补的函数，它可以通过多个数据集的插补来处理缺失值。
要将reformM的结果作为参数包含在aregImpute中，可以按照以下步骤进行操作：
- 首先，使用reformM函数将模型的系数矩阵转换为一个向量。例如，假设模型的系数矩阵为coef_matrix，可以使用以下代码将其转换为向量：
coef_vector <- reformM(coef_matrix)- 接下来，使用aregImpute函数进行多重插补。将转换后的系数向量作为参数传递给aregImpute函数。例如，假设要处理的数据集为data，缺失值的变量为missing_var，可以使用以下代码进行多重插补：
imputed_data <- aregImpute(data, y = missing_var, coef = coef_vector)在上述代码中，data是要处理的数据集，y是包含缺失值的变量，coef是转换后的系数向量。
通过以上步骤，你可以将reformM的结果作为参数包含在aregImpute中，实现多重插补的功能。
关于Hmisc包的更多信息和使用方法，你可以参考腾讯云的R语言开发文档中关于Hmisc包的介绍：Hmisc包介绍。
相关·内容
1回答
我正在学习在R中实现健壮的glms，但不知道为什么当我有一个模型，其中一些列由于共线性而被丢弃时，我无法从回归模型中预测值。特别地，当我使用预测函数来预测来自glmrob的值时，它总是给出所有值的NA。在使用glm从相同的数据和模型预测值时，我没有观察到这一点。我使用什么数据似乎并不重要--只要拟合模型中有NA系数(并且NA不是系数向量中的最后一个系数)，预测就无效。
这种行为适用于我尝试过的所有数据集和模型，在这些数据集和模型中，由于共线性关系，删除了内部列。我包括一个假数据集，其中有两列从模型中删除，这在系数列表中给出了两个NAs。glm和glmrob都给出了几乎相同的系数，而预测只适用于
浏览 3提问于2017-05-23得票数 0
回答已采纳
2回答
小鼠多次归责后logistic和序数模型的AUC (附R)
r、imputation、auc、r-mice
我在问一个问题，关于把变量包含到逻辑模型和序数模型中的可加性预测的好处。我用老鼠来计算丢失的协变量，并且很难找到方法来计算集合归并模型的AUC和R平方。有人有什么建议吗？
摘要读出只提供术语，估计，std.error，统计，df，p.value。
示例代码：
imputed_Data <- mice(Cross_sectional, m=10, predictorMatrix=predM, seed=500, method = meth)
Imputedreferecemodel <- with(imputed_Data, glm(Poor ~ age + sex + educat
浏览 13提问于2022-11-22得票数 3
2回答
用R中的保存估计模拟logistic回归
r、regression、extract、simulation、logistic-regression
我有点问题。我正在尝试开发一些代码，使我可以这样做: 1)运行一个logistic回归分析，2)从logistic回归分析中提取估计，3)使用这些估计创建另一个logistic回归公式，我可以在以后的模拟中使用原始模型。因为我是相对较新的R，我知道我可以通过索引来提取这些系数1×1，但很难将其“缩放”到具有不同系数数的模型中。我想知道是否有更好的方法来提取系数和建立公式。然后，我必须开发实际的变量，但是这些变量的开发必须足够灵活，以适应任意数量的变量和分布。这在Mplus中似乎很容易做到( Mplus手册中的示例12.7 )，但我还没有在R中找到这一点。据我所知，下面是代码：
#generat
浏览 12提问于2022-10-21得票数 2
回答已采纳
1回答
pool.compare生成不符合条件的参数错误。
r、r-mice
交替标题:模型矩阵和系数集显示不同的变量数。
我正在使用小鼠软件包R做一些分析。我想使用pool.compare()比较两种模型(保存在米拉对象中)，但我始终得到以下错误：
Error in model.matrix(formula, data) %*% coefs : non-conformable arguments
二进制运算符%*%表示R中的矩阵乘法。
表达式model.matrix(formula, data)生成“具有指定公式和数据的回归模型的设计矩阵”(来自model.matrix {stats}的R文档)。
在错误消息中，coefs是从est1$qbar中提取的，其中est1
浏览 3提问于2016-04-19得票数 1
3回答
从R中的lm中提取标准化系数
r、regression、lm、beta、standardized
我为愚蠢的question...but道歉我似乎找不到一个简单的解决方案
我想从一个拟合的线性模型中提取标准化系数(在R中)，必须有一个简单的方法或函数来做到这一点。你能告诉我这是什么吗？
编辑(在下面的一些评论之后)：我应该为我的问题提供更多的上下文信息。我在教一群心理学家的入门R工作坊。对他们来说，没有能力获得标准化系数的线性模型就像你根本没有运行模型一样(好吧，这有点夸张，但你明白了)。当我们做了一些回归后，这是他们的第一个问题，(我的错)我没有预料到(我不是心理学家)。当然我可以自己编程，当然我也可以寻找为我做这件事的包。但与此同时，我确实认为这是线性模型的一种基本和常见的必需功能，在
浏览 0提问于2014-06-19得票数 50
回答已采纳
2回答
如何求R中gradient的梯度
r、plot、line、gradient
我想在R中的散点图上获得线性模型的梯度/斜率的值
例如，我有一个图：
x <- 1:10
y <- c(0, 0, 1, 0, 1, 1, 0, 1, 1, 1)
plot(x, y)
abline(lm(y ~ x))
散点图现在有一条线，我需要知道它的梯度。理想情况下，我想要的是一个可以计算出各种x/y组合的梯度的函数，而不需要每次都绘制图形。我可以自己创建函数，但我不知道如何获得梯度。
浏览 2提问于2013-10-29得票数 3
回答已采纳
1回答
我想用测试我的模型的重要性。我读过这篇文章，因为我使用的是集群-健壮的模型，F测试不支持，我应该使用Wald测试。
我的脚本目前看起来是这样的，所有这些不同的选项都给了我更正的集群健壮se：
选项1:
m <- lm(y_var ~ var1 + poly(var2, 2) + quartier, data = df)
m_robust_clustered <- coeftest(m, vcov = vcovCL,
type = "HC1",
浏览 0提问于2021-11-12得票数 2
3回答
R中任意系数的预测()
r、glm、predict
我得到了由非r用户设置的logit模型的一些系数。我希望将这些系数导入r中，并在同一数据集(ROC和混淆矩阵)上生成一些拟合优度估计值，以对抗我自己的模型。我的第一个想法是将系数强制放入现有的GLM对象中，方法如下
摘要(Fit)$coefficients，1 <- y
或
摘要(Fit)$coefficients <- x
其中y和x是包含我试图用来预测和拟合的系数的矩阵，是以前创建的适合于数据集的虚拟glm对象。当然，这只会给我带来错误。
是否有任何方法将任意系数向量传递给预测()函数或指定模型中的系数？我可以通过将向量传递到GLM中的偏移量参数来强制这样做吗？
浏览 3提问于2014-09-06得票数 7
我刚刚在conda中安装了Rstudio和R，我正在尝试遵循本教程：
我一直跑到这一步
> completeData <- complete(imputed_Data,2)
我发现了一个错误：
(函数(类、fdef、mtable)中的错误:无法找到用于签名的函数‘完全’的继承方法“mids”、“数字”
它与我如何安装R有关吗？你知道怎么修吗？谢谢。
更多信息
> R.Version()$version.string
[1] "R version 3.4.3 (2017-11-30)"
> packageVersion("mice"
浏览 0提问于2018-03-23得票数 2
回答已采纳
1回答
混合glm零膨胀模型的Bootstrap方法
r、mixed-models、glmmtmb
我想使用m_F包引导混合glm零膨胀模型( glmmTMB )，但是尽管系数规范使用了coef或fixef，但我总是输出错误：
Error in bres[i, ] <- coef(bfit) :
incorrect number of subscripts on matrix
我的例子：
library(glmmTMB)
library(boot)
my.ds <- read.csv("https://raw.githubusercontent.com/Leprechault/trash/main/ds.desenvol.csv")
str(my.ds)
#
浏览 9提问于2022-07-05得票数 2
回答已采纳
1回答
获取两组数据之间所有可能的相关性
r、correlation
我正在尝试查找每个关联之间的几个信息 corr.test 我有两个数据集df1和df2 df1<- structure(list(col1A = c(1.64, 0.03, 0, 4.202, 2.981, 0.055,
0, 0.002, 0.005, 0, 0.002, 0.649, 2.55, 2.762, 6.402), col2A = c(2.635,
0.019, 0, 5.638, 3.542, 0.793, 0.259, 0, 0.046, 0.004, 0.017,
0.971, 3.81, 3.104, 5.849), col3A = c(0.91, 0.03
浏览 19提问于2019-01-30得票数 0
1回答
R中ARIMA输出的定义
r、definition、forecasting
我被要求尝试从一个ARIMA模型(Φ，θ)中找到参数，我最初以为它只是ARIMA模型的顺序，例如(1,0,1)(1,1,0)，但是我的经理说这不是他想要的。
我看了ARIMA模型的总结，在运行之后，我找到了phi，theta和delta的值，但是有不止一个值，我很不确定哪一个是我需要的。
下面是我如何找到这些值的方法：
M<-matrix(c("08Q1", "08Q2", "08Q3", "08Q4", "09Q1", "09Q2", "09Q3", "09Q4
浏览 2提问于2014-10-08得票数 2
回答已采纳
1回答
我试图在R中设置一个计算多项式的函数
P(x) = c1 + c2*x + c3*x^2 +.+ cn-1*x^n-2 + cn*x^n-1
对于x和集合系数c的各种值。
霍纳的方法是
集cn =10亿
对于i= n-1，n-1，…，2，1，集合bi = bi+1*x + ci
返回输出
到目前为止我所拥有的是：
hornerpoly1 <- function(x, coef, output = tail(coef,n=1), exp = seq_along(coef)-1) {
for(i in 1:tail(exp,n=1)) {
(o
浏览 2提问于2018-03-21得票数 0
回答已采纳
1回答
R中wald.test的传递假设
r、lag、lm、hypothesis-test
请原谅，因为我是这个论坛的新手。这项研究需要检查coefficients=0之和。测试可以使用eview进行，如c(2)+c(3)+c(4)=0，其中2是第二项的系数，因此是第四项的系数。相同使用R的代码是
require(Hmisc)#this package is used to generate lags
require(aod)#this package is used to conduct wald test
output<-lm(formula = s_dep ~ m_dep + Lag(m_dep,-1) + Lag(m_dep,-2) + s_rt
浏览 2提问于2015-09-21得票数 2
回答已采纳
1回答
为什么删除包含R中因素的Anova中的变量并不会降低自由度？
r、regression、anova、r-factor
我试图在不同的模型之间进行一次方差分析( ANOVA和线性回归之间的混合)，我遇到了一些问题。我想我把范围缩小到了关于方差分析的一个问题(或一些我不理解或做错的事情)：为了对两个模型进行比较，它们需要有不同的剩余Df (自由度)。
例如，让我们考虑R中的mtcar数据：
library(car)
test_data <- mtcars %>% mutate(factored_variable = as.factor(carb))
model_1 <- aov(drat ~ factored_variable , data = test_data)
Anova(model_1
浏览 3提问于2017-06-20得票数 3
回答已采纳
1回答
如何在svy_vglm中获取分类变量的全局p
r、survey、multinomial
我使用函数svyVGAM::svy_vglm运行一个带有调查权重的多项式模型： mmodel <- svy_glm(y~x1+x2+x3+x4..., family=multinomial, design=w_data) 其中x表示分类变量，其中一些变量具有三个或更多级别。通过模型汇总，我可以知道每个系数的p值，但我不知道如何获得全局变量的p值。 在其他上下文中，anova()，waldtest()，lrtest()，...可以使用，但它们似乎都不能与svy_vglm对象一起使用。Error: No tidy method for objects of class svy_vglm：t
浏览 49提问于2021-03-20得票数 1
回答已采纳
1回答
如何在texreg表中包含赔率比？
r、logistic-regression、texreg
所以我尝试在我的texreg表中显示赔率比：
library(texreg)
df <- mtcars
logit <- glm(
formula = am ~ cyl,
data = df,
family = binomial(link = "logit")
)
odds_ratio <- exp(coef(logit))
screenreg(logit, odds_ratio)
不幸的是，这会产生以下错误：
Error in screenreg(logit, odds_ratio) :
The 'file' arg
浏览 17提问于2020-05-14得票数 0
1回答
dplyr软件包中的多元线性回归
r、dplyr、regression、coefficients
假设我有20辆车，每辆车被测试了10次，在每一次测试中，有多个段是根据路径距离来定义的。
我想得到简单的线性回归斜率的回归对高度与路线的距离，为每个veh_id+test_id+segment组合。
在R中使用dplyr包，我认为这样的东西应该能起作用：
dataset <- dataset %>%
group_by (veh_id, test_id, segment) %>%
do(data.frame(mod = coef(lm(height ~ distance,data = .))[2])) %>%
left_join(dataset,.)
但是，错误报告为：
浏览 0提问于2017-06-15得票数 0
3回答
用相同的设计矩阵拟合R中的许多线性模型
r、parallel-processing、regression、mlm
对于神经成像应用程序，我正在尝试通过R(对lm的标准调用)中的最小二乘来拟合许多线性模型。假设我有一个设计矩阵X，这个设计矩阵在所有模型中都是相同的。正在拟合的数据(Y)将改变，因此所有拟合参数(例如，betas、p值、残差等)也将改变。
目前，我只是把它放在一个for循环中，所以它正在对lm进行数十万次调用。看起来肯定有更好的办法。
我认为计算量最大的部分是矩阵求逆。看起来这是通过lm.fit中的Fortran调用来处理的。
如果我手工做这个回归，我会做矩阵求逆，然后把它乘以各种数据集。事实上，当我有行为良好的设计矩阵(例如，所有连续值协变量)时，我已经编写了一个函数来实现这一点。然而，我真
浏览 2提问于2013-01-29得票数 2
回答已采纳
2回答
我已经拟合了一个累积链接模型，使用了R中的clm包中的ordinal函数，我需要从对象那里得到一个系数估计，以便进一步计算。所讨论的值是来自下面gist：中的表的
(对不起，我没有在这里输出表，但由于某种原因，我无法正确地格式化它。)
我使用str函数检查了对象，但找不到我要找的东西。
浏览 3提问于2015-02-22得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
