curl_setopt($curl,CURLOPT_URL,$_GET['url'])
curl_setopt($curl,CURLOPT_RETURNTRANSFER,true)
curl_setopt($curl,CURLOPT_USERAGENT,"Mozilla/4.0
(compatibleMSIE 6.0Windows NT 5.1SV1.NET CLR 1.1.4322.NET CLR 2.0.50727)")
$output = curl_exec($curl)
curl_close($curl)
print_r($output)
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中
