R help - Dec 2015 - The Birnbaum-Saunders Distribution simmulation Problem

Dear i optimized the gama and beta value using simmulation. if i run this code
it generate very small value of beta. Could any one help me in this regard. i
use gbs package to generate data.


gama=1.0
beta=1.3
n=25
iterCount=1000
for(i in 1:iterCount){
  x<-rgbs(n,gama,beta)
  P<-function(theta,x){
    n<-length(x)
    gama<-theta[1]
    beta<-theta[2]
    -n*log(5.013)+ n*(theta[1])^-2-n*log(theta[1])-n/2*log(theta[2])+
sum(log(theta[2]+x)-0.5/theta[1]^2*(x/theta[2]+theta[2]/x))}}
  P.out<-optim(theta<-c(gama,beta),ll.wd,x=x,method =
"Nelder-Mead",hessian=FALSE)
  gamhat<-P.out$par[1]
  betahat<-P.out$par[2]


	[[alternative HTML version deleted]]