similar to: lme, diagnostic plots and missing values

Dear all, This question is partly statistics and partly R and I apologise in advance for my (usual) verbosity! My data is a little more complicated that this suggests, but essentially I have proportion data from different studies (id), each from a specific country and region of the World. I would like to examine the variables that affect the proportion, but these factors are hierarchical. In case

Dear All, Having now successfully started using Sweave, I have just noticed an odd side effect with fonts. Using plain LaTeX, I have had no problems using \usepackage{palatino} or \usepackage{times} and the font correctly changes. However, once I convert the tex document into Snw, and run it through R CMD SWEAVE, LaTeX seems to ignore this package directive. The LaTeX log suggests the font

Hello all, I've been able to use R very successfully to run simple statistics and generate the plots I require. I've been evaluating Sweave, and have hit upon a small problem that I don't seem to be able to workaround. Sweave runs very well for single file latex documents, but I have a complex thesis made up of several parts and chapters. These are arranged with a master latex file

Dear all, I haven't seen this mentioned and presumably Dr. Deepayan Sarkar is too polite to advertise! I just received a flyer from Springer: A new book on Lattice released today!! http://www.springer.com/statistics/computational/book/978-0-387-75968-5?cm_mmc=NBA-_-Apr-08_UK_1753460-_-product-_-978-0-387-75968-5

It seems that Access needs that you surround the dates with a # symbol. You probably need something like. res <- sqlQuery(channel, "select * from test_table where market = 'atl-bos' and competitor = 'delta' and dd = #2007-11-20#") Hope this helps, Stefano -----Messaggio originale----- Da: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]Per

Since we're on the topic of book reviews, I just received Phil Spector's new R book called "Data Manipulation with R" and it is also quite a nice book. I haven't gone through it all and I won't give a detailed review but I have gotten a lot out of the first 100 pages that I have read. Note that I've been using R for almost 1.5 years so , for me, it's a

Dear all, I'm struggling with a plot and would value any help! I'm attempting to highlight a histogram and density plot to show a proportion of cases above a threshold value. I wanted to cross-hatch the area below the density curve. The breaks and bandwidth are deliberate integer values because of the type of data I'm looking at. I've managed to do this, but I don't think it

Dear All, Writing sooner than I thought I'd need to. I'm using R 2.4 on Mac OS X, with RODBC, PostgreSQL 8.1 and Actual's ODBC driver. I have all my data in Filemaker 8.5, but it is automatically exported into PostgreSQL for analysis as Filemaker's ODBC and JDBC access is awful, slow and has a tendency to crash. I have disability data where for each patient there is a survival

Dear all, I'm having difficulty getting access to data generated by survfit and print.survfit when they are using with a Cox model (survfit.coxph). I would like to programmatically access the median survival time for each strata together with the 95% confidence interval. I can get it on screen, but can't get to it algorithmically. I found myself examining the source of print.survfit to

Dear all, I have been stuck on this problem, am rather struggling and would appreciate some advice if anyone can help. I apologise if this is a bit long-winded, but I've tried to limit it to the bare essentials, but don't know how to make it more generic! I have some slightly odd real world data that I'm looking at representing number of positive diagnoses for different diseases,

I'm using lmer to fit mixed-effect logistic regression models. This is for a small data set. First, I fit a constant: Generalized linear mixed model fit using Laplace Formula: propm ~ (1 | study) Data: inducedSR71507.dat Family: binomial(logit link) AIC BIC logLik deviance 183.7 189.4 -89.84 179.7 Random effects: Groups Name Variance Std.Dev. study (Intercept)

Running the Chapter 4 examples in Pinheiro & Bates' "Mixed-Effects Models in S and S-PLUS" (2000), I get a message that the default optimizer doesn't converge, but using "optim" for the optimizer results in convergence: > > library(nlme) > > fm1IGF.lis <- lmList(IGF) > > fm1IGF.lme <- lme(fm1IGF.lis) > Error in lme.formula(fixed =

Dear useRs, I have written a function that implements a Bayesian method to compare a patient's score on two tasks with that of a small control group, as described in Crawford, J. and Garthwaite, P. (2007). Comparison of a single case to a control or normative sample in neuropsychology: Development of a bayesian approach. Cognitive Neuropsychology, 24(4):343?372. The function (see

Dear all, I'm still new to R, but have a fair experience with general programming. All of my data is stored in postgresql, and I have a number of R files that generate tables, results, graphs etc. These are then available to be imported into powerpoint/latex etc. I'm using version control (subversion), and as with most small projects, now have an ever increasing number of R scripts,

Dear all, I'm using the "survival" package with R 2.4.0 on Mac OS X 10.4.8. I have two core statistics books (one of which is Altman's medical stats book) which suggest showing the number of individuals at risk at different time intervals on the Kaplan-Meier curve. My plot shows two curves that later cross, because of one significant outlier. I have two queries: Is there an

I'd like to plot n (say n = 10) semi-transparent lines in such a way that if all n happen to exactly overlay each other, the resulting line is completely opaque. In principle, I believe that this is what setting alpha = 1/n should accomplish. In practice, different values of n produce different results. Consider the following function, which overlays n semi-transparent red lines: mylines

I'm not sure if the gee package is still actively maintained, but I for one find it extremely useful. However, I've come across a few infelicities that I'm hoping could be resolved for future versions. Hope it's okay to list them all in one post! They are: (1) AR(1) models don't fit when clustsize = 1 for any subject, even if some subjects have clustsize > 1. (2) If the

Hi everybody, I'm having trouble getting an eval-parse combination to work with lme in a for loop. Consider the code below. The call to lme outside the loop works. The call to aov inside the loop works, but the call to to lme inside the loop does not. The latter fails with Error in model.frame.default(formula = ~meanCol + Group + c1 + c3 + Subject, : variable lengths differ (found for

Dear list Maybe someone can help with the following problem: I have a source file containing the following code for plotting: pdf(file = 'data/mario/boxplot.pdf') bwplot(sig100 ~ target | group, data = fish, main = 'Boxplot der Signale in 100 Zellen nach Gruppe', xlab = '', ylab = 'Anzahl der Signale') dev.off() pdf(file = 'data/mario/xyplot.pdf')

Hi: I am working with someone remotely to allow them access to our data. The follow command using "download.file" works perfectly on my Mac: > > download.file(url="http://oceanwatch.pfeg.noaa.gov:8081/thredds/ > wcs/satellite/AG/ssta/14day? > request=GetCoverage&version=1.0.0&service=WCS&format=NetCDF3&coverage= >