[R] qq plot for 3-parameter lognormal distribution

2009-12-21 Thread reezwan you
Hi: Could you please guide me how to do a qq plot for a 3-parameter lognormal curve. Here is an example from VGAM package to fit lognormal curve. I need to do a qq plot to check the fit. Thanks. lambda = 4 y = lambda + rlnorm(n - 1000, mean=1.5, sd=exp(-0.8)) fit = vglm(y ~ 1, lognormal3,

[R] Fitting Distribution to Data

2009-11-30 Thread reezwan you
Hello everyone: I tried to fit a Beta distribution on a right-skewed dataset using: fitdistr(temp,densfun=beta,start=list(shape1=3,shape2=2)) To assess the fit, I proceeded as follows: Using distribution parameters from the sample resulting from fitdistr() function, I generated 1000 samples

[R] Bayesian Residual Analysis and GOF for Ordinal Regression Models with Latent Variable Formulation

2009-01-20 Thread reezwan you
Hi: I am wondering if someone could please guide me to some package or implementation of Bayesian Residual Analysis and GOF for Ordinal Regression Models with Latent Variable Formulation. I've chekced arm package which uses glm formulation (howeve, I'm looking for Latent Variable/Data

[R] MCMC for sampling from ordinal logistic regression

2008-10-23 Thread reezwan you
Hi: Is there a R-function that can generate samples from the posterior distribution of an ordered logistic regression model (just like MCMCoprobit from MCMCpack in R). I will greatly appreciate some guidance in this regard. Thanks Reez __

Re: [R] MCMC for sampling from ordinal logistic regression

2008-10-23 Thread reezwan you
: Thursday, October 23, 2008, 11:36 AM reezwan you reez88 at yahoo.com writes: Is there a R-function that can generate samples from the posterior distribution of an ordered logistic regression model (just like MCMCoprobit from MCMCpack in R). I am confused. What's missing in the output

[R] MCMClogit: using weights

2008-10-18 Thread reezwan you
Hi everyone: I am just wondering how can I use weights with MCMClogit function (in MCMCpack package). For example, in case of glm function as given below, there is weights option in the arguments. Aparently there is no option of using weights in MCMClogit. glm(formula, family = gaussian,