similar to: Preferable contrasts?

Dear All, Comparing linear mixed effect models in SAS and R, I found the following discrepancy: SAS R random statement random subj(program); random = ~ 1 | Subj -2*loglik 1420.8 1439.363 random effects variance(Intercept) 9.6033 9.604662

Summary: proc mixed vs. lme The objective of this summary is to help people to get more familiar with the specification of random effects with proc mixed or lme. Very useful are the examples of Ramon Littell's book: "SAS System for Mixed Models (1996)" (http://ftp.sas.com/samples/A55235) The same data set's are kindly made available by Douglas Bates in the

Dear Experts, May be the problem is still solved, however I tried to find the answer in the archives: I use: > R.version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 7.1 year 2003 month 06 day 16

Hi, I'm using: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.1 year 2002 month 06 day 17 language R and I would like to apply: > update.packages() trying URL

Dear R-helpers, First of all I have nothing against the French language! But now my problem, yesterday I installed R 2.1.1 and I had to experience that my RGui is speaking French. My windows locals is French (Switzerland). I'm used to English and I want to reset my RGui to English. I was seeking for the solution in the archives, however not successfully. By the way the searchable archives

Dear all,I am a new user of WinBUGS and need your help. After running the following code, I got parameters of beta0 through beta4 (stats, density), but I don't know how to get the prediction of the last value of h, the variable I set to NA and want to model it using the following code.Does anyone can given me a hint? Any advice would be greatly appreciated.Best

Dear all, I am reading with R some measurements vectors from an external device. The length of each vector (containing numbers) is not constant. I would like by doubling the numbers of elements I select to shrink the vector to a new one Let's call the vector I read from the device as oldvector and the shrinked version of it as newVersion I want to have newVersion[1]<-

Found this code: ftp://ftp.samba.org/pub/unpacked/ntdb/lib/ccan/ilog/ilog.c Tests show that it's faster to use the following code in FLAC__bitmath_ilog2_wide(): static const unsigned char DEBRUIJN_IDX32[32]={ 0, 1,28, 2,29,14,24, 3,30,22,20,15,25,17, 4, 8, 31,27,13,23,21,19,16, 7,26,12,18, 6,11, 5,10, 9 }; FLAC__uint32 v; int m;

Dear all, I have a question about linear models. Is there something comparable to the estimate statement in SAS? E.g. proc mixed data=test; class x1 x2; model y=x1|x2 /s; estimate 'x1' x1 -1 1 x1*x2 -0.5 -0.5 0.5 0.5; estimate 'x2' x2 -1 1 x1*x2 -0.5 0.5 -0.5 0.5; run; Thanks for the help, Dominik Dominik Grathwohl Biostatistician Nestl? Research

Dear All: Could you please tell me how I can draw figure formatted like figure 4.18 of MASS (4th) with the attached data set? Thanks Zhongming Yang --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: sample.txt Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050316/abfdb85e/sample.txt

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the

Hello, I'm using: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.1 year 2003 month 11 day 21 language R and I get the following error with: library(Design) df <- list(pre=c(0,, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1,

Hello, ? I believe that in R, it is not possible to analyze mixed effect-models when the distribucion is not gaussian (p.e. binomial or poisson), isn't? ? Somebody can suggest me alternative? ? thanks ? xavi ? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I have undertaken a small project of writing a Vorbis decoder completely from the spec, with the goal of catching any errors, both typographically and algorithmically, in the spec. My "reference" decoder is also being written in an extremely methodical style, with practically no algorithmic optimization, such that it will be clear that every step has been copied exactly from the spec.

Dear all, I have a large vector (lets call it myVector) and I want to plot its value with the logic below yaxis<-myVector[1] yaxis<-c(xaxis,mean(myvector[2:3]) yaxis<-c(xaxis,mean(myvector[4:8]) yaxis<c(xaxis,mean(myvector[9:16]) yaxis<c(xaxis,mean(myvector[17:32]) this has to stop when the new ..... yaxis<c(xaxis,mean(myvector[1024:2048]) will not find the correspondent number

Dear All,   I have difficulty in installing the package of Rcplex. Its description web is at http://cran.r-project.org/web/packages/Rcplex/index.html   First, I tried -------------------------------------------------------------------------------------------------- url <- "http://cran.r-project.org/src/contrib/Rcplex_0.3-0.tar.gz"; download.file(url, basename(url),

Hi There, i want to stream some pictures out of my application using the Theora codec. I can stream to a socket or a file. The current work based on the png2theora example. Something while enconding went wrong: I can see the whole video in VLC, but the complete duration of the file will not be displayed. If i play the video in Firefox, the first 2 seconds will not be played. I think the

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Is there a particular reason why CELT grabs 100KB of stack immediately? Is that really required or can that be trimmed down some/a lot? -a

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip