similar to: Why are looping variables not local?

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 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.

Hello, I've been looking in the on-line manuals and searching past posts but can't find an answer to this question. I'd like to define my own variance function in a GLM. The function glm(formula, family=quasi(var="var function")) lets me choose from a selection of built in variances, but I want to define my own function for the variance. Is there an S-plus

Dear R users: is there a package, similar to varclus in SAS or varclus in S, ported or written for R? Also, is there any other package in R that was designed for grouping the variables under different measures of distance (in cases where data is non-Gaussian, autocorrelated, and so on). Janusz. -- ** Janusz Kawczak ** ** UNC at Charlotte,

Hi, I want to perform goodness of fit test for multinominal distribution. The easiest way would be to use Chi2, but the measurment errors are not normally distributed. I thought about using some bootstraping method to perform the analysis. How can I do it in R? Tomek -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

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

Hello, I have two vectors x1 and x2 both in increasing order. I want to select the x1[j]th entry which is the max min of the x2[i]th entry. I can do this using if and for statements but is there a quick way to do it without running a loop? Thank you in advance, Pantelis

I have put a preliminary version of factanal in a package at http://www.stat.ox.ac.uk/pub/R/fa_0.1-1.tar.gz and would appreciate feedback. In particular 1) What else you might want. Please give precise references with algorithmic details, as that has been a major problem for me. 2) Some convincing examples of factor analysis. Most (literally) of the ones I have looked at are local maxima, not

I've compiled the 1.8.0 R on Solaris 2.9 and when trying to use plot command I get the Bus error and a core dump. Anybody experienced something like this? If yes, what should I patch? Janusz. [[alternative HTML version deleted]]

Dear List, I am looking for a possibility to fit a mixture model under R using maximum likelihood estimation. Venables and Ripley describe a solution working under S+ (in MASS, 3. ed., p. 263) which requires the D system function and deriv3. This solution does not seem to be portable to R or at least I do not realise how. Is there anyone who a) knows how one could make the MASS-method run under

When reading comma-delimited files as saved from a spreadsheet (unfortunately many of my scientific collaborators give me these) in read.table(), missing values are spotted most of the time. Unfortunately when comma is the first character on the line it gets it wrong. For example, reading the file 1,,3 ,5,6 ,8,9 with read.table("test.dat", header=F, sep=",") R gives

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

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

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

Is there anyone who succeeded in installing R-0.61 on a IBM RISC/6000 under AIX? Juergen. ======================================================================= Juergen Dippon Mathematisches Institut A, Universitaet Stuttgart Pfaffenwaldring 57, 70569 Stuttgart

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 ;-)

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

Has anybody ported the Classification and Regression Trees package to R? Stefano =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch

Hello, I'd like to do something like: n<-100 zz<-cbind(rnorm(n),rnorm(n),floor(runif(n)*3+1)) colors<-c("red","green","blue") plot(zz,col=colors(zz[3])) and have a matrix of scatterplots colored by class. The above does not work, of course, but I'm not sure exactly what function I'm looking for. Thank you for your time, Dave, -- Dave

Is there a way to "unplug" the windows graphics from R, so that it will work across a teletype telnet session or in a Linux box proving resistant to having windows added to it? Thanks Jonathan Myles -- 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