[R] Weibull simulation- number of items to replace is not a multiple of replacement length

[Phil Spector]

Perhaps the following will be instructive:

> mymat = matrix(1:6,3,2)
> p = 1
> mymat[p]
[1] 1
> mymat[p,]
[1] 1 4

When you index a matrix by a single subscript, it returns a
single element, corresponding to the columnwise vector representation
of the matrix.  I'm guessing you want to put the estimates into
the pth row (mymat[p,]) or pth column (mypat[,p]) of the matrix.

 					- Phil Spector
 					 Statistical Computing Facility
 					 Department of Statistics
 					 UC Berkeley
 					 spector at stat.berkeley.edu



On Thu, 16 Sep 2010, Halabi, Anan wrote:

> Hi,
> I write below code for simulation for weibull- estimating parameters by weibullMLE function,
> Although I define metrix for the variables still I got this message: number of items to replace is not a multiple of replacement length
> Any suggestion
>
>> est=matrix (NA, 2,2)
>> se=matrix (NA, 2,2)
>> for ( p in 1:2) {
> + sampleSize <- 20
> + shape.true <- 1.82
> + scale.true <- 987
> + sampWB <- urweibull(sampleSize, shape=shape.true, scale=scale.true, lb=0, ub=Inf)
> + sampWBmleWB <- weibullMLE(sampWB)
> + est[p]= sampWBmleWB$estimate
> + se [p]= sampWBmleWB$s
> + }
> Warning messages:
> 1: In est[p] = sampWBmleWB$estimate :
>  number of items to replace is not a multiple of replacement length
> 2: In se[p] = sampWBmleWB$s :
>  number of items to replace is not a multiple of replacement length
> 3: In est[p] = sampWBmleWB$estimate :
>  number of items to replace is not a multiple of replacement length
> 4: In se[p] = sampWBmleWB$s :
>  number of items to replace is not a multiple of replacement length
>> est
>         [,1] [,2]
> [1,] 1.507558   NA
> [2,] 1.525856   NA
>> se
>          [,1] [,2]
> [1,] 0.2513083   NA
> [2,] 0.2709073   NA
>
> ______________________________________________
> R-help at 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.
>