similar to: lmer applied to a wellknown (?) example: summary

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being used elsewhere in the litterature. Because this example is so thoroughly described, I thought

I am trying to code a steepest ascent algorithm to optimize parameters used in a survivor function type problem. My unknown parameters (alpha, Beta0, and Beta1) for which I have been able to optimize using Newton's method. I keep getting an error because my alpha becomes negative and I can't calculate the likelihood. Here is my log likelihood I am optimizing (in LaTex): l=\sum _{ i=1 }^{

Dear R users, I successfully installed SciViews the other day. However, when I try to run it now, the command/script window does not appear. Strange. Well, actually I see a tendency to a script window (in the lower part of the sciview window where it is suppose to be) during start up, but when the program is entirely open, the script window is gone. I have tried to uninstall and reinstall

Dear all, I have some csv-files (originating from Excel-files) containing empty cells. In my example file I have four variables of different classes, each with some empty cells in the original csv-file: > test <- read.csv2("test.csv", dec=".") > test id id2 x y 1 a 1 NA 2 b e NA 2.2 3 f 3 3.3 4 c g 4 4.4 > class(test$id) [1]

------- Forwarded message follows ------- From: Petr Pikal <petr.pikal at precheza.cz> To: Henrik Parn <henrik.parn at bio.ntnu.no> Subject: Re: [R] paired t-test. Need to rearrange data? Date sent: Tue, 08 Aug 2006 16:13:47 +0200 Hi Uff, it takes me a bit headache but this shall do it ?unstack ?table ?t new.list<-unstack(test.data,y~id)

Dear all, I run a linear model with three significant explanatory variabels x1: a factor with 4 levels x2 and x3: factors with two levels each x4: continuous model <- lm(y ~ x1 + x2 * x3 + x4) <> The data is not perfectly balanced between the different factor-combinations and I use treatment contrasts. <> With an aov.object, I assume I could have used model.tables(aov.object,

Dear all, I am trying to plot some data with differing range in y-values with type="b", adding error bars and break the y-axis into two parts, one lower part from 12 to 20, and one upper part from 34 to 40. I have tried to follow the basic ideas from the script provided here by Jim Lemon: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/56487.html My attempt looks like this:

Hi, I'm experiencing weird behaviour in Xapian/Omega. Here's my situation: * I'm indexing data through scriptindex, the output is as follows: ***************** Replace: 6706476 Replace: 6706477 Replace: 6706478 ***************** Since there is no output after entering several enters, I assume scriptindex is indexing data. * The search I perform on omega is as follows:

> -----Original Message----- > Date: Tue, 26 Apr 2011 13:35:20 +0100 > From: James Aylett <james-xapian at tartarus.org> > Subject: Re: [Xapian-discuss] Omindex: what are the default numbered > indexes? > To: <xapian at catcons.co.uk> <xapian at catcons.co.uk> > Cc: 'Xapian Discussion' <xapian-discuss at lists.xapian.org> > Message-ID:

Hello :-) How can the default omega URL be prettified? http://<host_ID>/cgi-bin/omega is working fine, giving us all omega's default CGI parameters. Now we want multiple databases which could be accessed using http://<host_ID>/cgi-bin/omega?DB=<index_ID> but this is starting to get messy. It will get messier when we start to customise templates with

Dear all, I have a question on handling of missing values in lmList. My data set have continuous predictor and response, x and y, and a grouping variable group.id. All these variables have NAs and the data set also has several other variables that also contains NAs. To create the lmList-object seems to work fine: y.list <- lmList(y ~ x | group.id, data=mydata, na.action=na.omit) However,

HI there, I am running Fedora 15. I run everything as SU by default. I am trying to get Xapian/Omega installed and running. I installed verion 1.26 of xapian from here: http://xapian.org/download I followed the example listed here: http://trac.xapian.org/wiki/OmegaExample I searched the archives and found a reference with some questions from other users so I'll answer

Hi guys, we are running a fileserver for a company with ~10 Mac OS Mavericks (not all up to date). The Samba server is running on Centos 5 with official Samba from Centos. We see big big lagging while browsing the fileserver. Some directories du to the activy of the company contains more than 4000 files (small one). Everything is fine if we connect with Windows or Linux machine, but with Mac Os

I'm using maildrop to filter mail like: xfilter "/usr/bin/spamc" xfilter "/usr/bin/bogofilter -ep" I want to change to deliver (so I can use sieve) but can deliver do this?

Dear R core team Using the following code produces an empty plot (similar to col = NA): > plot(1:9, col = factor(rep(1:3,3), labels = c("red", "blue", "black"))) My question: Shouldn't one get at least a warning (or an error) if one tries to use a factor as col argument? Thanks for an answer. Regards, Christoph Buser

Dear R gurus I've a question concerning the transformation of a new S4 class (that I defined) and the existing class "POSIXlt". I did the following: ## Definition of the new class: setClass("dtime", representation(data = "POSIXlt"), prototype(data = as.POSIXlt("2004/06/01"))) ## Transformation between the new class "dtime"

Dear R core team Please apologize for posting the same question twice on R-help and R-devel. Since I was not sure which list is appropriate I tried R-help (Tue Mar 22), but got no answer. Now I do not know if the formulation of my question was unclear or the question is not so easy to answer or to easy (what I do not hope). My problem: I create a new S4 class, containing one slot, data (of

Hy all, I apologize i've used rownames instead of colnames in my first exemple. (that's why i changed the mail object) I was on mars when i wrote my question... Every answer where correct ... but, because i've made a wrong question, people wern't able to understand... So finally i'm speaking about colnames : Let's be more precise: I've got a query that gives me for

Dear R-core There is a different output for the apply function due to the change of unlist as mentioned in the R news. Newly, applying as.factor() (or factor()) in str(dat <- data.frame(x = 1:10, f1 = gl(2,5,labels = c("A", "B")))) (d1 <- apply(dat,2,as.factor)) newly returns a character matrix while in R-2.3.1 the same command resulted in an integer matrix that was

Dear R core team I've a proprosal to improve the function interaction.plot. It should be allowed to use type = "b". This can be done by changing the function's header from function( , type = c("l", "p"), ) to function( , type = c("l", "p", "b"), ) Then it works. This type = "b" is useful, if