用R语言模拟weibull分布，后面总是提示ADR1、d1等，请哪位高手指点，不胜感激。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

> d1<-read.table("F://d1.txt",header=TRUE)
> d2<-read.table("F://d2.txt",header=TRUE)
> ds1<-data.frame(x=d1$x,y=d1$y)
> ds2<-data.frame(x=d2$x,y=d2$y)
> m1<-nls(y~I(exp(-Scale*x^Shape)),data=ds1,start=list(Scale=0.01,Shape=1),trace=T)
13.42884 :  0.01 1.00
3.991207 :  0.08704032 0.75903022
0.500503 :  0.02950857 1.56973971
0.07841429 :  0.03557713 1.62886837
0.07373831 :  0.03032661 1.71335289
0.0730652 :  0.03020684 1.72198921
0.07306428 :  0.03013387 1.72329074
0.07306428 :  0.03012837 1.72339020
0.07306428 :  0.03012794 1.72339785
> m2<-nls(y~I(exp(-Scale*x^Shape)),data=ds2,start=list(Scale=0.01,Shape=1),trace=T)
13.42884 :  0.01 1.00
3.991207 :  0.08704032 0.75903022
0.500503 :  0.02950857 1.56973971
0.07841429 :  0.03557713 1.62886837
0.07373831 :  0.03032661 1.71335289
0.0730652 :  0.03020684 1.72198921
0.07306428 :  0.03013387 1.72329074
0.07306428 :  0.03012837 1.72339020
0.07306428 :  0.03012794 1.72339785
> plot(d1$x,d1$y,pch=20,xlim=c(0,170),ylim=c(0,1),xaxs="i",yaxs="i",las=1,mgp=c(1.7,0.5,0),lab=c(10,10,12),xlab="Survival(cycle)",ylab="Survival Probablity",font.lab=2)
> points(d2$x,d2$y,pch=17)
> s<-seq(0,170,length=200)
> lines(s,predict(m1,list(x=s)),lty=2)
> lines(s,predict(m2,list(x=s)),lty=1)
> ADR1<-1-(sum(residuals(m1)^2)/(length(d1$x)-2))／(sum((dl$y-mean(dl$y))^2)／(1ength(dl$x)-1))
> ADR2<-1-(sum(residuals(m2)^2)/(length(d2$x)-2))／(sum((d2$y-mean(d2$y))^2)／(1ength(d2$x)-1))
> ADR<-c(ADR1,ADR2)
错误: 找不到对象'ADR1'
> Cor.eff<-c(cor(dl$x,d1$y,method="pearson"),cor(d2$x,d2$y,method="spearman"))
错误于inherits(x, "data.frame") : 找不到对象'dl'
>

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Weibull分布 Weibull 不胜感激 高手指点 Bull 不胜感激