请教for循环语句的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

sdata=read.table("D:/sdata.txt",header=T)
attach(sdata)

num=5
library(splines)
require(stats); require(graphics)
basis=bs(sdata$t, df=num)#产生5个B样条基

for (i in 1:num)
bs.i=basis[,i]

bs11=NULL
for (i in 1:num)
bs11=cbind(bs11,bs.i)    #请问这个循环哪里有问题，谢谢！

[/code]

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：for循环 循环语句 For Graphics require 请教 语句

你第一个for循环完了之后，bs.i等于basis[, 5]
然后第二个循环把它重复了5次。这应该不是你想要的结果吧。

我不懂B-spline，但是感觉你的两个for循环应该合并为
bs11<- NULL
for (i in 1:num)
    bs11 <- cbind(bs11, basis[,i])

哦？！不过不明白为何楼主这么写？既然bs11 = NULL，为何不干脆只用basis，还需要添加NULL到cbind去？得出的结果 bs11不就等于basis了吗？

1. >bs11 <- NULL
   
2. >for (i in 1:5) cbind(bs11,basis[,i])

复制代码

谢谢啊，现在可以了。我本想先看看bs.i，然后再把他们合在一起的，因为后面要分别用到这两部分。另外问一下，我有几个因变量，假设有3个，自变量是一样的，如何利用循环把得到的预测结果放在一起？

1. bs就是前面的结果
   
2. ww11=data.frame(cbind(ssq$s.1,ssq$s.2,ssq$s.3,ssq$s.4,ssq$s.5,ssq$s.6,ssq$s.7,bs11))#s.1,s.2,s.3,s.4,s.5,s.6,s.7怎么简写更方便呢？
   
3. library(pls)
   
4. pred=NULL
   
5. plsmod=NULL
   
6. 
   
7. for ( i in 1: ncol(sdata))
   
8. plsmod[i]=mvr(ww11[,i]~bs11, ncomp =num, data = ww11, scale = TRUE,method = "oscorespls")###问题就出在如何保存plsmod上？
   
9. validationplot(plsmod[i,])
   
10. pred[i,]=predict(plsmod[i,], comps = 1:num, newdata = predictdata) #预测。

复制代码

statchao 发表于 2011-5-14 16:14
谢谢啊，现在可以了。我本想先看看bs.i，然后再把他们合在一起的，因为后面要分别用到这两部分。另外问一下，我有几个因变量，假设有3个，自变量是一样的，如何利用循环把得到的预测结果放在一起？

1. bs就是前面的结果
   
2. ww11=data.frame(cbind(ssq$s.1,ssq$s.2,ssq$s.3,ssq$s.4,ssq$s.5,ssq$s.6,ssq$s.7,bs11))#s.1,s.2,s.3,s.4,s.5,s.6,s.7怎么简写更方便呢？
   
3. library(pls)
   
4. pred=NULL
   
5. plsmod=NULL
   
6. 
   
7. for ( i in 1: ncol(sdata))
   
8. plsmod=mvr(ww11[,i]~bs11, ncomp =num, data = ww11, scale = TRUE,method = "oscorespls")###问题就出在如何保存plsmod上？
   
9. validationplot(plsmod[i,])
   
10. pred[i,]=predict(plsmod[i,], comps = 1:num, newdata = predictdata) #预测。

复制代码

没有提供样本数据比较难帮到你， 只能建议试一试 lapply...

1. num = 1:5
   
2. ww11 <- structure(data.frame(lapply(num, function(i) ssq[[i]]$s)),names = names(ssq)) #由于不晓得ssq是list还是data.frame，所以...
   

复制代码

谢谢啊，s.1+s.2+s.3+s.4+s.4+s.5有没有简写的方式啊？

假设

1. s=cbind(s.1,s.2,s.3,s.4,s.5)

复制代码

则可以简写为

1. rowSums(s)

复制代码

1. good!

复制代码

我想运行x<- c(1:60)
for(i in 1:100){y[i]<- x+ rnorm(60,0,15)
i<- i+1}
想产生100组y，但是提示错误