Re: [R] coarsening levels problem

g # [1] other other akl other akl akl akl other other akl other other #[13] akl akl other akl other akl akl akl akl akl other other #[25] other akl akl akl akl other other other #Levels: akl other A.K. - Original Message - From: Murray Jorgensen To: r-help

[R] coarsening levels problem

Greetings from New Zealand. I hope that my difficulties are clear from the following output. Ideas gratefully received! Murray Jorgensen cc <- scan() 3000 3500 2500 2500 3000 3000 3000 3000 3000 2000 2000 2000 2000 4800 3000 4500 3000 2000 2000 4400 3000 3000 2500 3000 3000 2000 3000 2

[R] Converting Fortran or C++ etc to R

doing all this is so that I may try out some variants of the program. I expect this will be much easier to do in R! Thanks to all who replied. Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand

Re: [R] Converting Fortran or C++ etc to R

the Fortran (C/C++) code at every stage of the rewrite. Don't just write it all in one go and then hope it works! Small steps Barry -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zea

[R] Converting Fortran or C++ etc to R

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 Departme

Re: [R] Sweaving quotes

010, at 7:43 PM, Murray Jorgensen wrote: The significance code line to summary() applied to an lm() fitted model object is Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 The corresponding line in the LaTeX source produced by Sweave is Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘

Re: [R] Sweaving quotes

quotes -- since this is an issue for the R manuals, see our documentation for ways around this (such as package ae). On Wed, 28 Jul 2010, Marc Schwartz wrote: On Jul 28, 2010, at 7:43 PM, Murray Jorgensen wrote: The significance code line to summary() applied to an lm() fitted model obje

Re: [R] Sweaving quotes

Thanks, Marc, and also Jay Kearns. After experimentation I think useFancyQuotes = FALSE my be best as TeX style can look a bit funny for R output. Cheers, Murray Marc Schwartz wrote: On Jul 28, 2010, at 7:43 PM, Murray Jorgensen wrote: The significance code line to summary() applied to

[R] Sweaving quotes

Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: m...@waikato.ac.nzmajorgen...@ihug.co.nz Fax 7 838 4155 Phone +64 7 838 4773 wkHome +64 7 825 0441 Mobile 021 0200 8350

Re: [R] Sweave and scan()

. Thank for that. BTW my previous "Thanks for all your help" was meant to be "Thanks to all for your help". Cheers, Murray Duncan Murdoch wrote: On 28/07/2010 6:33 PM, Murray Jorgensen wrote: Omigod! The archival links shows that this was the same problem that caused me to give

Re: [R] Sweave and scan()

- c(64, 62, 66, 65, 62, 69, 72, 72, 70) part = c("Soprano","Soprano", "Soprano", "Alto","Alto","Tenor", "Tenor", "Bass","Bass") @ which does give the output I want (though in a different style) and

Re: [R] Sweave and scan()

Both suggestions generate similar errors to those of the original code. I would also be worried if the results would not puzzle my students. But thanks! Murray David Winsemius wrote: On Jul 26, 2010, at 11:54 PM, Murray Jorgensen wrote: I am introducing the scan() function to my class

[R] Sweave and scan()

ocessing 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 Mu

Re: [R] Html help

:13 a.m., Duncan Murdoch wrote: Murray Jorgensen wrote: I have just installed R 2.11.1 on my XP laptop. I like html help [...] How do I go about getting a local set of html help files? Since 2.10.0, HTML help is generated on demand. It doesn't go off your local computer, it works lo

[R] Html help

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 Jorg

[R] R or C++ on FreeNX servers

are and if so, how the parallelizing would be handled. Regards, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: m...@waikato.ac.nzFax 7 838

[R] R or C++ on FreeNX servers

o, how the parallelizing would be handled. Regards, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: m...@waikato.ac.nzFax 7 838 4155 Phone +64

[R] Help with package reshape

2604218 NA 11 NA 2691185 NA 12 NA 2753555 NA 11 NA 308 NA 948 NA 12 34928082511 12 11 Suggestions gratefully received! Regards, Murray Jorgense

