similar to: code for mixed model in R?

We just studied randomized block design analysis in my statistics class, and I'm trying to learn how to do them in R. I'm trying to duplicate a case study example from my textbook [1]: > # Case Study 13.2.1, page 778 > cd <- c(8, 11, 9, 16, 24) > dp <- c(2, 1, 12, 11, 19) > lm <- c(-2, 0, 6, 2, 11) > table <- data.frame(Block=LETTERS[1:5], "Score

Hi everyone, I have the following problem: I have three variables, 'group', 'city' and 'pressure' There is an interaction effect between group and city and I'd like to show this in an interaction plot: interaction.plot(group, city, pressure, type="b", col= c(1:2), leg.bty="o", leg.bg="blue", lwd=1,

Dear R-helpers: I am trying to fit a multivariate Cox proportional hazards model, modelling survival outcome as a function of treatment and receptor status. The data look like below: # structure of the data str(sample.data) List of 4 $ survobj : Surv [1:129, 1:2] 0.8925+ 1.8836+ 2.1191+ 5.3744+ 1.6099+ 5.2567 0.2081+ 0.2108+ 0.2683+ 0.4873+ ... ..- attr(*, "dimnames")=List of 2

Hi, I am doing a retrospective analysis on a cohort from a designed trial, and I am fitting the model fit<-glmD(survived ~ Covariate*Therapy + confounder,myDat,X=TRUE, Y=TRUE, family=binomial()) My covariate has three levels ("A","B" and "C") and therapy has two (treated and control), confounder is a continuous variable. Also patients were randomized to

Hi all, I dumped my schema from SQLServer 2000 without any problems. I then tried to "rake db_schema_import" on MySQL 4. To do this i changed my database.yml settings to use the mysql adapter (as it was set to use sqlsever before) and attempted the import. I haven''t changed a thing in the schema.rb script and, even though I do have generated a migration, I think I

Hi, >From idea to launch in less than 48 hours (those 48 hours including the new year party ;) , rails once again showed its effectiveness. Last friday I read the news about a guy who sold a million pixels for 1$ each on his homepage. After slamming my head during half an hour for not having thought of it [1], I thought I''d better spend some time on a rails project, and more as a

Is there anything in the Puppet/Passenger setup process that edits httpd.conf and sets the User/Group running apache to "puppet" ? A master I set up several months ago is configured that way. I do not recall doing it. But then the high voltage electroshock therapy does have its side effects :) “Sometimes I think the surest sign that intelligent life exists elsewhere in the

Hello, I need help in performing a Van_der_Waerden normal scores test in R. I have two arrays of scores(final on therapy scores from drug and placebo) and want to use the normal scores procdeure to test for significance. (observations are unequal in number - due to dropouts). Could you please help me out with the coding or let me know if there is a package that can be used (for example,

Hi, I am wondering if there is a way to change the value of the "record" option in a graphics device that is already open (and accepts this option). I don''t want to open a new device with, for example "dev.new(record=T)", but just want to change the settings of the current device. This can be done by pointing and clicking on the "history" tab of a

Hello Everyone, I have what I think is a complex text parsing task. I've provided some sample data below. There's a relatively simple version of the coding that needs to be done and a more complex version. If someone could help me out with either version, I'd greatly appreciate it. Here are my sample data. haveData <- structure(list(profile_key = structure(c(1L, 1L, 2L, 2L, 2L,

Hi, Can somebody advice on weird mysqlWriteTable bug. > mysqlWriteTable(conn, 'comparison',design2, row.names = F, overwrite=T) Error in mysqlExecStatement(conn, statement, ...) : RS-DBI driver: (could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near '"condition"

I'm using R on Windows v2.0.1 with the nlme package (v3.1-53) and am finding some unexpected discrepancies in the output of intervals.lme and coef.lme. I've included a toy dataset at the end, but briefly, the data are longitudinal data from couples in marital therapy. Each spouse's relationship satisfaction is measured 4 times; I've fit both linear and quadratic models to the

Hello Arnaud, On 05/04/14 22:01, Arnaud Quette wrote: > Philippe, Ivan, > > please see https://github.com/networkupstools/nut/issues/118 > and send back your data. I've pasted the requested files in github, but the auto-formatting may have damaged them, so here they are. Also, when I ran gen-snmp-subdriver.sh, it gave me the following message:

Philippe, Ivan, please see https://github.com/networkupstools/nut/issues/118 and send back your data. cheers Arno -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.alioth.debian.org/pipermail/nut-upsuser/attachments/20140405/682f10dd/attachment-0001.html>

Hi everyone, I just landed into the theora planet, as a game programmer, I searched for a free video fomat/codec and the theora choice became obvious. However I experienced rather bad performance (at least from a game programming point of view) After a couple a profiling, I discovered, as previous discused in a post found via Google, that the bottleneck is in the ogg library. An unsane part of the

Dear All, I am a bit struggling with the many packages for Markov models available in R. Apologies for now posting a code snippet, but I am looking for some guidance here. Please consider a set like the one below (which you can get with data<-read.csv('http://dl.dropboxusercontent.com/u/5685598/data_table.csv') ). ID therapy age1 age2 EFS 7308 ormo_lunga 78

I'll take a risk following Uwe's wonderful response to that advert. I'm looking for someone, perhaps particularly a stats student, who might want to do a piece of work with me on using R to present and analyse routine data that psychotherapists might submit on a cgi-bin interface and perhaps cross-referencing that against some largish referential data to which I have access.

COMPUTATIONAL ANALYSIS OF NEXT GENERATION SEQUENCE DATA (http://bioinformatics.nki.nl/vacancies.php) THREE BIOINFORMATICS POSITIONS PROJECT OUTLINE Genomic alterations are major determinant of responses to (targeted) therapies in cancer. In fact, the best positive and negative predictors of responses to targeted therapies are alterations in kinases or their direct downstream effectors. To

Lets say you had a pretty decent job, but your company mandated use of emacs. The poaching employer offers you a free copy of Textmate(tm). Would you jump ship solely in order to use Textmate at company B? Cheers, Dave -- Posted via http://www.ruby-forum.com/.

I''ve had a couple ideas for Ruby on Rails shirt graphics for a little while now: 1. "Fuck You" superimposed over, oh, probably a Ruby gem. Or Rails of some sort. 2. "Ruby on Rails" superimposed over rails of cocaine. Maybe also with a graphic of a crazed looking coder. Anybody got graphics? Or can make them? Heh! Joe -- Posted via http://www.ruby-forum.com/.