similar to: Using split.screen

I've got a colleague who's using a GLMM to analyse her data, and I've told her that she needs to include an offset. However, glmmPQL doesn't seem to allow one to be included. Is there anyway of doing this? Bob -- Bob O'Hara Rolf Nevanlinna Institute P.O. Box 4 (Yliopistonkatu 5) FIN-00014 University of Helsinki Finland Telephone: +358-9-191 23743 Mobile: +358 50 599

Moi! I have a problem where I want to find the ML estimate of a discrete parameter. I just want a function like optim that finds the max/min value for a function. Does anyone know of such a function for R? Thanks. Bob -- Bob O'Hara Rolf Nevanlinna Institute P.O. Box 4 (Yliopistonkatu 5) FIN-00014 University of Helsinki Finland Telephone: +358-9-191 23743 Mobile: +358 50 599 0540 Fax:

Moi! I was helping to teach a course on mixed models this week, and we came across a problem with coding more than one random effect in lme when they aren't nested. As an example, suppose we have an experiment where we sample moths from several populations, and place the moths on different trees, and measure a trait (in this case survival of offspring, but that's less important). We

I assume I'm missing something obvious here... The short form of my main question is: how do I extract variance components from an lme object? The longer form (plus optional supplementary question!): I'm looking at some quantitative genetics, and want to estimate two variance components so that I can then calculate a statistic called Qst from them. So I have this: reg1 <- lme(y ~

Hello from a presently sunny Helsinki! I've bee trying to repeat an analysis I did about 18 months ago (the reviewers of the paper want something adding to it). I'm using R1.5.0, but I couldn't see anything in the list of changes between this and 1.5.1 to suggest it would act any differently. The data is observational, on the changes in the population status of carabid beetles, and

Moi! I'm trying to add a rather long label to a y-axis, and it's so long that it won't fit into one line. However, things get strange when I try and split it over 2 lines. The problem seems to be the plus/minus symbol, which means I have to use expression(). I'm using R1.6.1 on Windoze 2000. As an example: thing1 <- expression(paste("log odds of survival probability (

I can only assume I'm betraying my ignorance here, but this is not what I would expect. I'm getting the following from a stepwise selection (with both step and stepAIC): > step(lm(sqrt(Grids)~ SE + Edge + NH), scope=~ (Edge + SE + NH)^2) Start: AIC= 593.56 sqrt(Grids) ~ SE + Edge + NH Df Sum of Sq RSS AIC <none> 2147.0 593.6 + Edge:NH 1

Currently I only have one machine backing up to another. I am backing up to a filesystem that is mounted as nosuid. The syncing process works perfectly except it fails on files that are suid. Example proc/ root/ tmp/ mknod tmp/mysql.sock : Invalid argument mknod tmp/orbit-daud/orb-12573780511489051058 : Invalid argument mknod tmp/orbit-daud/orb-139425434385535813 : Invalid argument mknod

Dear List, I have a couple of (~200) 3x3 transition probability matrices (ie each defines a Markov chain). They are all estimated from the same underlying process, so it ie meaningful to take their elemetwise mean and standard deviation. [1] First question: supposing that they are given in a list l, how do I get their elementwise mean and standard deviation? Fortunately, the mean of trans. prob.

Hola! Este mail llegó a r-help hoy, yo lo mandó también a esta lista. Es interesante con competiciones de predicción/clasificación! ¿Alguien que quiere cooperar? Kjetil ---------- Forwarded message ---------- From: Kjetil Halvorsen <kjetilbrinchmannhalvorsen@gmail.com> Date: Fri, Jan 29, 2010 at 12:19 Subject: Fwd: [R] Classification of supernovae - a challenge To: Kjetil Halvorsen

Full_Name: Roger Bivand Version: 1.7.1 OS: i686-pc-linux-gnu Submission from: (NULL) (129.177.30.18) identify() seg.faults when x is an empty list, when the n argument is given a positive value, avoiding the check for non-positive n (n <- length(x) when x is the non-existent x component of the empty list, and when neither x nor y are components of the list. suggested resolution: add test in

Dear list members, I've uploaded a new version of the Rcmdr package to CRAN. There are many additions and (I hope) improvements, as indicated in the relevant portion of the CHANGES file, reproduced below. As usual, comments, suggestions, and bug reports are appreciated. John ------------------------------------------- Version 0.8-4 o Minor bug fixes and additional input-error

A ZFS file system reports 1007GB beeing used (df -h / zfs list). When doing a ''du -sh'' on the filesystem root, I only get appr. 300GB which is the correct size. The file system became full during Christmas and I increased the quota from 1 to 1.5 to 2TB and then decreased to 1.5TB. No reservations. Files and processes that filled up the file systems have been removed/stopped.

The background: I'm trying to fit a Poisson-lognormal distrbutuion to some data. This is a way of modelling species abundances: N ~ Pois(lam) log(lam) ~ N(mu, sigma2) The number of individuals are Poisson distributed with an abundance drawn from a log-normal distrbution. To fit this to data, I need to integrate out lam. In principle, I can do it this way: PLN1 <- function(lam, Count,

Hola! I ask here since I learnt from this list that the LaTeX package listings should be good for typesetting R code. I encountered one problem: \begin{lstlisting} X %*% V \end{lstlisting} in the output the * in %*% disappears! same with %/%, etc, the / disappears. Any ideas? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. --

I encountered something strange with read.spss (package foreign, version 0.7 with R2.0.0 and version 0.8 with R2.0.1 beta, windows XP) I made a test file test.sav with SPSS version 11.5.1 containing only one numeric variable, with a value label for one value not occuring in the file. According to ?read.spss this should result in a factor, but it results in all NA. Using the argument

I get the following > library(foreign) Error in namespaceExport(ns, exports) : undefined exports: write.foreign Error in library(foreign) : package/namespace load failed for 'foreign' with rw2000dev as of (2004-09-17 Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

I just discovered that when the VGAM package (not on CRAN) is loaded, glm() doesn't work. This is because VGAM defines a family function() which gets found by glm() in place of the family function from stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I thought namespaces should protect us from this happening? Kjetil --

I am trying to install a local package and get this unexpected error: ---------- Making package UMSA ------------ adding build stamp to DESCRIPTION installing R files installing data files installing man source files installing indices Error: couldn't find function "na.omit" Execution halted na.omit of course is in package stats, and that is listed in the Depends field in

Hola! I had my hard disk and cpu fried, so had to reinstall debian from scratch, and now I have problems compiling R R-2.14.0. My system is wheezy with Xfce4 desktop. (-amd-64) ./configure make runs without problems, but make test dose not! Here is some of the output from make test: kjetil at kjetil:~/R/R-2.14.0$ make check make[1]: Entering directory `/home/kjetil/R/R-2.14.0/tests'