R help - Sep 2010 - Weibull simulation- number of items to replace is not a multiple of replacement length

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

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