similar to: unbalanced design

Kia ora list members: I'm having a little difficulty getting the correct standard errors from a glm.object (R 1.9.0 under Windows XP 5.1). predict() will gives standard errors of the predicted values, but I am wanting the standard errors of the mean. To clarify: Assume I have a 4x3x2 factorial with 2 complete replications (i.e. 48 observations, I've appended a dummy set of data at the

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

Dear R gurus, I have the following model: appcov.aov <- aov(yield ~ prevyield + trt + block) where prevyield is a continuous numeric covariate and trt and block are factors (yes, I did factor()!) Now, when I do a TukeyHSD, my diff's are all screwed up! For instance: treatment mean for treatmen "E" is 277.25 and for treatment "O" is 279.5, so I figure the diff O-E

Dear Speex community, I am having problems to get the Speex Decoding API to work properly in our project. I am working in an embedded project with a very limited RAM and Flash memory and would like to use Speex encoded files to save memory while having good voice sounds in our speakers. I have developed a Windows-based tool to convert WAV files to Speex files in packets of 160 bytes (we are

hi, everyone: i have a question on the reshape function. i have the following dataset : gene tissue patient1 patient2 patient3............. _________________________________________________ gene1 breast 10 20 50 gene2 breast 20 40 60 gene3 breast 100 200 300 which i hope to convert to the following format: gene patientID

Dear R users - I have a list of patient identifiers and diagnoses from inpatient admissions. I would like to reorganize the list, presently in a long format to a wide format in reshape, but in the absence of a "time" element, I am uncertain how to do this - any help greatly appreciated. ID Dx A nausea A diabetes A kidney failure A heart attack A fever B fever B

Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ######################################################> kfitm1 Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = "gauss"), data = kidney)

Hi Dear R-users I have a database with 18000 observations and 20 variables. I am running cox regression on five variables and trying to use survfit to plot the survival based on a specific variable without success. Lets say I have the following coxph: >library(survival) >fit <- coxph(Surv(futime, fustat) ~ age + rx, data = ovarian) >fit what I am trying to do is plot a survival

I have data I'd like to plot using the formula interface to boxplot. I call boxplot like so: with(mydata, boxplot(count ~ geno * tissue)) I get a boxplot with x axis labels like "wt.kidney". I would like to change the '.' to a newline. Where is this separator configured? Thanks, -Ed

Hello Everyone, I am having problem in defining specific axis for plotting a vactor. vecAVG <- c(0.2, 0.4, 0.6, 0.2, 0.4) names(vecAVG)<-c("brain","heart","kidney","lung","blood") par(mar=c(12,4.1,4.1, 2.1))

Hi, I am doing an analysis to see if these is tissue specific effects on the gene expression data . Our data were collected from 6 different labs (batch effects). lab 1 has tissue type 1 and tissue type 2, lab 2 has tissue 3, 4,5,6. The other labs has one tissue type each. The 'sample' data is as below:

Dear R-users, I have noticed small discrepencies in the reported estimate of the variance of the frailty by the print method for survreg() and the 'theta' component included in the object fit: # Examples in R-2.6.2 for Windows library(survival) # version 2.34-1 (2008-03-31) # discrepancy fit1 <- survreg(Surv(time, status) ~ rx + frailty(litter), rats) fit1 fit1$history[[1]]$theta

Dear All I am trying to do a repeated measures analysis using lmer and have a number of issues. I have non-orthogonal, unbalanced data. Count data was obtained over 10 months for three treatments, which were arranged into 6 blocks. Treatment is not nested in Block but crossed, as I originally designed an orthogonal, balanced experiment but subsequently lost a treatment from 2 blocks. My

Dear friends - We have 25 rats, 14 of these subjected to partial removal of kidney tissue, 11 to sham operation, and then followed for 6 weeks. So far we have data on 26 urine metabolites measured by NMR 7 times during the observation. I have smoothed the measurements by b.splines in fda including a roughness penalty, and inspecting the mean curves for nephrectomized and sham animals indicate

Dear all, in survival5, kidney data set, appears in help page: "survival5 does not reproduce the original analysis." What does it means? thanks in advance TCM -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20010927/e006b08d/attachment.html

It was recently pointed out on the lists that the S-PLUS Trellis suite has a function called reorder.factor that's useful in getting useful ordering of factors for graphs. I happily went ahead and implemented it, but it turns out that R (not S-PLUS) has a generic called reorder (with a method for "dendrogram"). Naturally, this causes R to think I'm defining a method for

Hello, R experts, I have a list called dataHP which has 30 elements (m1, m2, ..., m30). Each element is a 7x6 matrix holding yield data from two factors experimental design, with treatment in column, position in row. For instance, the element 20 is: dataHP[[20]] col1 col2 col3 trt1 trt2 trt3 [1,] 22.0 20.3 29.7 63.3 78.5 76.4 [2,]

I have two models called entity and user. The entities table has a column called users_id that contains the user id of the user that created the entity. In entity I have... has_one :user ... as I want to be able to show the user who created the entity from the entity object. But this produces the following error... Mysql::Error: #42S22Unknown column ''users.entity_id'' in

Hi, I have an n x m matrix of numerical observations. ie. stock prices I wish to convert the matrix of observations to a matrix of simple returns (by taking the differences between (column) observations.) Can any good soul suggest a function for this? -- View this message in context: http://r.789695.n4.nabble.com/R-beginner-matrix-algebra-tp4653335.html Sent from the R help mailing list

Hello listers, Does anybody know how to add text to an xyplot without whipping out the existing curve? That's all. For instance, Lets say you generate a graph like this A <- data.frame(x = rnorm(100), y = rnorm(100)) xyplot(y ~ x, data = A) How would you add 'Hello world'? I tried 6.02E23 different partial solutions found on the web and failed. I just need one EXAMPLE that