Hi,

I am almost getting there, but still have errors. Thanks for your help. 
I have tried improving but I get the following errors below:

>itn<-function(it){
+siml<-matrix(NA,ncol=5,nrow=it)
+for(g in 1:it){
+siml[g]<-optm(perm=20)[g]
+}
+siml
+}
>itn(3)
              [,1] [,2] [,3] [,4] [,5]
[1,] 0.8873775898   NA   NA   NA   NA
[2,] 0.0015584824   NA   NA   NA   NA
[3,] 0.0001414317   NA   NA   NA   NA


>itn<-function(it){
+siml<-matrix(NA,ncol=5,nrow=it)
+for(g in 1:it){
+siml[g]<-optm(perm=20)
+}
+siml
+}
>itn(3)
           [,1] [,2] [,3] [,4] [,5]
[1,] 0.8880941   NA   NA   NA   NA
[2,] 0.8869727   NA   NA   NA   NA
[3,] 0.8877045   NA   NA   NA   NA
Warning messages:
1: In siml[g] <- optm(perm = 20) :
   number of items to replace is not a multiple of replacement length
2: In siml[g] <- optm(perm = 20) :
   number of items to replace is not a multiple of replacement length
3: In siml[g] <- optm(perm = 20) :
   number of items to replace is not a multiple of replacement length


I expect something close to

         average           sd          se       min       max
         0.8881969 0.0008215379 0.000410769 0.8873842 0.8890167

        0.884659 0.0004215379 0.000410769 0.2342 0.676307

        0.8885839 0.0001215379 0.0002112 0.000082752992 0.8836337




Thanks fpr you help.


On 6/7/2013 5:24 AM, Berend Hasselman wrote:
> On 07-06-2013, at 10:59, Laz <lmra...@ufl.edu> wrote:
>
>> Dear R users,
>>
>> I am stuck here: My first function returns a vector of 5 values.
>> In my second function, I want to repeat this, a number of times, say 10
>> so that I have 10 rows and five columns but I keep on getting errors.
>>
>> See the code and results below:
>>
>> optm <-function(perm, verbose = FALSE)
>> {
>>    trace<-c()
>>    for (k in 1:perm){
>> trace[k]<-Rspatswap(rhox=0.6,rhoy=0.6,sigmasqG=0.081,SsqR=1)[1]
>>      perm[k]<-k
>>      mat<-cbind(trace, perm = seq(perm))
>>    }
>>    if (verbose){
>>      cat("***starting matrix\n")
>>      print(mat)
>>    }
>>    # iterate till done
>>    while(nrow(mat) > 1){
>>      high <- diff(mat[, 'trace']) > 0
>>      if (!any(high)) break  # done
>>      # find which one to delete
>>      delete <- which.max(high) + 1L
>>      mat <- mat[-delete, ]
>>      newmat<-apply(mat,2,mean)[1]
>>      sdm<-sd(mat[,1])
>>      sem<-sdm/sqrt(nrow(mat))
>>      maxv<-mat[1,1]
>>      minv<-mat[nrow(mat),1]
>>        }
>>    stats<-cbind(average=newmat,sd=sdm,se=sem,min=minv,max=maxv)
>>    stats
>> }
>>
>>> test<-optm(perm=20)
>>> test
>>          average           sd           se       min      max
>> trace 0.8880286 0.0009178193 0.0004589096 0.8870152 0.889241
>>
>>
>> itn<-function(it){
>> siml<-matrix(NA,ncol=5,nrow=length(it))
>>    for(g in 1:it){
>>     siml[g]<-optm(perm=20)
>>    }
>> siml<-cbind(siml=siml)
>> siml
>> }
>>
>>> ans<-itn(5)
>> Warning messages:
>> 1: In siml[g] <- optm(perm = 20) :
>>    number of items to replace is not a multiple of replacement length
>> 2: In siml[g] <- optm(perm = 20) :
>>    number of items to replace is not a multiple of replacement length
>> 3: In siml[g] <- optm(perm = 20) :
>>    number of items to replace is not a multiple of replacement length
>> 4: In siml[g] <- optm(perm = 20) :
>>    number of items to replace is not a multiple of replacement length
>> 5: In siml[g] <- optm(perm = 20) :
>>    number of items to replace is not a multiple of replacement length
>>> ans
>>            [,1]      [,2]      [,3]      [,4]      [,5]
>> [1,] 0.8874234 0.8861666 0.8880521 0.8870958 0.8876469
>>
>
> 1. Not reproducible code. Where does function Rspatswap come from?
>
> 2. You have several errors in function itn:
>      Argument it is a scalar: length(it) is 1. You need to do siml <- 
> matrix(NA,ncol=5,nrow=it)
> Next in the g-loop you want to fill row g so do: siml[g,] <- …..
>
> Finally why are you doing siml <- cbind(siml=siml)?
> Seems superfluous to me. Delete the line.
>
> Berend
>
>


        [[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.

Reply via email to