Search Archives

Hello, the multinom() procedure prints a lot of information during the fitting process. Is there a way to disable this "verbose" mode? Thanks, Bruno

Hello, I am looking for references about mixed models built on rank transformed data. Did anybody ever consider this topic? Thank you, Bruno ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Bruno L. Giordano, Ph.D. CIRMMT Schulich School of Music, McGill University 555 Sherbrooke Street West Montr...

Hello, Is there a way to compare the relative relevance of fixed and random effects in mixed models? I have in mind measures of effect size in ANOVAs, and would like to obtain similar information with mixed models. Are there...

Dear list, is anybody aware of R implementations of the overlapping clustering methods by Arabie and Carroll (ADCLUS, MAPCLUS and INDCLUS)? I found around only their Fortran implementation of the MAPCLUS model, which comes practically without any documentation. Alternatively, has...

Hello, with 4 different linear mixed models (continuous dependent) I find that my residuals do not follow the normality assumption (significant Shapiro-Wilk with values equal/higher than 0.976; sample sizes 750 or 1200). I find, instead, that my...

Hello, are there potential drawbacks in the use of rank (ordinal) variables as predictors in a "linear" mixed model for continuous dependents? I am trying this approach to create mixed effects models independent of the (monotonic) transform of the independent...

Dear R-list. I would need to know if the centroid of a cluster, whose coordinates are the average of the coordinates of the points in the cluster, is also the point that minimizes the sum of its distances from...

Dear list, I am trying to use lme to build the analogue of the following SAS MIXED random specification: random int+Variable1+Variable2 /subject = Subject group=Condition type=vc; which gives a Condition-blocked heterogeneity in the random effects variance...

Hello, I built 4 mixed models using different data sets and standardized variables as predictors. In all the models each of the fixed effects has an associated random effect (same predictor). What I find is that fixed effects with larger...

Hello, I looked for R packages which focused on order-restricted statistical inference, but I could find only the isoreg() function. I would need to test whether the means in my (repeated measures) data follow a given order, e.g...

Hello, I would greatly appreciate having some feedback on a problem. My data consist of two matrices of robust rank correlations between data collected in different conditions. The two matrices represent data from two between-subjects conditions. Each of the...

Dear list, In a mixed model I selected I find a non positive definite random effects variance-covariance matrix G, where some parameters are estimated close to zero, and related confidence intervals are incredibly large. Since simplification of the random...

Kappa() in library vcd kappa2() in library irr ckappa() in library psy cohen.kappa() in library concord...... Bruno ----- Original Message ----- From: "Philip He" <hydinghua at gmail.com> To: <R-help at stat.math.ethz.ch> Cc: "Rebecca Ding" <hua.ding...

I found this other more recent monograph: Silvapulle and Sen (October 2004), Constrained Statistical Inference: Order, Inequality, and Shape Constraints, Wiley-Interscience. http://ca.wiley.com/WileyCDA/WileyTitle/productCd-0471208272,descCd-reviews.html Thanks for the tip, Bruno ----- Original Message...

For ANOVA one option is the partial eta squared \eta^2_p: \eta^2_p=SSeffect/(SSeffect+SSerror) For multiple regression (continuous predictors) you might use the standardized parameter estimate, the regression coefficient you obtain standardizing the predictor: the larger...

Hello, I am writing down a general function to implement the bootstrapF method for repeated measures anova. I am passing the function several data frames: y = dependent subj = subject identifiers b = between-subjects factors (number = NB) w = within-subjects factors...

Hello, I am using logistic discriminant analysis to check whether a known classification Yobs can be predicted by few continuous variables X. What I do is to predict class probabilities with multinom() in nnet(), obtaining a predicted classification Ypred and...

----- Original Message ----- From: "Frank E Harrell Jr" <f.harrell at vanderbilt.edu> To: "Bruno L. Giordano" <bruno.giordano at music.mcgill.ca> Cc: <r-help at stat.math.ethz.ch> Sent: Friday, November 03, 2006 8:22 AM Subject: Re...

Well, just in case somebody is interested, the following R code gives the same estimates as the SAS code below: ######R code########### G.Data<-groupedData(Response~1|Subject,data=In.Data) G.Data$Condition<-as.ordered(G.Data$Condition) G...

Well, If posting a possible solution to one's own problem is not part of the netiquette of this list please correct me. Following Titus et al. (1984) one might use Cohen's kappa to have a chance-corrected measure...