similar to: automatic call generation for aov()

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 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?al, QC H3A 1E3 Canada http://www.music.mcgill.ca/~bruno/

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

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

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 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 ...

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

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

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'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 -

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 <-

Please see below Jon Flechsenhaar Boeing WNW Team Network Services (714)-762-1231 202-E7 -----Original Message----- From: Flechsenhaar, Jon J Sent: Wednesday, September 27, 2006 10:30 AM To: lartc@mailman.ds9a.nl Subject: FW: [LARTC] 2.6.14 - HTB/SFQ QoS broken? Please see below Jon Flechsenhaar Boeing WNW Team Network Services (714)-762-1231 202-E7 -----Original Message----- From:

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

Please see below Jon Flechsenhaar Boeing WNW Team Network Services (714)-762-1231 202-E7 -----Original Message----- From: Flechsenhaar, Jon J Sent: Wednesday, September 27, 2006 9:35 AM To: lartc@mailman.ds9a.nl Subject: [LARTC] 2.6.14 - HTB/SFQ QoS broken? Jody: I don''t know if you answer basic HTB questions but I''l try anyway. I''m implementing AF and EF with a

Hi all, I search the archive for finding a simple solution for using TukeyHSD with a multistratum aov result (a repeated emasure anova). The Question have been asked but I've found no clear answer. res<-aov(y~Fa*Fb+Error(Subject/(Fa*Fb)) ) I think that the problem is that res is an aovlist object instead of the "aov" object required by TukeyHSD. Is there an easy solution to

All: I have been told that HTB_HYSTERESIS might have some effect on rate calculations. This file is usually in /usr/src/linux/net/sched/sch_htb.c. If I change this file I have to re-compile. I am not quite sure how to do this. Can someone list the steps necessary to re-compile or point me to a doc that explains how to do so. Thanks. Jon Flechsenhaar Boeing WNW Team Network Services

Is it possible in a future version to be able to enable/disable mounts from inside the web interface the same way steamcast does now. I love icecast but i've been looking at steamcast and I like some parts of there web interface ie. bandwidth usage stats for each mount and in total for the server, and being able to enable/disable mounts. But I don't like the delay, it has that painful

Hi, I'd like to plot in in polar coordinates a line which is given as a vector of lengths and angles. library("plotrix")

I can use DSMARK to mark on the Egress side. Is there a way to mark/change the DSCP value of an incoming packet on the ingress side? Thanks. Jon Flechsenhaar Boeing WNW Team Network Services (714)-762-1231 202-E7