Checking on the function 'combinations' in gplot i find this: if (mode(n) != "numeric" || length(n) != 1 || n < 1 || (n%%1) != 0) stop("bad value of n")
So essentially, the argument 'length(data$times[data$ids == id[1]])' in the line 'combinations(length(data$times[data$ids == id[1]]), 2)' does not meet one of those four conditions. That's as much as i can tell you without some kind of reproducible code. 2010/9/16 trekvana <trekv...@aol.com> > > Hi all- > > this seems to be simple to figure out but since im new to writing functions > I dont know what is happening. Here is my code along with the error I am > receiving: > > semivario=function(data,ids,times,resids){ > id=unique(data$ids) > index=combinations(length(data$times[data$ids==id[1]]),2) > time=gamma=numeric(dim(index)[1]) > for (j in 1:dim(index)[1]){ > > > time[j]=abs(data$times[data$ids==id[1]][index[j,1]]-data$times[data$ids==id[1]][index[j,2]]) > > > gamma[j]=.5*(data$resids[data$ids==id[1]][index[j,1]]-data$resids[data$ids==id[1]][index[j,2]])^2 > } > > for (i in id[-1]){ > index=combinations(length(data$times[data$ids==i]),2) > time1=gamma1=numeric(dim(index)[1]) > for (j in 1:dim(index)[1]){ > > > time1[j]=abs(data$times[data$ids==i][index[j,1]]-data$times[data$ids==i][index[j,2]]) > > > gamma1[j]=.5*(data$resids[data$ids==i][index[j,1]]-data$resids[data$ids==i][index[j,2]])^2 > } > time=c(time,time1) > gamma=c(gamma,gamma1) > } > value=list(time=time,gamma=gamma) > return(value) > } > > and I would run it as > semivario(data=hope,ids=id,times=sctime,resids=scResid) > > What I want to do is everywhere we see the words data,ids,times,resids to > be > replaced with, for example, hope,id,sctime,scResid. Think of macros in SAS > where we use the & before every macro input variable name. Im not sure if > this is how R handles function inputs. > > Here is my error: > Error in combinations(length(data$times[data$ids == id[1]]), 2) : > bad value of n > > > Thanks > > -- > View this message in context: > http://r.789695.n4.nabble.com/Pesky-homemade-function-code-tp2541597p2541597.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. > [[alternative HTML version deleted]] ______________________________________________ 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.