楼主: 唐伯小猫
3188 12

[问答] 那位同学帮我解释一下这段代码?谢谢啦 [推广有奖]

  • 1关注
  • 16粉丝

VIP

学科带头人

12%

还不是VIP/贵宾

-

威望
0
论坛币
5371 个
通用积分
3.3476
学术水平
22 点
热心指数
28 点
信用等级
20 点
经验
36754 点
帖子
1427
精华
0
在线时间
1541 小时
注册时间
2005-10-3
最后登录
2025-9-18

楼主
唐伯小猫 发表于 2014-4-5 06:33:13 |AI写论文
30论坛币
附件1是 代码,2是代码运行完 结果

代码

的结果。因为刚刚接触R,所以很多不熟悉,恳请大家帮帮忙啊。谢谢啦

最佳答案

心若向阳,无畏悲伤。

沙发
lww1993 发表于 2014-4-5 06:33:14
RF<-randomForest(lat.max~.,data=dd,ntree=1000,importance=TRUE)
##建立随机森林,这是一个regression模型,生成1000棵分类树,importance=TRUE代表返回预测变量的重要性估计。

imp<-importance(RF)
##计算相应变量的重要性,值越大,代表当前的变量越重要;
##这里是回归模型,故第一个变量值如下计算:
##对于某个变量x,针对每棵树,用没有用到的值(out of bag)去计算MSE,然后置换out of bag样本中的x,
##计算MSE,这个时候MSE通常会变大,计算相应的差,对每棵树得到的差做平均,再除以其标准差,
##得到相应的值做为该变量的IncMSE
##第二个值如下计算:
##某个变量作为node时,计算其对node impurity的提高,我觉得是用的两个MSE的差。然后将所有值加起来,
##除以树的总棵数。

impvar<-rownames(imp)[order(imp[,1],decreasing=TRUE)]
##将预测变量名排序,根据其重要性的第一个指标值IncMSE,按从大到小排序;

op<-par(mfrow=c(3,3),mar=c(4,4,2,2))
##绘制图形,将整个界面等划分为3行3列,mar=c(4,4,2,2)以“下,左,上,右”的顺序指定边界大小,单位是英寸。
for(i in seq_along(impvar)){
  partialPlot(RF,dd,impvar[i],xlab=impvar[i],ylab="Longitude",ylim=c(26,30),main="")
}
##seq_along(impvar)相当于1:length(impvar)
##partialPlot中的参考文献对所画的函数描述不清,
##\tilde{f}(x) = \frac{1}{n} ∑_{i=1}^n f(x, x_{iC})
##我只是猜测,应该是固定除了关心的变量之外的变量,画出当前变量与响应变量之间的关系

藤椅
hugebear 发表于 2014-4-5 11:02:08
你把悬赏值当成代码售价?那谁还敢帮你看?

板凳
唐伯小猫 发表于 2014-4-5 14:27:03
我的天啊,昨天睡前晕了!!谢谢楼上的回复啊。马上更改
心若向阳,无畏悲伤。

报纸
唐伯小猫 发表于 2014-4-8 04:11:12
顶一下啦
心若向阳,无畏悲伤。

地板
唐伯小猫 发表于 2014-4-9 03:40:16
顶一下啦
心若向阳,无畏悲伤。

7
唐伯小猫 发表于 2014-4-9 03:41:34
顶一下啦
心若向阳,无畏悲伤。

8
windblood 发表于 2014-4-9 21:47:11
看不了附件图片

9
匡昊 发表于 2014-4-10 10:57:19
看不到代码

10
唐伯小猫 发表于 2014-4-10 13:55:45
上传,谢谢大家

Q2.png (54.03 KB)

Q2.png

Q1.png (63.38 KB)

Q1.png

心若向阳,无畏悲伤。

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-30 23:57