Re: [R] Generating by inverting function

2012-09-30 Thread Salma Wafi
Dear All, Thank you very much for your kind and response.   Actually, S(t) = exp(b* F(t)) represent one of survival models, where F(t) is known cumulative distribtion function. So, that is why I need to generate the survival times T based on this model. Then for generating time I have to invert

[R] Generating by inverting function

2012-09-30 Thread Salma Wafi
Hello, I am trying to generate random survival times by inverting the function,  S(t)= exp(b*F(t)), where b is constant and F(t) is some cumulative distribution function, let say that F(t) is cdf of normal distribution or any others distributions.   as we know that S(t) has uniform distribution

Re: [R] fitting lognormal censored data

2012-09-03 Thread Salma Wafi
Dear Berend Thanks for your response. I will write my code and estimate the parameters using functions optim() or nlminb(). For the code that I am trying to build for estimation the parameters using mle via EM ahgorithm, I will rewrite it and reformulated my question about, but before that I nee

[R] fitting lognormal censored data

2012-08-30 Thread Salma Wafi
Hi , I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likel

[R] Estimation parameters of lognormal censored data

2012-08-28 Thread Salma Wafi
Hi, I am trying to get the maximum likelihood estimator for lognormal distribution with censored data;when we have left, interval and right censord. I built my code in R, by writing the deriving of log likelihood function and using newton raphson method but my estimators were too high " overest