Re: [R] SE from nleqslv

2012-03-20 Thread FU-WEN LIANG
On Tue, Mar 20, 2012 at 4:24 AM, Berend Hasselman wrote: > > > On 20-03-2012, at 01:01, FU-WEN LIANG wrote: > > > Dear R-users, > > > > I use the "nleqslv" function to get parameter estimates by solving a > > system > > of non-linear e

[R] SE from nleqslv

2012-03-19 Thread FU-WEN LIANG
Dear R-users, I use the "nleqslv" function to get parameter estimates by solving a system of non-linear equations. But I also need standard error for each of estimates. I checked the nleqslv manual but it didn't mention about SE. Is there any way to get the SE for each estimate? Thank you very mu

Re: [R] Scale parameter in Weibull distribution

2012-03-06 Thread FU-WEN LIANG
hazard for T is 0.004. I'm wondering if it's possible to get the estimated baseline hazard by the output from survreg. Thank you very much. On Tue, Mar 6, 2012 at 5:06 PM, David Winsemius wrote: > > On Mar 6, 2012, at 5:53 PM, FU-WEN LIANG wrote: > > Hi all, >> >&

[R] Scale parameter in Weibull distribution

2012-03-06 Thread FU-WEN LIANG
Hi all, I'm trying to generate a Weibull distribution including four covariates in the model. Here is the code I used: T = rweibull(200, shape=1.3, scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3)) C = rweibull(n, shape=1.5, scale=0.008) #censoring time time = pmin(T,C) #observed time is m

[R] Generate a Weibull regression data

2012-02-22 Thread FU-WEN LIANG
Hi all, I'm trying to generate a Weibull distribution including four covariates in the model. Here is the code I used: # Generate survival time T = rweibull(200, shape=1.3, scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3)) C = rweibull(n, shape=1.5, scale=0.008) #censoring time time = pmin

Re: [R] Constraint on one of parameters.

2012-02-11 Thread FU-WEN LIANG
> HTH > Rubén > > -----Mensaje original- > De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] En > nombre de FU-WEN LIANG > Enviado el: jueves, 09 de febrero de 2012 23:56 > Para: r-help@r-project.org > Asunto: [R] Constraint on one of parameters. &

[R] Constraint on one of parameters.

2012-02-09 Thread FU-WEN LIANG
Dear all, I have a function to optimize for a set of parameters and want to set a constraint on only one parameter. Here is my function. What I want to do is estimate the parameters of a bivariate normal distribution where the correlation has to be between -1 and 1. Would you please advise how to