search for: myatt

All, I have updated my practical introduction to using R to analyse epidemiological data for the March 2001 Computer Software in Epidemiology course at the Nordic School of Public Health to include a short section on object orientation. You can get this from: http://www.myatt.demon.co.uk As a ZIP file containing MS Word '95 (DOC) and Adobe Acrobat v4.xx (PDF) version of the document as well as the sample data files used in the exercises. I place this material in the public domain. Feel free to make use of it as you will. Mark Myatt -- Mark Myatt -.-.-.-.-.-....

...Mean) * NSU chi2 <- sum(((as.vector(tab) - exp)^2) / exp) tab exp chi2 I have some small expected values and so need to collapse some categories. If anyone has code to do that already ... but my question is how would I do the same for a negative binomial distribution? Mark -- Mark Myatt -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 sta...

Dear all, Is there a function to generate a Poisson cluster process? I believe that S-plus has a function to generate some point processes, make.pattern(), which include a Poisson cluster process. Thanks in advance ------------------------ Takashi Mizuno zoono at sci.osaka-cu.ac.jp Plant Ecology Lab. Osaka City University ------------------------

This came up on r-help but indicates a bug. dnbinom(x,n,p) calls dbinom_raw(n-1,...) which returns 0 for n<1. -thomas ---------- Forwarded message ---------- Date: Thu, 08 Feb 2001 17:10:23 +0000 From: Yudi Pawitan <yudi@stat.ucc.ie> To: Mark Myatt <mark@myatt.demon.co.uk> Cc: R-Help <r-help@stat.math.ethz.ch> Subject: Re: [R] Goodness of fit to Poisson / NegBinomial Sorry, Mark, the program worked for Rw1011, but this is what happened in Rw1021: > dnbinom(0:5,.9,.4) [1] 0.4383833 0.0000000 0.0000000 0.0000000 0.0000000 0.0...

Hi, I have a set of script like this: ---------------------- summary(data) plot(time, users, type="o", xlab ="Time", ylab="Number of Users") ---------------------- I type "rterm --slave < script.r > output.doc" and the graphic doesn't actually in the output.doc file. How to make the graphic embedded in Word or Excel file? Thanks, Yu-Ling Wu

HI there, Well, I am stuck again. How do I turn off axis annotation? Oh, and thanks to Peter Dalgaard for his help with my last question. I have an observation that others there may find interesting. In the windows implementation of R, the tck option of par() does not seem to operate as the reference manual indicates. Apparently the length of ticks using this option is supposed to be

