similar to: lm with rank-deficient X matrix

Dear all, I think this question has been raised very recently but I think I didn't see repsonses on the list nor did I find the thread in the archives. (If it is there I'd appreciate directions ...) I am about to install Linux on a Compaq notebook. One important thing I want to do is being able to run R. I have fiddled with different distributions of Linux (LST, SuSE, Debian) and can

Dear all, well, today I just seem to run into these problems ... (I am still working on R 0.63.1 on a SuSE 5.2). I would like to have a data.frame that also contains character (and NOT factor) vectors. Two things that I tried (the second worked in R 0.62.2): > test.df <- data.frame (X= character (4)) > test.df$X [1] <- "a" Warning: invalid factor level, NAs generated or:

Dear all, Many thanks to all who always so quickly respond with helpful hints! Currently I am working with R 0.63.1 on SuSE 5.2 which I installed this morning within about half an hour: my deepest compliments to the developpers! (I had still 0.62.2 before.) This time I am trying the following: > xy <- list (a= c (1, 2, NA), b= c (2, 3)) > list.name <- "xy" > get

Dear all, I am trying to use ace on a 446 x 14 Matrix, with two categorical variables in the first two columns. If I try to run it I get: > aceiH <- ace (data.mat, resp, cat= 1:length(x.cat)) Error: NAs in foreign function call (arg 3) There is only one .Fortran call in ace () with the third argument being x = t(x). Thus I tried to debug ace, i.e. I went through ace until just before

Dear all, The following works with a simple result of lm but gives the error with a more complicated one (could this be due to singularities that occurr in the model?). > cooks.distance (X) Error: NAs in foreign function call (arg 6) > class (X) [1] "lm" But, actually, I would like to ask a more general question: Is there a place (I searched in the CRAN site) where I find

Dear all, is this a bug or do a miss something (this concerns version 0.63.1 on Linux)? > boxplot (c (1, 2, 3, 4, 5), ylim= c (-5, 5)) Warning: parameter "ylim" couldn't be set in high-level plot() function > boxplot (c (1, 2, 3, 4, 5), pars= list (ylim= c (-5, 5))) Warning: parameter "ylim" couldn't be set in high-level plot() function > and actually in

Dear all, in a current project I have a pretty huge .Rdata. Thus I was working with R --vsize 100 --nsize 1000000. Today when I tried to restart R I get the following error message: Error: a read error occured Fatal error: unable to restore saved data (remove .RData or increase memory) I increased memory up to --vsize 180 and --nsize 2000000, but the error reoccurrs. Is there a way to know

Dear all, I am working with R Version 0.62.2 (July 10, 1998) on a SuSE Linux 5.2. I have troubles in understanding the function image (). Can someone enlighten me? I'd like to draw a square image-plot (with every little part of the image a square as well) and add some (straight) lines to it. Also, I create an empty z matrix which I need to fill step by step. Thus I tried to find out how

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Hello I just installed openssh-3.7.1p1 on IRIX64 6.5 04100802 IP27. The compilation failed since the header file openbsd-compat/inet_ntoa.h was missing. I downloaded the package from different sources including ftp.openbsd.org, there was no difference. Since inet_ntoa.c has not been changed, copying the inet_ntoa.h from openssh-3.6p1 fixed the problem. The package should be fixed anyway.

I have just realised that I sent this to Per only. For those interested on the list: -----Original Message----- From: Gygax Lorenz FAT Sent: Tuesday, September 14, 2004 4:35 PM To: 'Per Tor??ng' Subject: RE: [R] glmmPQL and random factors Hi Per, > glmmPQL(Fruit.set~Treat1*Treat2+offset(log10(No.flowers)), > random=~1|Plot, family=poisson, data=...) > > Plot is supposed

Dear R users, I have not found anything on this in the archives. Does anyone know whehther the parameter use= is not functioning in cor or enlighten me what it is supposed to do? My R version is "R version 1.8.1, 2003-11-21" on Windows 2000. I am hoping to be able to update to 1.9.1 as soon as it has appeared (we are not allowed here to install software on our own and thus I am trying

Dear R users, Is there something like predict (..., type= 'response') for glmmPQL objects or how would I get fitted values on the scale of the response variable for the binomial and the poisson family? Any pointers are appreciated. Thanks, Lorenz - Lorenz Gygax, Dr. sc. nat. Tel: +41 (0)52 368 33 84 / lorenz.gygax at fat.admin.ch Tag der offenen T??r, 11./12. Juni 2004:

Dear R-mixed-effects-modelers, I could not answer this questions with the book by Pinheiro & Bates and did not find anything appropriate in the archives, either ... We are preparing a short lecture on degrees of freedom and would like to show lme's as an example as we often need to work with these. I have a problem in understanding how many dfs are needed if random terms are used for

Dear all, I am doing some data handling and tabualtions and then I am using glmmPQL to fit a binomial hierarchical model (which I assign to a new object). If I do a save.image () or a q ('yes') after this, I get the following warning messages: Warning messages: 1: namespaces may not be available when loading 2: names in persistent strings are currently ignored Sometimes R seems to

We have a samba share which is a NFS mounted filesystem with user quota. The samba server runs SuSE 9.2, the samba version is 3.0.20b. The filesystem is nfs mounted from a NetApp fileserver. If a user on a windows (xp or 2k) client copies some files onto his samba share and needs more space than the quota allows, the copy operation ends without any error message! As I have seen in the samba

Hi all, I am not a very proficient R-user yet, so I hope I am not wasting people?s time. I want to run a linear mixed model with 3 random factors (A, B, C) where A and B are partially crossed and C is nested within B. I understand that this is not easily possible using lme but it might be using lmer. I encountered two problems when trying: Firstly, I can enter two random factors in lmer but

I don't think, this has been answered: > I'm trying to run a 3-way within-subject anova in lme with 3 > fixed factors (Trust, Sex, and Freq), but get stuck with handling > the random effects. As I want to include all the possible random > effects in the model, it would be something more or less > equivalent to using aov > > > fit.aov <- aov(Beta ~ >

Dear all, I do not seem to grasp how contrasts are set for ordered factors. Perhaps someone can elighten me? When I work with ordered factors, I would often like to be able to reduce the used polynomial to a simpler one (where possible). Thus, I would like to explicetly code the polynomial but ideally, the intial model (thus, the full polynomial) would be identical to one with an ordered factor.