cat函数用于输出到外部文件:
cat("av3=function(x){n=length(x)m=mean(x)d3=(x-m)^3sum(d3)/n} ",file="aa.r")
cat("任意字符",file=“文件名”)
之后要运行这个R就写一句source("文件名")就好了
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中