R语言保存循环结果

Python032

R语言保存循环结果,第1张

k <- list()

for(i in 1:1000)

{

  k[[i]] <- nn2()

}

newdata=c()                        #1

for(i in 1:1000)

{

#方法一：三次样条法

library(splines)

m1 <- lm(h~bs(a,df=3),data=k[[i]])

#预测百分位数值

new <- data.frame(a=7:20)

cs.p <- predict(m1, new)

#均方差

mse.cs <- sum( (st$p50-cs.p)^2 )/14

#最大范数误差

mne.cs <- max(abs(st$p50-cs.p))

newdata<-rbind(newdata,mse.cs)        #2

print(newdata)                        #3

}

aa<-mean(newdata)          #4

新建newdata来保存循环的结果，以便对循环的结果进行后续操作比如求均值并保存在aa中

png("绘制图形.png")

{

multiplot绘图语句；

}

dev.off()

运行这段代码，就可以在当前目录下，生成一个PNG图形。

以下为示例：

a<-rnorm(50)

a<-sort(a)

pdf("d:/demo.pdf")

plot(a,type="b",col="red")

hist(rnorm(40))

dev.off()