similar to: Unexpected behaviour of read.table

Dear R-users I am introducing my collegue to R and she is interested to use R to perform Polytomous Logistic Regression called also Multi-category logistic regression. Is there any program in R doing this? Thank you so much for any feedback. I take this opportunity to thank people who were kind to answer my query about how to include a postscript file generated from R, in a latex file and

Has anyone tried to build R on Suns by linking to the versions of BLAS and LINPACK in the Sun Performance Library rather than the standard versions included in the distribution? On the face of it this would seem to be a highly desirable thing to do since these versions of the routines are more efficient and can exploit parallelism on multiprocessor machines. David Clayton MRC Biostatistics

Hello folks, Any ideas how to do this? data.frame is a data frame with column names "x1",...,"xn" y is a response variable of length dim(data.frame)[1] I want to write a function function(y, data.frame){ lm(y~x1+...+xn) } This would be easy if n was always the same. If n is arbitrary how could I feed the x1+...+xn terms into lm(response~terms)? Thanks Richard -- Dr.

Dear everybody! Are there algorithms for multinominal logit/probit available for R? Is it my fault that I cannot find these in CRAN? Has somebody programmed these? with best wishes Ott Toomet Ott.Toomet at mail.ee -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

> Date: Tue, 7 Nov 2000 17:34:37 +0000 (GMT) > From: Jonathan Myles <mylesj at icrf.icnet.uk> > > I've now got a colleague at MRC-BSU to compress it, and email it as an > attachment to my UNIX machine, where I uncompressed it and then ftp'ed > it in ascii format to my laptop. So it doesn't seem to be an ascii/bin > issue ... Is it possible to get a binary

Dear Wine-friends, I was wondering if any of you would have a clue around why WinBUGS (http://www.mrc-bsu.cam.ac.uk/bugs/), a nifty Markov chain Monte Carlo sampler widely used in Bayesian statistical modelling, no longer works when run through any version of WINE newer than 1.1.12. I have encountered this issue of my machines at home (which runs on Zenwalk 6.0) and work (Mandriva 2008.1) a

Is there any way of running R in ESS using -n and -v arguments other than the default? Thanks Jonathan -- Dr. Jonathan Myles e-mail:jonathan.myles at mrc-bsu.cam.ac.uk MRC Biostatistics Unit Tel. 01223 330371 Institute of Public Health FAX 01223 330388 University Forvie Site Robinson Way CAMBRIDGE CB2 2SR

Hi Just one of those niggles....... I've just been trying to plot a filled circle. I thought that this would do it plot(1,1,type="n") points(1,1,pch=1,bg="blue",cex=5) #bg: background ("fill") color for open plot symbols But I need to do this instead points(1.2,1,pch=19,col="blue",cex=5) Am I misunderstanding the "bg" option in the points

Is there any way of getting just the LaTex output which R CMD Rd2dvi generates on the way to produces a .dvi file, which could then be incoorporated into another document? Thanks, Jonathan -- Dr. Jonathan Myles e-mail:jonathan.myles at mrc-bsu.cam.ac.uk MRC Biostatistics Unit Tel. 01223 330372 Institute of Public Health FAX 01223 330388 University Forvie Site

Hello, I have found a problem with using the glm function with a gamma family. I have a vector of data, assumed to be generated by a gamma distribution. The parameters of this gamma distribution are estimated in two ways (i) using the glm() function, (ii) "by hand", using the optim() function. I find that the -2*likelihood at the maximum found by (i) is substantially larger than that

Hello, I've had a quick look on the list and can't find an answer to this niggle. Whilst debugging some code, I noticed that looping variables in R are not local to the loop as seams common in procedural languages. For example consider the following piece of code: for(i in 1:3){ cat(i,"") for(i in c("a","b","c")){ cat(i,"")

Hello, In WinBUGS 1.4 manual (http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/manual14.pdf), the gamma density is presented as dgamma(r,mu) where r and mu are the shape and rate parameters, respectively. In JAGS (rjags) manual version 1.0.2, May 9, 2008 (http://www-fis.iarc.fr/~martyn/software/jags/jags_user_manual.pdf), on page 26 the gamma density is presented as dgamma(mu,r) instead of dgamma(r,mu).

Hi I'm using kronecker() with a matrix and a vector. I'm interested in the column names that kronecker() returns: > a <- matrix(1:9,3,3) > rownames(a) <- letters[1:3] > colnames(a) <- LETTERS[1:3] > b <- c(x=1,y=2) > kronecker(a,b,make.dimnames=TRUE) A: B: C: a:x 1 4 7 a:y 2 8 14 b:x 2 5 8 b:y 4 10 16 c:x 3 6 9 c:y 6 12 18 > The

Hi I am trying to use the deriv and eval functions to obtain the value of a function , say "xi-(alpha0+alpha1*gi)" , differentiated with respect to alpha0 and alpha1, in the following way # for gi = 0 > dU1dtheta <- deriv(~ xi-(alpha0+alpha1*gi), c("alpha0","alpha1")) > eval(dU1dtheta) (Intercept) -0.2547153 attr(,"gradient")

Hi Does anyone know of any R coding/functions for MCMC approaches? I am currently using BUGS but I wonder if the bazaar has produced anything? I think I am pushing BUGS to it's limit and possibly past it at the moment. John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Thanks to Vin Everett and David Clayton <david.clayton at mrc-bsu.cam.ac.uk>, an RPM package of R-0.62.4 for Red Hat Linux/SPARC 5.1 is now available on CRAN. This means that RPMs are available for all three platforms (Intel, Alpha, SPARC) currently supported by Red Hat. Martyn -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-announce mailing list -- Read

Thanks to Vin Everett and David Clayton <david.clayton at mrc-bsu.cam.ac.uk>, an RPM package of R-0.62.4 for Red Hat Linux/SPARC 5.1 is now available on CRAN. This means that RPMs are available for all three platforms (Intel, Alpha, SPARC) currently supported by Red Hat. Martyn -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-announce mailing list -- Read

I want to cite R for an article. What is the right thing to do? Shall I include "Notes on R" or is there a better way? -- mailto:koloska at rcs.urz.tu-dresden.de http://rcswww.urz.tu-dresden.de/~koloska/ -- -- right now the web page is in german only but this will change as time goes by ;-)

Dear R users, I'd like to announce a new package on CRAN called ``denstrip''. It implements ``density strips'' and other graphical methods for illustrating and comparing distributions in a compact fashion. Posterior distributions of parameters are often summarised using point and line drawings of means and credible intervals. This is common, for example, in multiple

Dear R users, I'd like to announce a new package on CRAN called ``denstrip''. It implements ``density strips'' and other graphical methods for illustrating and comparing distributions in a compact fashion. Posterior distributions of parameters are often summarised using point and line drawings of means and credible intervals. This is common, for example, in multiple