[R] Help with package reshape

NA 11 NA 2691185 NA 12 NA 2753555 NA 11 NA 308 NA 948 NA 12 34928082511 12 11 Suggestions gratefully received! Regards, Murray Jorgensen -- Dr Murray

Re: [R] two-factor linear models with missing cells

+ Biv" now refers to the mean of the (3,4) cell rather than the that of the (1,4) cell making the (3,4) cell a sort of 'honorary' member of the first row. This could have been done to the (2,4) cell but I guess the rule is to drop the cell with the highest sum of row and column number

[R] two-factor linear models with missing cells

combinations and how does lm() decide what parameters to report? Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: m...@waikato.ac.nzFax 7 838 41

[R] two-factor linear models with missing cells

combinations and how does lm() decide what parameters to report? Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: m...@waikato.ac.nzFax 7 838 41

[R] Tick label orientation

dj=TRUE) axis(1,at=seq(0,24,4)) would give me horizontal tick labels. It doesn't. What would? Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: m...@waikato.ac.nz

Re: [R] sampling from Laplace-Normal

: Theory and Methods, 33(8), 1733-1753. Perhaps you can get around your problem using this representation or alternatively reproduce the problem with one of the component rvs. Cheers, Murray Jorgensen Hi, I have to draw samples from an asymmetric-Laplace-Normal distribution: f(u|y, x, beta

Re: [R] step() and stepAIC()

Prof Brian Ripley wrote: On Sun, 12 Oct 2008, Murray Jorgensen wrote: The birth weight example from ?stepAIC in package MASS runs well as indeed it should. However when I change stepAIC() calls to step() calls I get warning messages that I don't understand, although the output is si

[R] step() and stepAIC()

: using type="numeric" with a factor response will be ignored (and three more the same.) Checked with 2.7.0 on Os X and 2.6.2 on Windows. Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton,

Re: [R] lattice plotting character woes

mbol",plsy) xyplot(y ~ x |Region) # only 1,2,3,4 used as plotting symbol # I expected 1:6 Mark Leeds has pointed out that whatever numbers appear as plotting characters on the screen, they are replaced by circles when you save a pdf via pdf() xyplot() dev.off() I have

Re: [R] lattice plotting character woes

ts,pch) { panel.xyplot(x,y,pch=pch[subscripts]) } )) cheers, Rolf ## Attention: This e-mail message is privileged and confidential. If you are not the intended recipient please delete the

[R] lattice plotting character woes

,2] > zz = as.character(patches$Buffer.zone/5) > plsy <- trellis.par.get("plot.symbol") > plsy$pch = zz > trellis.par.set("plot.symbol",plsy) > xyplot(y ~ x |Region) > table(zz) zz 0 1 2 4 19 20 20 20 I would dearly like to have my plotting symbol

Re: [R] Extracting formula from an lm object

Thanks to Professors Ripley and Turlach for their help. Brian's solution was simple and elegant but Berwin's was also informative. Murray Jorgensen Prof Brian Ripley wrote: On Sun, 24 Aug 2008, Berwin A Turlach wrote: G'day Murray, On Sun, 24 Aug 2008 18:34:39 +1200 M

[R] Extracting formula from an lm object

y code I will fit another linear model of the form z ~ A + B + x for some other response variable z. I would be grateful for any suggestions of a nice way to do this. Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics,

[R] Programming Concepts and Philosophy

rogram rather than using "10" and "120" throughout the script. This sort of stuff is not specifically R but can be a problem for students with little programming background. I am happy to summarise responses. Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.wa

[R] Recommended Packages

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, N

Re: [R] [OT] good reference for mixed models and EM algorithm

2007. useful although it takes a Bayesian rather than EM approach. Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED]Fax 7 838 4155 P

[R] Reading R-help Digests with Mozilla Thunderbird

text in a message. [I'm using Windows XP]. Does anyone know how to turn off this Thunderbird feature? Also what would be the best way to raise this problem with the Thunderbird developers? Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.htm