similar to: contrasts in lm and lme

On Thu, Jan 31, 2013 at 10:29 AM, Weijia Wang <wwang.nyu@gmail.com> wrote: > Hi, > > > > I have a new question about subsetting in R. > > > > Say we have this data frame: > > > > PT_ID Blood_Pressure OBS_TYPE > > 92 1900 90.0 DBP > > 94 1900 90.0 DBP > > 174 2900 140.0 SBP > > 176 2900

Dear colleagues, I had a question wrt the car package. How do I evaluate whether a simpler multivariate regression model is adequate? For instance, I do the following: ami <- read.table(file = "http://www.public.iastate.edu/~maitra/stat501/datasets/amitriptyline.dat", col.names=c("TCAD", "drug", "gender", "antidepressant","PR",

Tian It appears the attachment might not have worked so I'll embed Bill's message at the end. Peter Alspach > -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Peter Alspach > Sent: Thursday, 17 August 2006 8:02 a.m. > To: T Mu; R-Help > Subject: Re: [R] confusing about contrasts concept

Dear R help, I have the following data frame: structure(list(prochi = c("ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1"), date_1st_event = structure(c(14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784 ), class = "Date"),

I try to move from SPSS to R/S and am trying to reproduce the results of SPSS in R. I calculated a one-way anova with "spk" as experimental factor and erp as depended variable. The result of the Anova are the same concearning the mean square, F and p values. But I also wanted to caculate the contr.sdif(4) contrast on spk. The results are completely different now. I hope anybody can

When I use the postscript with the ComputerModernItalic family, minus signs are not properly displayed, but are changed into capital gammas. The attached test-postscript.pdf file illustrates the problem, which I created starting with the following R code: > postscript("test-postscript.eps", family = "ComputerModernItalic") > plot(seq(-5,5,1),seq(-5,5,1)) > dev.off()

Dear ecology fellows, I tried to implement Morisita-Horn distance (instead of Bray that is in the current version) in the code for the Simper analysis in vegan. I would be very grateful if someone can check if the code is right. function (comm, group, ...) { if (any(rowSums(comm, na.rm = TRUE) == 0)) warning("you have empty rows: results may be meaningless")

Hello, I have a problem concerning estimation of GLMM. I used methods from 3 different packages (see program). I would expect similar results for glmm and glmmML. The result differ in the estimated standard errors, however. I compared the results to MASS, 4th ed., p. 297. The results from glmmML resemble the given result for 'Numerical integration', but glmm output differs. For the

I am trying to construct a graph of 6 box plots of blood pressures. I want them to be on a single set of axis and I want the SBP to be ontop of the DBP. I have an array bp with the data in it and I tried a[1,]<-c(145,60,147,62,140,57) a[2,]<-c(160,75,160,74,160,70) a[3,]<-c(140,55,140,65,142,55) boxplot(data.frame(a), main = "Blood Pressures", at=c(1,1,2,2,3,3),

Looking for help kind of in a hurry. I've been searching google but not finding any options. Is there any way to fix missing /dev paths to luns without rebooting? For example, see the output from lsscsi below. The only way I know to fix this is with a reboot, but I REALLY Need to avoid that if possible. Thanks James [2:0:1:150] disk DataCore Virtual Disk DCS - [2:0:1:151]

Hi, I am using a function 'polygenic' in the package called 'GenABEL'. -------------------------------------------------------------------------------- > h2.gs <- polygenic(b.dbp~age+age2+age3+sex+b.bmi, kin=gs.gkin, data=gs) ------------------------------------------------------------------------------- ## and I got the following error message: Error in svd(X): infinite or

I tried to use R version of package I noticed the original MatLab Pckage is much better documented. I had a look at the R demo code "gtm_demo" and found that variable Y is used in advanced of being created: I wrote my own few lines as follows: inDir <- "C:/Documents and Settings/Monville/Alanine Dipeptide/DBP1/DHA" setwd(inDir) T <-

Sorry for the late reply. Thanks, Eli. That was really helpful. On Fri, Apr 12, 2019 at 9:31 PM Eli Friedman <efriedma at quicinc.com> wrote: > “--entry main” is not going to do what you want; the entry point for a > Linux program is not equivalent to the C “main”. > http://dbp-consulting.com/tutorials/debugging/linuxProgramStartup.html > has a description of how a C program

Dear R-devel'ers Below is an outline for a set of routines to improve support for multivariate linear models and "classical" repeated measurements analysis. Nothing has been coded yet, so everything is subject to change as loose ideas get confronted by the harsh realities of programming. Comments are welcome. They might even influence the implementation... -pd General

Hi R-users I spend a lot of time searching on the web but I didn?t found a clear answer. I have some doubts with 'frailty' function of 'survival' package. The following model with the function R ?coxph? was fitted: modx <- coxph(Surv(to_stroke, stroke) ~ age + sbp + dbp + sex + frailty(center,distribution = "gamma", method='aic'), data=datax) Then I get

Hi folks. I am trying to get dBpoweramp 13.5 to work on wine 1.2. It doesn't recognise my cdrom (/dev/sg3). I tried pretty much everything to get it working. I tried it also on different machines. Here is the discussion and a more detailed description what I tried so far at the dBpoweramp forum: http://forum.dbpoweramp.com/showthread.php?p=101496#post101496 EAC ( in ASPI mode) is working

Dear All, I am going through a worked example provided by Harrell, Lee and Mark (1996, Stats in Medicine, 15, 361-387). I know that the code provided is for S-PLUS and R but the languages don't differ enough for this to be a problem. I am using the Hmisc and Design libraries and have used the following code (as shown in the example provided in the referenced paper): '%in%' <-

Full_Name: Kem Phillips Version: 2.7.1 (2008-06-23) OS: Windows Xp professional Submission from: (NULL) (98.221.200.108) The Hmisc function describe fails, giving the error message: Error in formatDateTime(dd, atx, !timeUsed) : could not find function "format.dates" Loading the chron package, where function dates apparently resides, does not fix the problem. Note

Hi, I have been wrestling with ZFS issues since yesterday when one of my disks sort of died. After much wrestling with "zpool replace" I managed to get the new disk in and got the pool to resilver, but since then I have one error left that I can''t clear: pool: data state: ONLINE status: One or more devices has experienced an error resulting in data corruption.

odfWeave is throwing up a cryptic error, after successfully processing a file with a large number of figures (~30) and many \Sexpr{} calls. The error is (at least to me) cryptic. I am not sure where to look to correct this problem. Here's the error: Post-processing the contents Error in .Call("RS_XML_Parse", file, handlers, endElementHandlers, as.logical(addContext), :