Bootstrap   发布时间:2022-04-18  发布网站:大佬教程  code.js-code.com
大佬教程收集整理的这篇文章主要介绍了R中的Bootstrap多项式回归大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。
我试图在R中引导一个简单的多项式回归,我收到一个错误

@H_301_11@

Error in is.data.frame(data) : object ‘d’ not found

真正奇怪的是,我使用相同的代码(调整到这个特定的问题),就像在启动包at Quick-R的教程中一样,当我使用不同的函数(如lm())时,相同的代码也有效.当然,我做了一些愚蠢的事情,但我看不出是什么.如果有人可以提供帮助,我会非常感激.

这是一个例子:

require(foreign)
require(nnet)
require(boot)

# an example for multinomial logistic regression
ml = read.dta('http://www.ats.ucla.edu/stat/data/hsbdemo.dta')
ml = ml[,c(5,7,3)]

bs <- function(formula,data,indices) {
    d = data[indices,] # allows boot to select sample
    fit = multinom(formula,data=d)
    s = summary(fit)
    return(list(fit$coefficients,fit$standard.errors))
}

# 5 replications
results = list()
results <- boot(
    data=ml,statistic=bs,R=5,parallel='multicore',formula=prog~write
)

解决方法

错误发生在summary()部分,multinom()返回的对象也没有系数和standard.errors.似乎,summary.multinom()依次计算你的数据中的hessian d,由于某种原因(可能是一个范围问题)无法找到.快速解决方法添加Hess = TRUE:

bs <- function(formula,indices) {
  d = data[indices,] # allows boot to select sample
  fit = multinom(formula,data=d,Hess = TRUE)
  s = summary(fit)
  return( cbind(s$coefficients,s$standard.errors) )
}

# 5 replications
results = list()
results <- boot(
  data=ml,formula=prog~write
)

大佬总结

以上是大佬教程为你收集整理的R中的Bootstrap多项式回归全部内容,希望文章能够帮你解决R中的Bootstrap多项式回归所遇到的程序开发问题。

如果觉得大佬教程网站内容还不错,欢迎将大佬教程推荐给程序员好友。

本图文内容来源于网友网络收集整理提供,作为学习参考使用,版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ:384754419,请注明来意。
标签: