similar to: glmm multinomial?

There's something (probably a lot) missing in my understanding of plotmath. The LaTeX code for what I'm trying to produce is \log(-\log(\hat R)) My (probably hopeless) attempt at it uses nested group()s: expression(plain(log) * group("(", -plain(log)*group("(", \hat(R), ")"), ")")) Can anyone do it right for me, please? Thanks, Ted. Dr E.A.

I'll be teaching a graduate course on (the analysis of) experimental designs next year, using R. Does anyone know of a suitable textbook? (Venables and Ripley MASS Ch6 on Linear Models covers roughly the right material, but at a level that is way too difficult for my students.) Sorry if this has been asked before. I've looked through the archives and found Julian Faraway's book, but

Andy, Thanks. Note that DeVeaux's book, at www.williams.edu/Mathematics/rdeveaux/book/book.html, is marked "For use only by students of Math 244 Spring 2001, Williams College" and "DRAFT: do not distribute or copy". Ted. Dr E.A. Catchpole School of Maths & Stats University College, UNSW Institute of Maths & Stats Australian Defence Force Academy

There appears to be a bug in the R xfig() driver. When I run a simple example, eg > data(InsectSprays) > boxplot(count ~ spray, data = InsectSprays, col = "lightgray") the boxplot is fine. Doing the same thing after > xfig(file='test.fig') and then opening in Xfig (ver3.2 patchlevel 2) on my Linux box produces boxplots where the boxes, although correctly shaded in,

There appears to be a bug in the R xfig() driver. When I run a simple example, eg > data(InsectSprays) > boxplot(count ~ spray, data = InsectSprays, col = "lightgray") the boxplot is fine. Doing the same thing after > xfig(file='test.fig') and then opening in Xfig (ver3.2 patchlevel 2) on my Linux box produces boxplots where the boxes, although correctly shaded in,

I would like to do a logistic multiple regression on a binary variate, with what Genstat calls an offset, ie a regression variable with the coefficient constrained=1. I'm far from home and have Splus but not R here, but any information from either source would be very welcome. I'd hate to have to go back to using Genstat . . . Thanks, Ted. --------------- Dr E.A. (Ted) Catchpole

I have downloaded and tried to install JGR - Java GUI for R - Version 1.5. I followed all the instructions (at least in my best resolution) and installed Java latest version. R installation complete fine. Has anyone succeded in installing JGR on SuSE 10.3 ? In the following I hade detailed the failed installation, As root /i did the following: linux-Mimin:/usr/local/bin # ./R R version 2.6.0

Dear useRs, ================================================= Search for your R routines directly from Firefox! >> http://addictedtor.free.fr/rsitesearch ================================================= Mango Solutions, providers of R and S-Plus consulting, application development and training, are happy to announce the first release of the R Site Search extension for Mozilla Firefox.

Hi All, Here 's what I got using dgamma function : > nu<-.2 > nu*log(nu)-log(gamma(nu))+(nu-1)*log(1)-nu*(1) [1] -2.045951 > dgamma(1,nu,nu,1) [1] 0.0801333 > dgamma(1,nu,nu,0) [1] NaN Warning message: NaNs produced in: dgamma(x, shape, scale, log) Could anyone tell me what is wrong here ? I am using R-2.4.1 on windows XP. Thanks a lot.

"&&" seems to behave strangely and gives different results from "&" e.g. in a data frame selection (regardless whether terms are bracketed)? ===========Script======================= test=data.frame(gender=c("F","M","M","F","F"),side=c("R","L","R","L","R")) test

Dear All, Is there a conmand to calculate the multiplication of the elements in a vector? For example: a=c(1,2,3,4) I want to get 1*2*3*4=24. Because the dimension of the vector is high, I want to know there is a command for this task in R? Thank you so much! Kang __________________________________________________ [[alternative HTML version deleted]]

Hello Listers, I am trying to sample a vector to create a new one of sample length, witha sum equal to the sum of the initial vector: initial = 10, 30, 10 (sum=50) sample example = 5, 35, 10 (sum=50) or 25, 15, 10 (sum=50), etc ... My problem is to control the sum, so it stays constant. Any suggestions would be very helpful ! Thank you in advance, Eric

Dear List, I want to let R calculate the time (run-time) it requires to run a self-written simulation function. I tried as follows: it enables me to see the starting and finishing time points. ############################# "sim.result" <- function(nsim, ...){ Starting <- date() ... # calculations # final.result <- ... # Output for display # cat("# of Iterations used

Hi all, I have a problem on how R calculates Percentiles : Suppose I have following data set: > data1 [1] -16648185 -14463457 -14449400 -13905309 -13443436 -13234755 -12956282 -11660896 [9] -10061040 -9805005 -9789583 -9754642 -9562164 -9391709 -9212182 -9151073 [17] -9092732 -9068214 -8978151 -8943912 -8761890 -8632106 -8541580 -8501249 [25] -8234466 -8219015

Dear R-users, I haven?t found rules for logical operators. I need to select data according the following rule: Condition A & (Condition B | Condition C) How should I write it?? Is Condition A & Condition B | Condition C correct or will R execute (Condition A & Condition B) | Condition C ? Thanks for your help. Delphine Fontaine Delphine Fontaine Statistician Statistics

Hi there, I want to integrate a function over an irregular polygon. Is there any function which can implement this easily? Otherwise, I am thinking of divide the polygon into very small rectangles and use "adapt" to approximate it. Do you have any suggestions to get the fine division? Any advice is appreciated. Haiyong

Despite lots of investigation, I haven't found any R packages might be suitable for the following problem. I'd be very grateful for suggestions. I have three-way nested data, with a series of measures (obs) taken in quick succession (equal time spacing) from each subject on different days. The measures taken on the same day are temporally correlated, so I'd like to use an AR1

Hi again! Well, I just realised that I am running an *unstable* version of R (if such exists) from early January, that I picked up for some weird reason (the appearance of manova??). In fact the starting message is R : Copyright 2001, The R Development Core Team Version 1.3.0 Under development (unstable) (2001-01-22) I have been tried to reproduce the same example in R 1.2.3 for Windows

I am trying to run a GLMM on some binomial data. My fixed factors include 2 dichotomous variables, day, and distance. When I run the model: modelA<-glmmPQL(Leaving~Trial*Day*Dist,random=~1|Indiv,family="binomial") I get the error: iteration 1 Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 >From looking at previous help

Dear list members, In re-running with GLMM() from the lme4 package a generalized-linear mixed model that I had previously fit with glmmPQL() from MASS, I'm getting a warning of a convergence failure, even when I set the method argument of GLMM() to "PQL": > bang.mod.1 <- glmmPQL(contraception ~ as.factor(children) + cage + urban, + random=~as.factor(children) + cage +