找回密码
 立即注册
查看: 5120|回复: 0

求大神指导呀?如何实现散点图,拟合曲线以及置信线在同一张图里

[复制链接]
发表于 2014-4-24 09:18:50 | 显示全部楼层 |阅读模式
如何实现散点图,拟合曲线以及置信线在同一张图里,跪求大神指导接下来的代码如何写呀

boot<-function(x,pr,fboot,B)
{
    Tboot=0
    for(i in 1:B)
     {   Y=sample(x,20,prob=pr,replace=TRUE)
         Tboot[i]=fboot(Y)
     }
  Tboot
}
Data=matrix(0,ncol=20,nrow=10)

  Data[1,]=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1)
  Data[2,]=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,2)
  Data[3,]=c(0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,2,2,2,2,3)
  Data[4,]=c(0,0,0,0,0,0,0,1,1,1,2,2,2,2,2,3,3,4,4,5)
  Data[5,]=c(0,0,0,1,1,2,2,3,3,3,3,3,4,4,4,4,5,6,7,7)
  Data[6,]=c(1,2,2,2,3,3,4,4,4,4,4,5,5,5,5,5,6,6,7,8)
  Data[7,]=c(3,3,3,5,5,6,6,7,7,7,7,8,8,8,8,9,9,9,12,13)
  Data[8,]=c(3,4,6,7,8,8,8,8,9,9,9,10,10,10,11,11,11,12,13,16)
  Data[9,]=c(5,6,7,8,8,9,9,9,9,10,10,10,10,11,11,12,12,13,16,16)
  Data[10,]=c(5,6,6,7,8,8,9,10,11,11,11,11,12,12,13,13,13,15,16,18)

T=matrix(0,ncol=10000,nrow=10)
for(i in 1:10)
   { x=Data[i,]

pr=c(1/20,1/20,1/20,1/20,1/20,1/20,
         1/20,1/20,1/20,1/20,1/20,1/20,
         1/20,1/20,1/20,1/20,1/20,1/20,1/20,1/20)

  tb=boot(x,pr,mean,10000)

    T[i,]=tb  
   }

A=matrix(0,ncol=10000,nrow=10)

A[1,]=T[1,]/0.5 A[2,]=T[2,]/1 A[3,]=T[3,]/2 A[4,]=T[4,]/4 A[5,]=T[5,]/6 A[6,]=T[6,]/8 A[7,]=T[7,]/10 A[8,]=T[8,]/12 A[9,]=T[9,]/14 A[10,]=T[10,]/16 x<-c(T[1,],T[2,],T[3,],T[4,],T[5,],T[6,],T[7,],T[8,],T[9,],T[10,]) y<-c(A[1,],A[2,],A[3,],A[4,],A[5,],A[6,],A[7,],A[8,],A[9,],A[10,]) fit<-lm(y~x+I(x^2))



回复

使用道具 举报

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

本版积分规则

Archiver|手机版|小黑屋|R语言中文网

GMT+8, 2024-11-25 10:07 , Processed in 0.019247 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表