首页 文章
  • 4 votes
     answers
     views

    如何在没有模型对象的情况下从ns样条参数进行预测

    我有一个适合R的glm的系数,我想预测一组新数据的预期值 . 如果我有模型对象,这将很简单,使用predict() . 但是,我现在不在现场,出于数据机密性的原因,我不再拥有模型对象 . 我只有使用summary(model)生成的摘要对象,它包含模型系数 . 使用系数来预测简单模型的预期值很容易 . 但是,当模型包含三次样条ns()时,我想知道如何执行此操作 . 当模型还包括分类变量时的任何快捷...
  • 1 votes
     answers
     views

    如何在回归模型后预测R中的新数据?

    我想根据我创建的回归模型绘制并生成新数据: Trend <- seq(1:length(Daten$TSLAC)) TISPYNVDA_TSLA <- lm(TSLAC ~ Trend + SPYC * NVDAC, data = Daten) summary(TISPYNVDA_TSLA) plot(predict(TISPYNVDA_TSLA), type = "l...
  • 0 votes
     answers
     views

    R:绘制多元回归的预测结果

    我想观察一个治疗变量对我的结果Y的影响 . 我做了一个多元回归: fit <- lm (Y ~ x1 + x2 + x3) . x1 是治疗变量, x2 , x3 是控制变量 . 我使用了保持 x2 和 x3 的预测函数 . 我绘制了这个预测函数 . 现在我想在我的情节中添加一行类似于简单的回归 abline ,但我不知道该怎么做 . 我想我必须使用line(x,y),其中 y = p...
  • 3 votes
     answers
     views

    如何获得lmer对象的置信区间?

    我试图获得混合模型预测的置信区间 . 预测函数不输出任何置信区间 . 使用merTools包中的predictInterval函数建议的StackOverflow答案很少,以获得间隔,但这两个函数的预测估计值之间存在差异,我试图在下面的图中进行比较 . 有人可以告诉我这里我做错了什么吗?此外,我尝试构建的实际模型类似于下面的代码段中显示的模型,除了拦截之外我只有没有固定效果组件 . library...
  • 2 votes
     answers
     views

    通过交互进行GAM预测

    使用广义加法模型预测的常规是什么,包括与R库gam的交互? library("gam") x <- data.frame(a=runif(100,1,10), b=runif(100,1,10)) x$y <- x$a*x$b res <- gam(as.formula("y ~ s(a) + s(b)"), data=x[1:90,])...
  • 1 votes
     answers
     views

    回归 - 样本外预测

    我试图找出如何处理我的预测问题,我不确定我的理解是否在这个领域是正确的,所以如果有人可以帮助我真的很好 . 首先,我的目标是用回归预测时间序列 . 我没有使用ARIMA模型或其他启发式模型,而是专注于机器学习技术,如随机森林回归,k-最近邻回归等回归 . 以下是数据集的概述: Timestamp UsageCPU UsageMemory Indicator Dela...
  • 1 votes
     answers
     views

    glmer AICcmodavg的预测值

    我知道可以用AICcmodavg得到预测值(原始尺度〜概率)和SE用于固定效果,但是我没有成功地尝试......有人可以帮助我吗?提前致谢 library(lme4) (gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binom...
  • -3 votes
     answers
     views

    GLM后预测R [关闭]

    我训练了我的模型6000训练样本(glm)然后,我尝试预测200000行的矢量,但结果我只收到了6000行 . 我将这些参数用于函数predict(): predict( object = model_ppp2, newdate = Model_education, type = c("link", "response&...
  • 2 votes
     answers
     views

    R中的predict和glm.predict出错

    The Problem 我在R中训练了一个线性回归来预测数据框 data 中的变量 city . 此trainig是在数据子集上完成的,由 train.index 指定 . model = glm('data[, this.target] ~ data$city', data = data, subset = train.index) 我试图在保持数据上测试这个模型,该数据由 test.ind...
  • 2 votes
     answers
     views

    predict()错误:替换的行数多于数据

    当我使用 predict() 时,我真的可以't figure out why I'收到错误 . 我查了这篇文章,但我仍然得到同样的错误predict() . 我将数据框分成两部分(1. Train,2 . Test) . 我在火车上运行了一个逻辑模型并将其应用于测试,但是我遇到了错误 . 这是代码: train=rteam[which(rteam$season!="A"),]...
  • 0 votes
     answers
     views

    如何预测在R中使用glm?

    所以我有一个像这样定义的glm oring.glm = glm(oring.data$Damaged ~ oring.data$Temp, data = oring.data, family=binomial) 数据看起来像这样 Oring Temp 1 15 0 20 1 30 我想预测在特定温度下Oring会发生什么 我试过这样做 logodd...
  • 6 votes
     answers
     views

    如何预测随机和固定效应模型?

    我刚刚从STATA更改为R,并且在实现RATA等效的STATA命令 xtlogit,fe or re 和 predict 时遇到了一些麻烦 . 我可以请求一些帮助来调整以下场景: data <- read.table("http://people.stern.nyu.edu/wgreene/Econometrics/healthcare.csv",header=TRUE,...
  • 1 votes
     answers
     views

    glm`prepected()`错误:没有适用于'predict'的方法应用于类“list”的对象

    我在控制输入预测函数的对象类型时遇到问题 . 这是我生成 glm 对象的简化函数 . fitOneSample <- function(x,data,sampleSet) { #how big of a set are we going to analyze? Pick a number between 5,000 & 30,000, then select that many ...
  • -2 votes
     answers
     views

    从logit模型预测

    我正在运行以这种方式定义的logit模型: diversity_model <- glm(booking_bool ~ df$var_distance + df$var_price + df$var_prop_review_score + df$var_starrating + srch_hits + min_rating + ma...
  • 1 votes
     answers
     views

    让其他人看到e1071 :: predict.svm()概率的不连续性?

    我一直在R中使用e1071 :: svm(...,probability = TRUE)来拟合二进制SVM分类器,然后使用predict.svm()来获得训练样本和测试样本的概率 . 当我将概率转换为log(赔率)并将其与决策值进行对比时,我发现预测中存在不连续性: Plot of log(odds) = log(prob/(1-prob)) vs. Decision Values 只要概率低于0...
  • 0 votes
     answers
     views

    用R中的GAM包预测命令

    我试图在“gam”包中使用gam函数的预测命令 . 我发现当我有两个或更多样条线时,我无法预测新样本,即使新样本是原始数据或原始数据的子集 . 下面代码中的最后两行给出了错误 . 编辑:gam.data是gam包附带的数据集,用于gam和predict.gam的示例 . 它是一个带有列的data.frame:x,y,z,f,probf,ybin和ybin 2 - 全数字 . - require...
  • 0 votes
     answers
     views

    ggplot2图中的置信区间和使用R中的预测函数获得的值不相同

    这是我的数据(例子): EXAMPLE<-data.frame( X=c(99.6, 98.02, 96.43, 94.44, 92.06, 90.08, 87.3, 84.92, 82.14, 79.76, 76.98, 74.21, 71.03, 67.86, 65.08, 62.3, 59.92, 56.35, 52.38, 45.63, 41.67, 35.71, 30....
  • 1 votes
     answers
     views

    如何在R中用固定效应绘制Logit回归的预测概率?

    我是R的全新手 . 我估计有以下logit等式: allAM <- glm (AM ~ VS + Prom + LS_Exp + Sex + Age + Age2 + Jpart + X2004LS + X2009LS + X2014LS + factor(State), family = binomial(link = "logit"), data = mydata) ...
  • 1 votes
     answers
     views

    在lme4中模拟和预测之间的区别?

    lme4中的“模拟”和“预测”之间有什么区别(或者通常在混合效果模型中)? lme4文档只说: 模拟:模拟来自“merMod”拟合模型对象的响应 . 预测:预测值的数字向量 . lme4小插图说: 拟合方法检索以随机效应的所有模式为条件的拟合值;预测方法默认返回相同的值,但也允许以不同的随机效果集为条件进行预测 . 例如,如果re.form参数设置为NULL(默认值),则预测以模型中的所有随...
  • 1 votes
     answers
     views

    将csv导入OrderedDict并使用回归进行预测

    我 Build 了一个回归模型来预测来自5个变量(5列)的能量(1列)...我使用我的实验数据来训练和拟合模型并且它的得分很好...... import numpy as np import pandas as pd import matplotlib.pyplot as plt data = pd.read_csv('new.csv') X = data.drop(['E'],1) y = d...
  • 0 votes
     answers
     views

    H2O无人驾驶AI - 可视化和预测错误

    我的应用程序无人驾驶AI有问题 . 上传数据集没问题,但可视化或预测对我不起作用 . 附件中的错误消息 . 感谢您的解决方案 . predict error visualize error
  • 4 votes
     answers
     views

    svm预测错误

    我训练了一个svm模型 . 我想测试它,但我在predict()函数中遇到错误 . 为简单起见,我在这里将测试和训练数据拆分为非随机70/30分割 . library(e1071) train <- mydata[1:9731, ] test <- mydata[(9731+1):13901, ] mysvm <- svm(formula = outcome ~ BW +...
  • 0 votes
     answers
     views

    预测在R中使用多个变量

    我的R课程有点问题 . 我制作了一个以下数据集: 现在,我将使用以下命令基于此数据集绘制值: plot(x ~ Group.1, data = jarelmaks_vaikelaen23mean, xlab = "Vanus", ylab = "PD", main = "Järelmaks ja väikelaen") 之后...
  • 2 votes
     answers
     views

    用R中的plm函数预测

    我想知道是否可以使用来自R中的plm包的plm函数来预测新的预测变量数据集 . 我用以下方法创建了一个模型对象: model <- plm(formula, data, index, model = 'pooling') 现在我希望从新数据集中预测一个尚未用于估计模型的因变量 . 我可以通过使用模型对象中的系数来完成它,如下所示: col_idx <- c(...) df <-...
  • 0 votes
     answers
     views

    用回归量预测SVM模型

    我预测2016年两个月的电力负荷的SVM模型 . 训练数据集包含2012年至2016年的数据 . 我在我的数据集中添加了一组回归量,但是在预测它时会返回2012年的数据预测 - 2016年而不是2016年的第一个月 . 我做错了什么? 使用的代码: svm.model <- svm(NLtraints~NLxreg, data = NLtraints, type = "nu-reg...
  • 2 votes
     answers
     views

    关于使用张量流LSTM用时间序列数据预测下一个值的实际项目

    最近我研究了一个问题来预测第二天的网络视图 . 我选择了RNN-LSTM模型但效果并不理想 . 我的原始数据大约是200天,每天有1440点(每天有1440分钟,每分钟有一个值)=> [200 * 1440, 1] . 在特征工程之后,我将1-feature(仅为Web视图值)扩展为 8-features (当天的索引(范围从0到1439),Web视图值,is_weekday(0,1),...
  • 0 votes
     answers
     views

    保存和加载线性模型系数以进行预测的更好方法

    我正在努力为一年的数据集中的50000个客户制作每小时系数 . (365行* 28列) 我想保存这些系数,以便稍后在另一个R代码文件中进行预测 . 目前,我正在为客户使用保存功能保存24小时模型列表 . 所以,50000 Rda文件(每个7mb) . 然后,分别加载(加载函数)它们以使用R中的预测函数进行预测 . 这样效率不高,现在我想为一百万客户做这件事,这需要花费大量的时间和空间 . 有没有...
  • 0 votes
     answers
     views

    用不同的功能形式手动预测

    我有一个数据帧,系数来自glm(下面是 betas ) . 数据框包含协变量标签,协变量形式和估计值 . 形式是线性(Li),平方/二次(Sq)和log(Ps) . betas <- structure(list(CovGen = c("A", "B", "C", "D", "E", &quo...
  • 2 votes
     answers
     views

    在R中多次插补后计算预测平均值(或预测概率)和SE

    我想计算预测值和标准误差,但我不能简单地使用predict(),因为我正在使用15个多重插补数据集(生成Amelia包) . 我在每个数据集上运行回归模型 . 然后,使用使用鲁宾规则的Amelia函数mi.meld()将结果组合成一组模型系数和标准误差 . 示例数据和代码: dd<-list() for (i in 1:15){ dd[[i]] <- data.frame( Age...
  • 0 votes
     answers
     views

    R中的Predict()函数 . 如何使用它来预测因变量

    我有一个关于如何使用函数predict()的问题 . 我有一个包含n行和10列的数据集 . 第一列是因变量,其他变量是独立的 . 我在第一个变量上有50%的缺失数据,即x1,其他变量被完全观察到 . 我想通过使用以下模型中的对应案例和回归系数来预测x1(缺失部分): lm(new_A[,1]~new_A[,2]+new_A[,3]+new_A[,4]+new_A[,5]+new_A[,6]+new...

热门问题