It's not immediately evident to me just looking at this, but I do have
one suggestion. Run this again and see if you get the same error and
if you do, use this command to help you find where -
traceback()
traceback is extremely helpful when you are debugging these things, it
will go through your function and find where your function crashed if
not why exactly.
Give it a try.
A
--
Adrienne Wootten
Graduate Research Assistant
State Climate Office of North Carolina
Department of Marine, Earth and Atmospheric Sciences
North Carolina State University
On Tue, Jan 18, 2011 at 11:32 AM, wangwallace <talentt...@gmail.com> wrote:
>
> I was running a sampling syntax based on a data frame (ago) of 160 rows and
> 25 columns. Below are the column names:
>
>> names(ago)
> [1] "SubID" "AGR1" "AGR2" "AGR3" "AGR4" "AGR5" "AGR6" "AGR7"
> "AGR8"
> [10] "AGR9" "AGR10" "WAGR1" "WAGR2" "WAGR3" "WAGR4" "WAGR5" "WAGR6"
> "WAGR7"
> [19] "WAGR8" "WAGR9" "WAGR10" "ocbi" "ocbo" "cwbi" "cwbo"
>
>> dim(ago)
> [1] 160 25
>
> Here below is the syntax:
>
>> crossed1<-function(df){
> + s1<-df[sample(nrow(df),16,replace=F),]
> + s2<-t(apply(s1[,2:11],1,sample,10))
> + s2<-data.frame(s2)
> + s3<-df[-s1$SubID,]
> + ind<-t(replicate(nrow(s3),sample(10)))
> + s4<-matrix(nrow=nrow(s3),ncol=10)
> + rownames(s4)<-rownames(s3)
> + s4<-data.frame(s4)
> + for(i in seq(nrow(s3))){
> + s4[i,1]<-s3[i,1+ind[i,1]]
> + s4[i,2:10]<-s3[i,11+ind[i,2:10]]
> + }
> + colnames(s2)<-colnames(s4)
> + rbind(s2,s4)[order(as.numeric(rownames(rbind(s2,s4)))),]
> + }
>> result1<-vector("list",1000)
>> for(i in 1:1000)result1[[i]]<-crossed1(ago)
>
> These syntaxes worked out perfectly. I successfully drew 1000 random
> samples.
> I applied the syntax above to another data set called agr, which basically
> has the same column names:
>
>> names(agr)
> [1] "SubID" "AGR1" "AGR2" "AGR3" "AGR4" "AGR5" "AGR6" "AGR7"
> "AGR8"
> [10] "AGR9" "AGR10" "WAGR1" "WAGR2" "WAGR3" "WAGR4" "WAGR5" "WAGR6"
> "WAGR7"
> [19] "WAGR8" "WAGR9" "WAGR10" "ocbi" "ocbo" "cwbi" "cwbo"
>
> agr and ago are basically the same except that agr has more rows than ago:
>
>> dim(agr)
> [1] 300 25
>
> here below is the same syntax I've used for agr:
>
>> crossed1<-function(df){
> + s1<-df[sample(nrow(df),30,replace=F),] # note I have changed 16 to 30
> + s2<-t(apply(s1[,2:11],1,sample,10))
> + s2<-data.frame(s2)
> + s3<-df[-s1$SubID,]
> + ind<-t(replicate(nrow(s3),sample(10)))
> + s4<-matrix(nrow=nrow(s3),ncol=10)
> + rownames(s4)<-rownames(s3)
> + s4<-data.frame(s4)
> + for(i in seq(nrow(s3))){
> + s4[i,1]<-s3[i,1+ind[i,1]]
> + s4[i,2:10]<-s3[i,11+ind[i,2:10]]
> + }
> + colnames(s2)<-colnames(s4)
> + rbind(s2,s4)[order(as.numeric(rownames(rbind(s2,s4)))),]
> + }
>> result1<-vector("list",1000)
>> for(i in 1:1000)result1[[i]]<-crossed1(agr)
> Error in xj[i] : only 0's may be mixed with negative subscripts
>
> However, I got an error message shown above. It is kinda weird. Could
> anybody please explain what I've did wrong? Many thanks in advance!! :)
> --
> View this message in context:
> http://r.789695.n4.nabble.com/error-message-tp3223412p3223412.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
