Ok here's the lines of codes:

sIr=replicate(50000,ts(arima.sim(list(ar=0),n=150,innov=rgev(150,xi=0.2,
mu=inP[1,5],sigma=inP[1,4]),start.innov=rgev(150,xi=0.2,
mu=inP[1,5],sigma=inP[1,4]))))

sIrr=data.frame(cbind(sIr))

EsIpr=sapply(sIrr[1:100,1:20000],function(m) pargevn(lmom.ub(m)))

kIpr=EsIpr[3,]#try to call for kappa
Error in EsIpr[3, ] : incorrect number of dimensions

#it turn out that the Eslpr is a list. It is kinda weird because previously
I run the same routine for

Eslr=sapply(sIrr[1:150,20000]), function (m) pargevn(lmom.ub(m))# it works
out fine when I call
kI=EsIr[3,]

Thanks!


On Tue, Oct 16, 2012 at 12:35 PM, Muhammad Rahiz
<muhammad.rahiz@gmail.com>wrote:

> Maybe something along the following lines?
>
> m <- rep(list(rnorm(10)),3) # dummy data
>
> data <- c()
> for (a in 1:length(m)) {
> data[a] <- m[[a]][1] }
>
> Muhammad
>
>
> On Tue, Oct 16, 2012 at 1:17 PM, Balqis <aehan3616@gmail.com> wrote:
>
>> Hi all,
>>
>> I have a list of 20000 data, and the list look like below. I wonder what
>> is
>> the simplest way to extract 'kappa' value (or 'xi' or 'alpha' for the
>> matter)  from each of the data. How can I simply code it without having to
>> change the list to a dataframe first? Many thanks!
>>
>>
>> $X19997
>>          xi       alpha       kappa
>> 784.7718640 165.4065141  -0.2709599
>>
>> $X19998
>>          xi       alpha       kappa
>> 860.7832604 215.1144980  -0.1769741
>>
>> $X19999
>>          xi       alpha       kappa
>> 839.7692675 196.0103874  -0.3848809
>>
>> $X20000
>>          xi       alpha       kappa
>> 847.6117040 212.3279185  -0.2576029
>>
>> regards,
>> aehan
>>
>>         [[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.
>>
>
>
>
> --
> Muhammad Rahiz
> (za) +27 071 719 0104
> (skype) muhammad.rahiz
>
>

	[[alternative HTML version deleted]]


