[R] A printing macro

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

[R] stepAIC for overdispersed Poisson

of observations so this should not be problematical) and then use stepAIC with scale = phi. Should this be OK? Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED

Re: [R] A printing macro

(G,dvar), list(dvar=yn))) } There are other ways, such as forming the cross-classification table, setting its dimnames and passing that to chisq.test. On Mon, 13 Nov 2006, Murray Jorgensen wrote: I am exploring the result of clustering a large multivariate data set

[R] Recoding categorical variables

]] - i rbind(y,z) (the recoded y is stored in z) Murray Jorgensen -- 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 Phone +64 7 838

Re: [R] Recoding categorical variables

Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Murray Jorgensen Sent: Wednesday, October 25, 2006 7:13 PM To: r-help

Re: [R] Pooled Covariance Matrix

Thank you, Professor Ripley. Murray Jorgensen Prof Brian Ripley wrote: On Wed, 20 Sep 2006, Murray Jorgensen wrote: I am in a discriminant analysis situation with a frame containing several variables and a grouping factor, if you like: set.seed(200906) exampledf - as.data.frame(matrix

[R] Pooled Covariance Matrix

seemed to imply that doing this was a bit subtle even for ncol=1. Hence I will risk a question to the list. Thanks for any help, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email

[R] unquoting

of the variable's distribution. OK, OK, I could just go AF[,2] [1] 3 4 but once I thought of unquoting I have some sort of inner need to know how to do it! Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton

Re: [R] unquoting

.) get() is just such a function. Murray Jorgensen Prof Brian Ripley wrote: ?get I really think this has nothing to do with `quoting', rather to do with evaluating variables from their names. At first I though you were looking for noquote(), which does unquote in the conventional sense

Re: [R] Fitting models in a loop

Thanks to all for their help. I am busy today but tomorrow I will have time to digest all the feedback and follow up if necessary Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand

Re: [R] Fitting models in a loop

just getting across the elementary ideas. Murray Jorgensen -- 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 Phone +64 7 838 4773

[R] Fitting models in a loop

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: [EMAIL PROTECTED]Fax 7 838 4155 Phone +64 7 838 4773 wkHome +64 7

[R] [Fwd: Re: Parameterization puzzle]

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 [EMAIL PROTECTED] To: Prof Brian Ripley [EMAIL PROTECTED

Re: [R] Parameterization puzzle

*** age1:SmokeYes -1.3755 0.4340 -3.169 0.00153 ** which is just what I wanted. Cheers, Murray Jorgensen Prof Brian D Ripley wrote: On Fri, 21 Jul 2006, Berwin A Turlach wrote: G'day all, BDR == Prof Brian Ripley [EMAIL PROTECTED] writes: BDR R does not know that poly(age,2

[R] Parameterization puzzle

? Can I re-express mod2.glm so that this term does not appear? 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: [EMAIL PROTECTED]Fax 7

[R] princomp and eigen

0.02261302 0.1478618 [4,] 0.1221215 -0.0986234 0.81407655 -0.5591414 I would have expected the princomp component standard deviations to be the square roots of the eigen() $values and they clearly are not. Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff

[R] Extracting Variance components

$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: [EMAIL PROTECTED]Fax 7

Re: [R] Frequencies into Poisson responses

Actually I had just meant that because I had some factors I had a conceptual table. Using ftable() was just my way of getting the factors into data frame form. But thank you for showing me that as.data.frame() does exactly what I want. Murray Jorgensen Prof Brian Ripley wrote: In fact you

[R] Frequencies into Poisson responses

(ftabc) A1 - gl(na,nb,n) B1 - gl(nb,1,n) C1 - gl(nc,na*nb,n) required - data.frame(A1,B1,C1,freqs) required 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: [EMAIL

Re: [R] Pasting data into scan()

:\\Files\\Teaching\\STAT321\\tensile.dat) Read 30 items strength - scan(C:\\Files\\Teaching\\STAT321\\tensile.dat,sep=\t) Read 30 items I will send the file to BR (but not to the list!) Murray Jorgensen Prof Brian Ripley wrote: No tabs are being echoed, so your cut/paste is removing the tabs. What

[R] Pasting data into scan()

messages. Help! 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 Phone +64 7 838 4773 wkHome +64 7 825 0441

[R] Pasting data into scan() - oops!

I forgot to mention that I am using Windows XP. Original Message Subject: Pasting data into scan() Date: Tue, 02 May 2006 11:55:03 +1200 From: Murray Jorgensen [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch The file TENSILE.DAT from the Hand et al Handbook of Small Data Sets

Re: [R] vector-factor operation

Thanks to Berwin Turlach and Petr Pikal for tapply(vec, Fac, mean)[Fac] and Gabor Grothendieck Thomas Lumley for ave(vec,Fac) . Looking at the code for tapply and ave I guess that the latter is to be preferred. Murray Jorgensen Gabor Grothendieck wrote: Look at ?ave ave(vec, Fac) ave

[R] vector-factor operation

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: [EMAIL PROTECTED]Fax 7 838 4155 Phone +64 7 838 4773 wk

[R] Nested error structure in nonlinear model

comments welcome. I regret that I cannot go into much detail about the actual problem. -- 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

[R] nonlinear model: pseudo-design matrix

at the point. (I have figured out how to extract the gradient information from an nls fitted model using the nlsModel part, but I wish to implement a score test, so I need to be able to extract the information at points other than the mle.) Thanks, Murray Jorgensen -- Dr Murray Jorgensen http

Re: [R] Huber location estimate

Prof Brian Ripley wrote: On Thu, 22 Dec 2005, Murray Jorgensen wrote: We have a choice when calculating the Huber location estimate: set.seed(221205) y - 7 + 3*rt(30,1) That's Cauchy, BTW, a very extreme case. Sure, the sort of situation where one might want a robust estimator

Re: [R] Huber location estimate

D'oh! Apologies for wasting everybody's time! Murray Martin Maechler wrote: Murray == Murray Jorgensen [EMAIL PROTECTED] on Thu, 22 Dec 2005 22:13:45 +1300 writes: Murray Prof Brian Ripley wrote: On Thu, 22 Dec 2005, Murray Jorgensen wrote: We have a choice when

[R] Huber location estimate

and Methods}, volume = {A6(9)}, number = {}, pages = {813-827}, year= {1977} } as saying that the two methods were equivalent. Obviously they aren't quite. Comments welcome. Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department

[R] Simple introduction to Bayesian Statistics

too sophisticated starting from chapter 2. Try Introduction to Bayesian Statistics William M. Bolstad ISBN: 0-471-27020-2 by my colleague Bill Bolstad. This is written for a course targeting bright first or second year students and assumes very little background. Cheers, Murray Jorgensen -- Dr

[R] Extracting variance components in lme

? 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: [EMAIL PROTECTED]Fax 7 838 4155 Phone +64 7 838 4773 wkHome +64 7 825 0441Mobile

[R] Extracting Variance Components fro lme

It seems that what I need to get the within group component is as.numeric(VarCorr(fm1Rail.lme)[2,2]) thanks to Bert Gunter and Peter Alspach. Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton

Re: [R] Extracting variance components in lme

Woops! I should have written: as.numeric(VarCorr(fm1Rail.lme)[1,2]) for the within component. -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED

[R] Downloading zip files

I do find it tedious having to right-click Save link as on every individual file from CRAN. I'm sure someone knows a faster way to do it. Cheers, Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato

[R] Error in integrate

(A[j],\n) } -- 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 Phone +64 7 838 4773 wkHome +64 7 825 0441Mobile 021 1395 862

[R] R for teaching multivariate statistics (Summary)

://www.statslab.cam.ac.uk/~pat/AppMultNotes.ps.gz Many thanks to these respondants for their help. Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED

[R] R for teaching multivariate statistics

used or recommended. If you reply directly to me I will summarize to reduce list traffic. 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: [EMAIL PROTECTED

[R] Coarsening Factors

(sites,regions) region [1] III I I II I I I III I III I II III III II I III I III [20] II III I III II III III III II I I Levels: I II III However this seems like using a sledgehammer to crack a nut. Can someone suggest a simpler way to do this task? Murray Jorgensen

Re: [R] Sorting Text Frames

Jumal tervüt teilee Murray Jorgensen -- 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 Phone +64 7 838 4773 wk Home +64 7 825

[R] Sorting Text Frames

: 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

[R] nlme SASmixed in 2.0.1

PROTECTED], REML = REML) : Unable to find the argument object in selecting a method for function deviance Everything works fine under 1.8.1 and plot(PBIB) is of trellis style, which is what I think the authors intend. Cheers, Murray Jorgensen -- Dr Murray Jorgensen http

Re: [R] Order of messgaes/ missing messages

loss on its path through the network. It's statistical, you see! Murray Jorgensen hadley wickham wrote: Note that if this bothers you then you could try a web-based email systems (e.g. yahoo, hotmail, etc.) to receive your r-help messages since the path to any of them would be independent of your

Re: [R] Association between discrete and continuous variable

I'm wondering if mutual information al la Cover Thomas (1991, Ch 2) is not the killer association measure for all types of random variables? Murray Jorgensen PS Yes, this is probably OT! Richard A. O'Keefe wrote: What's the reommended way, in R, to determine the strength of association between

[R] Nested source()s

is fair enough when you think about it. Still it's something that you might want to do. Does anyone know how to get achieve the substance of what nested source() commands would give you? Murray Jorgensen -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics

Re: [R] text editor for R

I tried R-WinEdt a few years ago, but as I remember it interfered with my usual use of WinEdt which is as a front end to MiKTeX. Is there a way to use WinEdt both ways? Murray Jorgensen Marc Schwartz wrote: On Wed, 2004-07-07 at 17:47, Yi-Xiong Sean Zhou wrote: Hi, What is the best text

Re: [R] text editor for R

Oh, yes. I think I did do something like this, but for some reason the two incarnations bothered me. Murray Jorgensen David Scott wrote: On Thu, 8 Jul 2004, Murray Jorgensen wrote: I tried R-WinEdt a few years ago, but as I remember it interfered with my usual use of WinEdt which is as a front

Re: Summary [R] R 1.9.0, special characters in variable names.

This is not what I would call a summary. A summary should: 1. State the original question. 2. Give a pre'cis of the responses. Murray Jorgensen Sixten Borg wrote: Summary: The locale setting in the operating system seems to be involved in what confused me a little bit. Thank you all

Re: [R] MCLUST Covariance Parameterization.

-fit. I think that this sort of approach is likely to be more effective than fitting the fully unstructured model directly. Murray Jorgensen [EMAIL PROTECTED] wrote: Hello all (especially MCLUS users). I'm trying to make use of the MCLUST package by C. Fraley and A. Raftery. My problem

Re: [R] Problem with logtrans, from library MASS

Thanks for the help. The variables do all come from a frame but with various transformations and manipulations. I prefer not to stick them all into a new frame just to call one function. Thanks again, Murray Jorgensen Prof Brian Ripley wrote: On 14 May 2004, Peter Dalgaard wrote: Murray

[R] Problem with logtrans, from library MASS

(price ~ doors+CC+KW+KG+LENGTH, alpha=seq(-1,0,length=30)) Error in eval(expr, envir, enclos) : Object price not found Comments welcome, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email

Re: [R] Re: Mozilla Misbehavior

It was a proxy settings problem. All fixed now. Sorry to disturb the list. Murray Jorgensen [EMAIL PROTECTED] wrote: Hi Murray: Perhaps the errors are thrown by Mozilla because it encounters an encrypted/secure page and does not know what to do (..:)..)? Would it help if you went to Edit

Re: [R] Modalwert

is using some non-standard Microsoft extensions of html? Or is it because my Mozilla 1.6 browser is wrongly configured? Why are the mail archives on secure pages anyway? Murray Jorgensen Thomas Petzoldt wrote: Sonja Dornieden wrote: Hai - kann mir jemand sagen, wie ich den Modalwert in R berechne

[R] Changing Gui preferences

want to get rid of! 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 Phone +64 7 838 4773 wk+64 7 849 6486 homeMobile

[R] Rgui front-end has encountered a problem and needs to close

but the thing worked OK. But back home tonight I get the same problem! The only difference between home and work is that at work I'm connected to a LAN. Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton

[R] Scanning tab-separated numbers

XP Release 1.8.1 I can't seem to scan in tab-separated numbers even when I try to tell R to expect that. (this may be related to the Sweave problem I mentioned a few days ago.) Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University

[R] Sweave and sep = \t

=) f \end{Sinput} \begin{Soutput} [1] 8 8 8 8 8 \end{Soutput} \end{Schunk} But I want the sep = \t to be left as is. 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

Re: [R] references on cluster analysis

I don't really believe that there is any satisfactory definition of the true number of clusters let along a procedure that would reliably find it. Murray Jorgensen Martin Maechler wrote: Back from my vacation, I haven't seen an R-help answer on this (Christian, where have you been

[R] In praise of -

of the scalar fred. Just put - fred at the start of the code block, type any number, then paste the code into the console. 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

Re: [R] In praise of -

you need to do fred - 3 paste fred - 8 paste etc Not a big saving, but some of us are lazy! -- 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

Re: [R] Sweave problem

usage, intended for use only at a keyboard. On Thu, 5 Feb 2004, Murray Jorgensen wrote: 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

[R] Sweave problem

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

Re: [R] Finding Sweave.sty and other problems

with file paths. On Thu, 29 Jan 2004, Murray Jorgensen wrote: Hi, I've just tried to run example-3 from Friedrich Leish. I'm using R 1.8.1 and MiKTeX 2.2 on Windows XP. I go === library(tools) Sweave(example-3.Snw) Writing to file example-3.tex Processing code chunks ... 1

[R] Finding Sweave.sty and other problems

for 1000 14 hyphenation exceptions out of 607 23i,1n,17p,117b,40s stack positions out of 1500i,500n,5000p,20b,32768s No pages of output. === Any comments welcome! Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato

[R] Assignments in loops

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

Re: [R] Assignments in loops

Thanks to Andy, Peter and Roger for drawing my attention to assign(), which is just what I needed and works fine. Murray At 14:11 30/12/2003 +1300, Murray Jorgensen wrote: Greetings all. Any help with the following would be appreciated. I want to create a data frame for each file in a directory

Re: [R] missing data and completed missing data

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 Phone +64 7 838 4773 wk+64 7 849 6486 homeMobile 021 1395 862

Re: [R] mclust - clustering by spatial patterns

- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED

[R] reverse lexicographic order

Hi all, I have some email addresses that I would like to sort in reverse lexicographic order so that addresses from the same domain will be grouped together. How might that be done? Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics

Re: [R] reverse lexicographic order

,''), function(x) paste(rev(x),collapse=''))),...)] } revsort(as.character(1:20)) [1] 10 20 1 11 2 12 3 13 4 14 5 15 6 16 7 [16] 17 8 18 9 19 The ... args are given to order(), so na.last=FALSE and decreasing=TRUE are possibilities. Duncan Murdoch -- Dr Murray Jorgensen http

Re: [R] mode

the mode. At 16:39 12/12/2003 +1300, Murray Jorgensen wrote: Opps! This is what I should have written: The mode of a data vector x might be defined as the limit of m_p as p tends to zero from above and where m_p is the m minimizing sum(abs(x - m)^p). I would not expect the mode so defined

Re: [R] mode

anyone know how to do it? Thanks I don't think the mode of a sample from a continuous random variable is well defined. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Dr Murray Jorgensen http

Re: [R] mode

https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- 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 Phone +64 7 838 4773 wk

[R] Character graphics

are enabled. Professional Graphics are disabled. Use the GPRO command to enable Professional Graphics. Has anyone tried to get similar unprofessional displays out of R? Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics

Re: [R] add a point to regression line and cook's distance

Not a good idea, unless the regression function is *known* to be linear. More likely it is only approximately linear over small ranges. Murray Jorgensen Wiener, Matthew wrote: If you know that the line should pass through (0,0), would it make sense to do a regression without an intercept? You

Re: [R] sampling without repetition

PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help 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 Phone +64 7

[R] Plotting lm() attributes

the removed elements would often have nonmissing values. 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 Phone +64 7 838 4773 wk

[R] ETA for 1.8.1 ?

When is version 1.8.1 likely to be released? 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 Phone +64 7 838 4773 wk

[R] detecting singular matrices

,] -0.6159646 -0.08408654 -0.7468756 [3,] -0.9599328 0.59351831 0.3734378 Which seems to be preferable output. 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

[R] R-project [.com?] [.net?]

that it is not unusual for names close to those of popular sites to be used. It is good that they use a different language or there might well be confusion. Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New

[R] Suming logical vectors

FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE [37] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE [49] FALSE FALSE sum(sort(IATmedian)[0:50]==0) [1] 2 sum(sort(IATmedian)==0) [1] 2 sum(IATmedian==0) [1] NA -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz

Re: [R] Suming logical vectors

Blush! Not too hard to explain at all! Deepayan Sarkar wrote: Your IATmedian has some NAs (which are removed by sort) ? On Thu, 18 Sep 2003, Murray Jorgensen wrote: Can anyone explain the following? [R 1.6.0 Windows XP, yes I will upgrade soon.] Murray sort(IATmedian)[0:50]==0 [1] TRUE

[R] Quit asking me if I want to save the workspace!

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: [EMAIL PROTECTED

Re: [R] Quit asking me if I want to save the workspace!

Rafael A. Irizarry wrote: you can type this: q(no) see the help file for q Still more work than two mouse clicks. -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: [EMAIL PROTECTED

Re: [R] Quit asking me if I want to save the workspace!

Ah! now that tells me what I want to know. I was trying to type C:\Program Files\R\rw1071\bin\Rgui.exe --no-save instead of C:\Program Files\R\rw1071\bin\Rgui.exe --no-save into the Target box. Silly me! Jason Turner wrote: On Wed, 2003-09-17 at 14:26, Murray Jorgensen wrote: Rafael

Re: [R] R tools for large files

to see if I can build up my random subsets of files this way. BTW, I now estimate the file at about 100,000 lines so indeed, it is not all that large! Murray Jorgensen Prof Brian Ripley wrote: On Mon, 25 Aug 2003, Murray Jorgensen wrote: At 08:12 25/08/2003 +0100, Prof Brian Ripley wrote: I

Re: [R] R tools for large files

not be put down to laziness or rudeness! Cheers, Murray Jorgensen PS Give my regards to Chris Hennig. Martin Maechler wrote: Hi Murray, from reading your summarizing reply, I wonder if you missed the most important point about connections connection := generalization of file): Once you open() one

[R] R tools for large files

in a file - form a data frame by selecting all cases whose line numbers are in a supplied vector (which could be used to extract random subfiles of particular sizes) Does anyone know of a package that might be useful for this? Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz

Re: [R] R tools for large files

Could you be more specific? Do you mean the chapter on connections? Ko-Kang Kevin Wang wrote: Hi, Have you looked at R Data Import/Export? On Mon, 25 Aug 2003, Murray Jorgensen wrote: __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch

Re: [R] R tools for large files

, Andrew C. Ward CAPE Centre Department of Chemical Engineering The University of Queensland Brisbane Qld 4072 Australia [EMAIL PROTECTED] Quoting Murray Jorgensen [EMAIL PROTECTED]: I'm wondering if anyone has written some functions or code for handling very large files in R. I am working

Re: [R] R tools for large files

responses interspesed below. BTW, I forgot to mention that I'm using Windows and so do not have nice unix tools readily available. On Mon, 25 Aug 2003, Murray Jorgensen wrote: I'm wondering if anyone has written some functions or code for handling very large files in R. I am working with a data

[R] Interlacing two vectors

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

Re: [R] r-help using random generating

/mailman/listinfo/r-help -- 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 Phone +64 7 838 4773 wk+64 7 849 6486 homeMobile

[R] help with likelihood contour plot

) [1] TRUE Thanks, 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 Phone +64 7 838 4773 wk+64 7 849 6486 homeMobile

[R] lda on curves

, ...) { if (is.null(class(x))) class(x) - data.class(x) UseMethod(lda, x, ...) } which isn't very helpful. Any ideas about how to perform an unscaled canonical variates analysis? Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics

Re: [R] google

- toysrus.com is second. __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand

Re: [R] order() on vector with Inf's

[EMAIL PROTECTED] wrote: On Thu, 23 Jan 2003, Murray Jorgensen wrote: [useful suggestions omitted] llbet[llbet==Inf] - -Inf should correct the problem, but why not correct the calculation? Not a bad idea, but the numerical difficulties always occur well away from the likelihood maximum

[R] order() on vector with Inf's

contains some NAs which I wish to ignore. Unfortunately my llbet contains Inf's, which are not points of infinite likelihood, if anything they should be -Inf's. Anyway na.last does not seem to help me with these. What should I do? Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz