search for: jorgensen

...scanexamp.tex Processing code chunks ... 1 : echo term verbatim Error: chunk 1 Error in parse(text = chunk) : unexpected numeric constant in: "height = scan() 64 62" > Comments would be appreciated. (And thanks to Ross Darnell for a lot of help on another list.) Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 0200 8350

...ctor [Fac], returning a vector of the same length as vec. After a while I realised that lm1 <- lm(vec ~ Fac) fitted(lm1) did what I want. But there must be another way to do this, and it would be good to be able to apply other functions than mean() in this way. Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862

...-em-all may not be viable much longer. Has any thought been given to dividing "contributed" binaries into a recommended set, perhaps a couple of hundred, and the remained. That way one could install the recommended ones routinely and add in the others as required. Any comments? Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862

How do you stop R from putting up a dialog box when you quit Rgui? (I use Windows and I never save workspaces that way) Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk +64 7 849 6486 home Mobile 021 1395 862

...ariables x it would be helpful if the function printed "x" rather than the formal argument "dvar". I have a vague idea that things like deparse() and substitute() will come into the solution but I have not yet come up with the right incantation. Any help appreciated! Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862

I want to interlace two vectors. This I can do: > x <- 1:4 > z <- x+0.5 > as.vector(t(cbind(x,z))) [1] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 but this seems rather inelegant. Any suggestions? Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk +64 7 849 6486 home Mobile 021 1395 862

...(MASS) > data(wtloss) Warning message: Data set `wtloss' not found in: data(wtloss) > data(abbey) Warning message: Data set `abbey' not found in: data(abbey) And yet all the .rda files for the MASS datasets are in D:\Program Files\rw1022\library\MASS\data Help, anyone? Murray Jorgensen, Department of Statistics, U of Waikato, Hamilton, NZ -----[+64-7-838-4773]---------------------------[maj at waikato.ac.nz]----- "Doubt everything or believe everything:these are two equally convenient strategies. With either we dispense with the need to think." http://www.stats.waikat...

Bother! This cold has made me accident-prone. I meant to hit Reply-all. Clarification below. -------- Original Message -------- Subject: Re: [R] Parameterization puzzle Date: Fri, 21 Jul 2006 19:10:03 +1200 From: Murray Jorgensen <maj at waikato.ac.nz> To: Prof Brian Ripley <ripley at stats.ox.ac.uk> References: <44C063E5.3020703 at waikato.ac.nz> <Pine.LNX.4.64.0607210716270.12611 at gannet.stats.ox.ac.uk> Apologies for a non-selfcontained example. Here is what I should have sent: pyears <- sc...

how about this: tapply(vector,list(factor1,factor2),function)? -----Original Message----- From: owner-r-help at stat.math.ethz.ch [mailto:owner-r-help at stat.math.ethz.ch]On Behalf Of Murray Jorgensen Sent: 24. september 2000 02:46 To: R-help Subject: [R] Crosstabulation I can't seem to find a function in R similar to Splus crosstabs() for creating a multi-way table from factors and a count vector. Murray Jorgensen -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

Here is the file minimal.Snw: \documentclass[a4paper]{article} \title{R tips and tricks} \author{Murray Jorgensen} \usepackage{Sweave} \begin{document} \maketitle \section*{Entering data from a single variable} The following data are transformed tensile strength measurements on polyester fibres. They may be found on the file \texttt{TENSILE.DAT}. We may enter this data into R using the \texttt{scan} command....

I'm going to try my hand at converting some Fortran programs to R. Does anyone know of any good articles giving hints at such tasks? I will post a selective summary of my gleanings. Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 0200 8350

...lay a few polynomial regression fits I can do something like for (i in 1:6) { mod <- lm(y ~ poly(x,i)) print(summary(mod)) } Suppose that I don't want to over-write the fitted model objects, though. How do I create a list of blank fitted model objects for later use in a loop? Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862

...es=c("character",rep("numeric",7)), comment.char = "#") } I need something like an unquote() function that will allow the "brand" on the LHS of the assignment to be treated as an object name instead of a character string. Murray Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk +64 7 849 6486 home Mobile 021 1395 862

...tation of likelihood ratio and score test based confidence intervals in generalized linear models" by Juha Alho (1992) (Statistics in Medicine, 11, 923-930). Being lazy, I thought that I would ask if anyone else on the list has already done this? (or something equivalent) Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk +64 7 849 6486 home Mobile 021 1395 862

...0916 Random effects: Formula: ~1 | Subject (Intercept) Residual StdDev: 2.06847 0.7800331 [...etc...] I can extract the estimate of the variance component \sigma (0.7800331) via sigma <- fm1OrthF$sigma How do I extract the other component \sigma_b (2.06847) ? Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862

...mn: V1 145 Cahuila (United States) 862 Paipai (Mexico, USA) 863 Paipai (Mexico, USA) 864 Paipai (Mexico, USA) 311 Eyak (Alaska) Does anyone know what is going on? Tusen tak in advance, in fact 1254 tak in advance! Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 1395 862

Hai - kann mir jemand sagen, wie ich den Modalwert in R berechne?! IRgendwie finde ich den Befehl nicht.... greetz und herzlichen Dank Sonja

...C:\apps\R\R-2.11.1\doc\html\packages.html does not include all the packages that I have installed and linked. I don't want to read my html help from the web because sometimes I am off-line or on a slow connection. How do I go about getting a local set of html help files? Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz majorgensen at ihug.co.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 0200 8350

Has anyone got a cute way of compiling a table of the available plotting symbols? I've been repeatedly pasting in cno <- cno + 1 cno plot(x,y,pch=cno) which is tedious, and besides I have to note down the symbols by hand. Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax +64-7 838 4155 Phone +64-7 838 4773 home phone +64-7 856 6705 Mobile +64-21 139 5862 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

...appear to be in the proper directories.] I've done everything except modify the rprofile or the options(). Is that necessary, and if so how exactly do I do that? And, if I do modify them how do I restore them back to their default settings? Any help would be much appreciated. Thanks, Jeff Jorgensen jcjorgensen at wisc.edu