similar to: dfs in lme

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:

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

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

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

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

> I tried to copy a 8GB Xen domU disk image from a zvol device > to an image file on an ufs filesystem, and was surprised that > reading from the zvol character device doesn''t detect "EOF". > > On snv_66 (sparc) and snv_73 (x86) I can reproduce it, like this: > > # zfs create -V 1440k tank/floppy-img > > # dd if=/dev/zvol/dsk/tank/floppy-img

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 would like to fit an lm in which a subset of the explanatory variables are linearly dependent. Thus I would like to include the restriction that all betas of these variables sum to 1. Is there a way to this in R? Or is this what happens automatically if I set singular.ok= T (which is the default, I believe)? Thanks for your help. Lorenz --

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

server: ms 2003 with ads client: debian 3.1/samba 3.0.14 smb.conf: .. [global] workgroup = SP-GRUPPE password server = 10.85.117.150 realm = SP-GRUPPE.DE encrypt passwords = no server string = %h server (Samba %v) obey pam restrictions = yes passdb backend = tdbsam, guest passwd program = /usr/bin/passwd %u passwd chat = *Enter\snew\sUNIX\spassword:* %n\n *Retype\snew\sUNIX\spassword:* %n\n .

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, 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 R-Users, I have a newbie-question about producing several plots with five variable pairs within in one code (loop). Problem: I define two objects, erveryone has five vectors (columns /variables): dimensionen <- data.frame(meinspss$attr_diff_gesamt, meinspss$finanz_diff_gesamt, meinspss$leist_diff_gesamt, meinspss$soz_diff_gesamt, meinspss$wert_diff_gesamt) gruppe <-

Hi, I've got a browsing problem, I'm searching the error for a long time now so I think betimes it's a bug ;-) For debugging purposes, I've left away everything in the network that's not necessary. There's now only one XP-Prof workstation (domain member) and one samba server (3.0.22, acting as PDC) anymore. Bind9 is the local DNS server; Forward and PTR entries are