求助：怎样求门限值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

高手们，在Splus软件中，Peak-over-Threshold（POT）方法中，什么命令可以求门限值？怎样可以得到广义帕累托分布的参数？万分感激！！！！！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：门限值 Threshold 广义帕累托分布 帕累托分布 hold 求助 限值

library(segmented)
olm<-lm(walevel~day)

#specify number and starting values for the breakpoints..
oseg<-segmented(olm, seg.Z=~day, psi=c(20,50,80))
plot(oseg,add=TRUE,col=2)
oseg$psi

#or you can avoid to specify starting values for psi
oseg<-segmented(olm, seg.Z=~day,
psi=NA,control=seg.control(stop.if.error=FALSE,K=30))

plot(oseg,add=TRUE,col=2)
oseg$psi