similar to: gamlss.dist

Dear All: I have succeeded in fitting a GAMLSS.dist model to growth data I am working with it. My aim is to create a matrix of predicted percentiles and the corresponding the fitted model's sigma mu nu by agebins. Q: How do it generate these parameters as in L M S per Cole and Green 1992? Here are my working codes. Name of fitted model is gamlssfit > Agebin<-seq(6,36,6)

Hi Folks! Just started using the gamlss package and I tried a simple code example (see below). Why the negative sigma? John > y <- rt(100, df=1)> m1<-fitDist(y, type="realline")Warning messages:1: In MLE(ll3, start = list(eta.mu = eta.mu, eta.sigma = eta.sigma, : possible convergence problem: optim gave code=1 false convergence (8)2: In MLE(ll4, start = list(eta.mu =

Dear R-helpers, How do I transform v <- c('insd-otsd', 'sppr-unsp') into c('insd--otsd', 'sppr--unsp') ? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903

Hi, I'm a bit puzzled by the gamlss fitting of exponential and lognormal data. Gamlss seems to think that exponentially distributed data fits better with a lognormal distribution, and vice versa. For example, X <- rexp(1000) Gexp <- gamlss(X~1,family=EXP) # X~1 is X tilde 1 GAMLSS-RS iteration 1: Global Deviance = 2037.825 GAMLSS-RS iteration 2: Global Deviance = 2037.825 Glno

How does one define PDFs as yet undefined in R, such as the ex- gaussian, the sum of two RVs, one exponential, one Gaussian? The PDF would then be the convolution of an exponential PDF, dexp(), and a normal, dnorm(). Kindly cc me in your reply to r-help. Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Dear R-listers, Can someone suggest a method for generating a finite Bernoulli sequence that is likely to have a given point-biserial correlation with an existing Bernoulli sequence? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall

Dear r-helpers, I downloaded http://popgen.unimaas.nl/~jlindsey/rcode/rmutil.tar (it was originally .tgz, but got unzipped by my browser). Can anyone give me detailed instructions on installing this and Lindsey's other packages on R version 2.4.0 (2006-10-03)---(powerpc- apple-darwin8.7.0, locale: C)? _____________________________ Professor Michael Kubovy University of Virginia

Dear R-helpers, How do I tell panel.polygon what greoup and panel it applies to whithin xyplot(y ~ x | c, groups = g panel = function(x, y, groups, ...){ panel.polygon(x = xpol[c, g], y = ypol[c, g], default.units = 'native') panel.xYplot(x, y, groups, ...) llines(x = c(1, 6), y = c(-24.283333, 35.941667), lwd = 2, lty = 3, col = 4) } x[c, g] and y[c, g] describe

Dear R-helpers, I would like the bands in an xYplot() to be gray without dark edges. I've tried playing around with col.bands, but that seems not to affect the edges, although changin col.fill does change the fill as I would expect. I would appreciate pointers. _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box

Dear r-helpers, After updating to R version 2.2.1, 2005-12-20, powerpc-apple- darwin7.9.0 (from binary), JGR wouldn't launch. I then deleted everything I could find related to JGR, and reinstalled it from the latest binary. Both R and JGR are installed in /Applications/ \~LocalApps/mathStat/. Advice on how to identify the source of the problem? My machine: Machine Name:

Dear r-helpers, When use latex() on a matrix, I set the option dcolumn = T. As a result, in front of each tabular I get \newcolumntype{.}{D{.}{.}{-1}} The LaTeX compiler complains about these multiple redeclarations. _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room

Dear Friends, Please note that in the following CI lower > CI higher: > require(lmer) > require(gmodels) > fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy) > ci(fm2) Estimate CI lower CI upper Std. Error p-value (Intercept) 251.66693 266.06895 238.630280 7.056447 0 Days 10.52773 13.63372 7.389946 1.646900

Dear R-helpers, This produces a red title. How do I make it black, without making the labels black? legend(title = 'Channel Integration', 0.7, 0.3, c('loudness','luminance'), col = c('red', 'blue'), text.col = c('red', 'blue'), lty = 1, xjust = 0.5, bty = 'n') _____________________________ Professor Michael Kubovy University

Dear R-helpers, I need to generate all the binary sequences of length n (here n = 8) that start with 1 and have no fewer than two of each digit, and are not cyclic permutations of each other. Here is what I have done: len <- 8 df <- as.data.frame(numeric(2^(len - 1)) %o% numeric(len)) require(partitions) for (i in 1:2^(len - 1)) df[i, ] <- binary(i, dim = len)[[1]] df <-

Here's a small dataset: type <- c('hierarchical','partial','single','complete','single +hierarchical','single+partial','partial+hierarchical','single+partial +hierarchical') freq <- c(1455,729,688,65,29,28,16,17) lodds <- log(freq/(3027 - freq)) dotplot(type~lodds) I would like to have the x-axis have ticks at nice,

Dear r-helpers, After successfully running require(nlme) vfr.lmL <- lmList( estimate ~ (slant + respType + visField + hand)^2 | subject, vfr ) pairs(vfr.lmL, id = 0.01, adj = -0.5) # Pinheiro & Bates (p. 141) produces the following error: Error in sprintf(gettext(fmt, domain = domain), ...) : object "form" not found Any guesses as to what I may have done wrong?

df1 <- data.frame(a = LETTERS[1:2], b = LETTERS[3:4], c = 1:2) I am looking for an idiom that swaps the elements of df$a and df$b when (e.g.) df$c == 2, resulting in df2 <- data.frame(a = LETTERS[c(1, 4)], b = LETTERS[c(3, 2)], c = 1:2) _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA

Hi, I'm playing around with gamlss and don't entirely understand the sigma result from an attempted lognormal fit. In the example below, I've created lognormal data with mu=10 and sigma=2. When I try a gamlss fit, I get an estimated mu=9.947 and sigma=0.69 The mu estimate seems in the ballpark, but sigma is very low. I get similar results on repeated trials and with Normal and

Hi everyone, I''m trying to use function vcov to extract the covariance matrix from a GAMLSS object. But I''m getting some strange errors and I was hoping someone could help me out? Vcov works with the same model for lm and glm objects, but not gamlss objects. I''ve searched various help sites to no avail. Its very possible the reason is that vcov failed though,

Hi, I am running a negative binomial model using Gamlss and when I try to include random effect, I get the following message: Warning messages: 1: In vcov.gamlss(object, "all") : addive terms exists in the mu formula standard errors for the linear terms maybe are not appropriate 2: In vcov.gamlss(object, "all") : addive terms exists in the sigma formula standard