...g code here cat("Error", "\n") } } } I tested this by steering R from a simple console (e.g. using something like: echo "<r-command> > TESTFIFO This works OK but where has stdout gone? Mark -- Mark Myatt -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 sta...

A few days ago I loaded R version 1.4.1 to run on Windows 98 and haven't been able to find information on the GUI approach in the FAQs, R Reference Index, or Search Engine . Edit has an item that allows changing GUI preferences, and the Target path in the icon ends in Rgui.exe, so it appears that R for Windows does include using the GUI approach. In S-Plus I wrote a application with

For some mysterious reason I do not believe this was delivered to the list. Sorry if it turns out to be a duplicate.---Bob >> Subject: optimal R windows machine I have the happy task of purchasing a new windows desktop PC. I have a limited but not unreasonable budget and access to discounts and so forth. What would be folks' suggestions of priorities and requirements be if I will

...a function to compute the Odds Ratio and its confidence intervall, from >a logistic model (glm(.......,family=binomial....). I've written my own, but >certainly someone did a better job. I show a simple function to do this in my introductory notes available from: http://www.myatt.demon.co.uk basically it is: lreg.or <- function(model) { lreg.coeffs <- coef(summary(salex.lreg)) lci <- exp(lreg.coeffs[ ,1] - 1.96 * lreg.coeffs[ ,2]) or <- exp(lreg.coeffs[ ,1]) uci <- exp(lreg.coeffs[ ,1] + 1.96 * lre...

> Ok...So I actually *typed* the table with WordPad, using the Tab button >in order to separate the columns, and saved it in rich text format. As you >might have guessed, it didn´t work. I got a: > Error in count.fields(file, sep, quote, skip, blank.lines.skip) : > string terminated by newline or EOF RTF is not what is meant by text. What I am writing is text. Your

Hi all. A simple question. Is there a function to compute the Odds Ratio and its confidence intervall, from a logistic model (glm(.......,family=binomial....). I've written my own, but certainly someone did a better job. Thank you in advance, Stefano *********************************************** Stefano Calza Istituto di Statistica Medica e Biometria Universit? degli Studi di Milano Via

On 3 Apr 01,, R-help Digest wrote (re: R-help Digest V2 #382): > - ----- Original Message ----- > From: "Mark Myatt" <mark at myatt.demon.co.uk> > To: "Greg Trafton" <trafton at itd.nrl.navy.mil> > Cc: <r-help at stat.math.ethz.ch> > Sent: Monday, April 02, 2001 9:31 PM > Subject: Re: [R] PICT output? > > > > 2) Is there a way to get an R plot to come ou...

There is now an updated version of "Using R for Data Analysis & Graphics" on CRAN, under Documentation | Contributed. Data sets and scripts are available from my web page, to which there is a link on the relevant CRAN page. John Maindonald email: john.maindonald at anu.edu.au phone : (6125)3473 fax : (6125)5549 Statistical Consulting Unit, Room 1194,

...for each combination x and y and then add z to that. I think that is what you want. Use: lset(col.whitebg()) Before plots to get rid of the grey background ... better, IMHO) I think the problem posting to R-Help might be down to your sending an attachment. Best wishes, Mark -- Mark Myatt -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 sta...

...b) can handle cluster analysis with 2000 records (mclust, preferred, although I'd take kmeans if nec'y). Splus 5.0 on the Manchester Cray has been falling over on this, much to my great annoyance, no matter how large a workspace I give mclust, hclust or kmeans, etc. mark at myatt.demon.co.uk noted: >I don't think that there are any "graphic- >interface-only" versions of SPSS. well, what we have (seems to be version 6.14) is definitely pull-down menus, the lot -- like using Micro$oft Excel -- but slower than treacle on an Artic day to run. >Of cou...

...t stat.ufl.edu> # tot <- sum(tab) rp <- apply(tab, 1, sum) / tot cp <- apply(tab, 2, sum) / tot mu <- tot * ( rp %o% cp ) norm(tab - mu) / sqrt( mu * ( (1 - rp) %o% (1 - cp) ) ) } "twoWay" <- function(x=NA, y=NA, userDefined=NA ) { # adapted from Mark Myatt <mark at myatt.demon.co.uk> if (is.na(userDefined)){ result <- chisq.test(table(x,y)) } else{ result <- chisq.test(userDefined) } # print(result) obs <- result$observed exp <- result$expected chi.table <- ((obs - exp)^2)/exp row.sum <- apply(obs,1,sum...

I have been inputting my data into text files using a text editor on a Windows 2000 PC. I would like to create a front end form that does validation on the fly to reduce the errors in my input and make it possible for family members that are not into computers to help me out with data input tasks. I assume that my form program would translate the input into a text file which I would then read

Dear All, Does anyone know of any GIS packages for R? I know that there are GIS facilities for S-Plus (eg Arc Info), and am looking for freeware equivalents. Best wishes, Sam Brown ----------------------------------------------------------- Sam Brown sam at zoo.cam.ac.uk Department of Zoology tel +44 1223 330933 Downing Street fax +44 1223 334466

I performed an aov() analysis and got the following results: Df Sum Sq Mean Sq F value Pr(>F) block 1 0.0040 0.0040 0.3282 0.5672 Residuals 269 3.2766 0.0122 Can anyone tell me how to extract the F value column and Pr(>F) column from the summary output of aov analysis? Many thanks in advance, Liqing Eco. Evol. Biol. UCIrvine