similar to: Rank transformation and the linear mixed model

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 (absolute) standardized parameter estimates have also a higher estimate of the related random effect. In other words, the higher the average of the absolute

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 then compute the percentage P(obs) of objects classified the same in Yobs and Ypred. My problem now is to figure out whether P(obs) is

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. A<B=C<D. I took a look at the monograph by Barlow et al. (1972) on this topic and found that for my case the null hypothesis is always A=B=C=D. This might be

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-covariance matrix. Needless to say, I have a hard time in specifying Condition-specific heterogeneities in the variance-covariance

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 (number = NW) after grouping of all these variables in a single data frame the aov() call looks like this:

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 anybody used this Fortran program, and is willing to share some knowledge about its use? Thank you, Bruno

Dear R users, I m quite a novice in using R for factor analysis and I would need some help to choose the right function. I have a contingency table and I would like to perform a Correspondence analysis on this table, followed by a hirarchical clustering of my variables projected in on the first principal components. Here are my question : - what is the more appropriate function to do so ...

Hi, In view of the problems I was having with DomU network timeout after a live migration (I posted that problems here a while ago but never got anything except from private emails) I finally updated my Debian/Squeeze Dom0s last night to a new kernel, from 2.6.32-23 to 2.6.32-26. Now live migration just hangs...Any ideas? Xen-related Debian packages (all from Debian repositary, except drbd

It would be great to know if and where an R code for Common Principal Component Analysis is available. Thanks, Daniel Daniel Berner Redpath Museum & Dept. of Biology McGill University 859 Sherbrooke St. W. Montreal, QC, H3A 2K6 Canada Phone: 514-398-4086 ext. 00908 Fax: 514-398-3185 Email: daniel.berner at mail.mcgill.ca

I am playing with the a 1-way anova with and without the "-1" option. I have a simple cooked up example below but it behaves the same on a more complex real example. From what I can tell: 1) the estimated means of the different levels are correctly estimated either way (although reported as means with the -1 and as contrasts without the -1 as expected) 2) the residuals are

Working on Windows I have had to deal with CSV files that, unfortunately, contain embedded Control-Zs, i.e. ASCII character 26 in decimal, and the readLines() function in R on Windows (2.15.2 and 3.0.0) appears to truncate at the control-Z. There is no problem at all on Ubuntu Linux with R 3.0.0. Am I mistaken or is this genuine? # Create a small file with embedded Control-Z h3 <-

Hi all, What is (or how to I change) the maximum width of a pdf image. I am trying to make a pdf image: pdf(file="out.pdf",width=200,height=20) and I can't get the image to be more that 200 inches wide i.e. I get an empty out.pdf . Is there a solution around this. Thanks. Vince -- +-----------------------------------------------------------+ | Vincenzo Forgetta

Hello ... I tried to generate .fig figures with R, using the xfig() function ...When I open the figures using xfig software under linux ... the foreground color change strangely ... whereas when I display the same figure in R or when I saved it in using postscript() function there is no problem. Any idea about the behavior of the colors when using xfig ?? Thanks in advance Christophe Grova

I've downloaded R for windows (9.0.1) and it is great! I've converted all my lecture notes for my GLM course to run on R (they are available on my web page below). I must admit I particularly like the default contrast options, which are identical to GLIM. Also I like the gl function - very useful! I have a couple of questions/bugs: 1. predict.glm doesn't work, but predict.lm does -

Hello, I'm stumped. I'm trying to join a machine to a domain and I get this : # samba-tool domain join campus.mcgill.ca DC -Udavid.salib at CAMPUS.MCGILL.CA Finding a writeable DC for domain 'campus.mcgill.ca' Found DC DC01.campus.MCGILL.CA Password for [david.salib at CAMPUS.MCGILL.CA]: workgroup is CAMPUS realm is campus.MCGILL.CA checking sAMAccountName Adding

Hi everyone, I have built an R package and for the sake of speed I have decided to rewrite some part of the code in C++. In my original R code I use the pnorm, qnorm, rnorm, pgamma, dgamma, rgamma, rbeta and runif function. First I was thinking in going with the boost libraries, but I noticed the functions described above are available within the R.h header file (or is it Rmath.h?). So my

Hello everyone! First of all, please note that I'm working under Windows 7. I have written a Gibbs sampler in R and I'm now in the process of translating it in C++ to increase the speed. I'm relatively new to C++ and this question may be trivial, but my search so far have been unsuccessful. In my Gibbs sampler I am using some basic R functions (like the rep function for example)

Hello, I need to perform a mixed-model (with nesting) MANCOVA, using Type III sums of squares. I know how to perform each of these types of tests individually, but I am not sure if performing a mixed-model MANCOVA is possible. Please let me know. Erika <>< <>< <>< <>< <>< <>< <>< Erika Crispo, PhD candidate

Dear all, I'm trying to create my own package by following for instance http://epub.ub.uni-muenchen.de/6175/2/tr036.pdf In step 6 (page 19) it says that the command 'R CMD build <mypackage>' should result in a tar.gz file. However, after running the command I end up with a tar file only (no gz extension). I also get the following warnings in the command prompt C:\Program

Hi, I'm having some problems installing packages from the source files on Windows, using R CMD INSTALL pkg. I'm running WindowsXP, and I've followed the instructions as per the README.packages file from the R installation. I ran into a hitch, with the install failing following a "hhc: not found" warning. I figured out that this was related to the html help workshop.