全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2250 9
2021-10-05
本人小白,没接触过编程,大佬们能讲的细一点么?这是我瞎写的代码,也能运行
A <-NULL###结果矩阵
for (i in 1:4263) {###循环4263次
        REI <- mydata[,c(1,24)]###从整个数据集中提取需要的两列
        REI <- REI[1:1000+i,]###分别取数据集(两列的那个)的前1001行、1002行...直至全部数据集(一共5263行)
        a <- fitgpd(REI$实体指数损失率...,threshold = 5)$pat###运行一个函数,分别取出里面的一个参数
        A <- rbind(A,a)###将每次的参数作为一条数据追加到结果矩阵中
}


但是出了一些问题,按理说A矩阵里的最后一个数不应该是全部数据集运行函数后的参数么?
A矩阵最后一个数是0.024,如下
结果矩阵
但我将全部数据运行后,得到的参数是0.0188,如下
全部数据运行函数结果
这是咋回事啊?我该咋整啊?感谢大佬们指点
附件列表

实体总指数.xlsx

大小:12.58 MB

 马上下载

全部数据集

实体总指数(损失率).xlsx

大小:1.07 MB

 马上下载

全部数据

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2021-10-5 17:55:32
复制代码
我用自己编的数据是没问题的,如何还不行,请把配图把问题描述清楚一些,并提供一些数据
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2021-10-5 19:37:17
s609078902 发表于 2021-10-5 17:55
我用自己编的数据是没问题的,如何还不行,请把配图把问题描述清楚一些,并提供一些数据
感谢大佬抽空解答,但是还是不行,数据集跟图上传了,您的方法也试了,也不是我想要的结果,甚至跟全部数据集运行函数后的param值对不上。这个函数在POT包里,要不麻烦您运行下?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2021-10-6 10:04:05
Aixir 发表于 2021-10-5 19:37
感谢大佬抽空解答,但是还是不行,数据集跟图上传了,您的方法也试了,也不是我想要的结果,甚至跟全部数 ...
没看到有实体指数损失率这个变量呀
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2021-10-6 10:51:14
Aixir 发表于 2021-10-5 19:37
感谢大佬抽空解答,但是还是不行,数据集跟图上传了,您的方法也试了,也不是我想要的结果,甚至跟全部数 ...
复制代码
1. 问题出在 1:1000+i 的运算优先级上,写成 1:(1000+i) 这样就没问题了;
2. 执行 for 循环之前应该提前留好足够的空间储存结果,像我这里这样,可以提高执行速度;
3. 你注释给得太多了,反而会影响阅读代码的速度;
4. 赋值 REI 的两行语句可以合并成一行,也能提高点运行速度。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2021-10-6 10:58:44
Aixir 发表于 2021-10-5 19:37
感谢大佬抽空解答,但是还是不行,数据集跟图上传了,您的方法也试了,也不是我想要的结果,甚至跟全部数 ...
复制代码
1. 问题出在 1:1000+i 的运算优先级上,写成 1:(1000+i) 就没问题了;
2. 执行 for 循环之前应该预留足够的空间给结果, 像我这里这样,可以提高运行速度;
3. 你的注释太多了,会影响阅读代码的速度;
4. 复制 REI 的两行语句可以合并成一行,也可以提高运行速度。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群