Displaying 1 result from an estimated 1 matches for "estiid50".
2012 Jun 27
2
how to apply the same function to multiple data set
...rs,
I'm trying to repeat the same procedure to 1000 data set. I know this is
very easy, but I got stuck finding the right and fastest way in running it.
IID50=Riidf[1:50,1:1000] #where IID50 is a dataframe consist of 1000 time
series(as column) and 50 time scales (row).
#what I tried to do:
estIID50=rep(NA,1000)
for (i in 1:1000)
estIID50[i]=pargev(lmom.ub(IID50[1:50,i]))
#warning message
In estIID50[i] = pargev(lmom.ub(IID50[1:50, i])) :
number of items to replace is not a multiple of replacement length
#pargev is a function from lmomco package. I would like to apply it to the
1000 set o...