[R] Status of the "bs" Package

[krusty the klown]

Perhaps the problem stands on the rbs function which generates random samples
from the Birnbaum-Saunders distribution:
> library(bs)
> set.seed(1)
> x<-rbs(n=1000,alpha=0.5,beta=1.0)
> # sample mean
> mean(x)
[1] 1.117749
> # expected value
> beta*(1+alpha^2/2)
[1] 2.125
> # so different!

-- 
View this message in context: http://r.789695.n4.nabble.com/Status-of-the-bs-Package-tp827806p3033546.html
Sent from the R help mailing list archive at Nabble.com.