Re: [R] off-topic: better OS for statistical computing

. Linux. It's best for ***everything***. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped

Re: [R] computing distance in miles or km between 2 street addre

On 7/09/2007, at 10:17 AM, (Ted Harding) wrote: On 06-Sep-07 18:42:32, Philip James Smith wrote: Hi R-ers: I need to compute the distance between 2 street addresses in either km or miles. I do not care if the distance is a shortest driving route or if it is as the crow flies. Does

Re: [R] Pie Chart Labels

On 5/09/2007, at 9:37 AM, Greg Snow wrote: The best option is to use a bar chart or dot chart instead of a pie chart. Right on, Red Freak!!! :-) cheers, Rolf Turner

Re: [R] Help on inverse distribution

On 5/09/2007, at 11:16 AM, Lisa Hu wrote: To make it specific, I need to simulation Y with inverse beta distribution, that is, Y~inverseF(X), where F is the CDF of beta distribution. THANKS On 9/4/07, Lisa Hu [EMAIL PROTECTED] wrote: Dear All, I need to use the inverse of some

Re: [R] Excel

? cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing list

Re: [R] Excel

cogently expressed! cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help

Re: [R] a new-bie question about obtaining certain value from the print out

and indicates a certain lack of insight which *could* do harm in more complex situations. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti

Re: [R] combining datasets by row

() for(i in 1:count) rslt[[i]] - read.csv(paste(etc.)) do.call(rbind,rslt) } cheers, Rolf Turner ## Attention:\ This e-mail message

[R] Linear modelling confusion.

, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] piecewise linear approximation

is a vector of candidate values for the ``join point'', e.g. b - seq(lo,hi,length=100), where lo and hi are suitably chosen lower and upper bounds on the location of the join point. cheers, Rolf Turner bsr - function(x,y

Re: [R] Excel

on Excel, available at the URL: http://www.stat.uiowa.edu/~jcryer/JSMTalk2001.pdf Executive summary: Friends don't let friends use Excel for statistics. cheers, Rolf Turner

Re: [R] Choice of data.frame column by index?

On 24/08/2007, at 12:51 PM, Moshe Olshansky wrote: This won't work since it produces a matrix (try this). On the contrary, Patrick's solution is correct. I tried it. It works just fine. cheers, Rolf Turner What should work

Re: [R] small issue with densityplot

On 22/08/2007, at 8:44 AM, Douglas Bates wrote: In some ways it might be better to remove the dependence on coda and write the classes and methods in the lme4 package. That way I can use S4 classes and irritate all those people who rail against S4 classes and methods (and you know who you

Re: [R] (Most efficient) way to make random sequences of random sequences

On 22/08/2007, at 8:49 AM, Emmanuel Levy wrote: Hi, I was wondering the what would be the (most efficient) way to generate a sequence of sequences, i mean: if I have 1,2 and 3. I'd like to generate a sequence of length N*3 (N ~ 1,000,000 or more) Where random permutations of the

Re: [R] tackle memory insufficiency for large dataset using save() load()?

, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] tackle memory insufficiency for large dataset using save() load()?

! cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing

Re: [R] Does anyone.... worth a warning?!? No warning at all

'' --- as Mathew *initially* (and quite reasonably) expected it to do. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti

Re: [R] time series periodic data

ed., Wiley, 2000) for a good introduction to complex demodulation. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped

[R] Covariance of data with missing values.

that anyone may be able to give. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped

Re: [R] Covariance of data with missing values.

on unobserved). Read the documentation for the various function *very* carefully! Drop me a line if you want more info. Best wishes, Ted. Also, consider mlest() in the mvnmle package. On 15-Aug-07 21:16:32, Rolf Turner wrote: I have a data matrix X (n x k, say) each row of which

Re: [R] CDF of a gamma distribution

On 15/08/2007, at 9:53 AM, Jake Verschuyl wrote: I have some bird flight height data that follows a gamma distribution. The data (x) goes from 0 to 700 meters (n=1055). The calculated parameters calculated from the fitdistr(x) are (shape = 5.1379, rate = 0.017541), and therefore

Re: [R] deldir package - voronoi

cleverer than I will come to your rescue. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped

Re: [R] Colours in R, was: deldir package - voronoi

On 13/08/2007, at 11:27 AM, Wolfgang Huber wrote: Dear Rolf and Binabina perhaps this is of use to some: Colour for Presentation Graphics. Ross Ihaka. www.stat.auckland.ac.nz/~ihaka/colour/color.pdf Choosing Color Palettes for Statistical Graphics Achim Zeileis and Kurt Hornik.

Re: [R] small sample techniques

to Nair's original question; it is clear that he is interested in a two-independent-sample t-test. cheers, Rolf Turner ## Attention:\ This e-mail

Re: [R] manipulating arrays

?append cheers, Rolf Turner On 28/07/2007, at 2:12 AM, Henrique Dallazuanna wrote: Hi, I don't know if is the more elegant way, but: X-c(1,2,3,4,5) X - c(X[1], 0, X[2:5]) -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 27/07/07

Re: [R] spatstat - Fitting a Strauss model with trend determined by kernel density smoother

family form --- and hence is not estimated by ppm(), at least not directly. cheers, Rolf Turner ## Attention:\ This e-mail message is privileged

Re: [R] for loop doesn't stop with upper loop value

, Rolf Turner ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing list https

Re: [R] generating correlated Bernoulli random variables

and p2) in order that the values of q1 and q2 lie between 0 and 1. I.e. for some values of rho your goal will be impossible to achieve. cheers, Rolf Turner

Re: [R] generating correlated Bernoulli random variables

Typo in my previous email: sqrt((1-p1)*(1-p2)) should have read sqrt(p1*(1-p1)*p2*(1-p2)) I think! cheers, Rolf ## Attention:\ This e-mail message is privileged

Re: [R] Extracting sums for individual factors in data frames

?tapply ## Attention:\ This e-mail message is privileged and confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] Catenating factors.

I was recently asked by one of new colleagues how to combine, or catentate, or concatentate two factors. If x and y are factors, doing c(x,y) appears to coerce x and y to numeric mode before catenating them. So what does one do if one wants the result to be a factor whose levels are the union

[R] quartz() on MAC OSX

much hand-cuffed by the officious ITS policies here as to what I can install on my Mac. (Effectively, nothing.) ``cheers'' (Hah!), Rolf Turner [EMAIL PROTECTED

Re: [R] quartz() on MAC OSX

On 22/5/07 6:48 AM, hadley wickham [EMAIL PROTECTED] wrote: Hi Rolf, Two possible solutions: * DISPLAY=0.0.0.0:0 R - and then X11() should work without having to use xterm * install.packages(CarbonEL); library(CarbonEL); quartz() It is clear that life is determined to frustrate

Re: [R] Re : Generate a random bistochastic matrix

is that it would be a hard problem. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] generate random numbers that sum up to 1

As I have previously asked, in response to a similar question: Is this a homework problem? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch

Re: [R] Block comments in R?

I must say I agree with Richard O'Keefe who wrote: I wrote: R documentation comments really belong in [.Rd] files where help() can find them. Barry Rowlingson [EMAIL PROTECTED] replied: R documentation comments belong in .Rd files at the moment, but how joyous would

Re: [R] ifelse(logical, function1, function2) does not work

method that I wrote got dispatched, but the function method didn't? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help

Re: [R] setting the environment variable R_LIBS

/library export R_LIBS HTH. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] How can I generate this numbers

This sounds so simple it must be a homework problem, no? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list

Re: [R] How can I generate this numbers

Ricardo Rios wrote: Hi Rolf Turner, I have a statistical model, it model need this numbers for calculate the probability. This numbers must be random. For example I need that magicfunction(3) [1] 0.3152460 0.5231614 0.1615926 magicfunction(3) [1] 0.6147933 0.3122999 0.0729068

Re: [R] Print/Save/Cat/Write list

?sink __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] unvector ?

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] Matrix multiplication using apply() or lappy() ?

, Rolf Turner __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

[R] Local library under Windoze.

, what's the point of having this property for folders?) Can anyone enlighten me? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch

[R] Summary --- Local library under Windoze.

by the *relative* path name, with no problem. Thanks to all. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list

Re: [R] write a table to file with unequal length of lists

) cheers, Rolf Turner [EMAIL PROTECTED] Original message: Dear R helpers, To illustrate my problem, here is a simplified example. I want to write a table to a file similar to: x a 1

[R] Summary and thanks: .Rprofile under Windoze.

. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] .Rprofile under Windoze.

! '' as it should. I would appreciate enlightenment. Ta. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Adding a footnote in plot-window in R

?legend __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] fitting a gaussian to some x,y data

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] Creating a pixel image

, apparently even to the cognoscenti. cheers, Rolf Turner [EMAIL PROTECTED] P. S. A library, as someone --- Martin Maechler??? --- has repeatedly tried to explain to the list

Re: [R] model = F causing error in polr()

Try ``model = FALSE'' rather than ``model = F'' and see if it makes a difference. You make have an unwanted variable named ``F'' lurking somewhere. (In general it is a *bad* idea to use ``F'' when ``FALSE'' is intended.) cheers, Rolf

Re: [R] OT UNIX grep question

interpretation of the dollar sign you have to use *single* quotes. grep '^dog$' /usr/share/dict/words *does* work. (Try it!) cheers, Rolf Turner [EMAIL PROTECTED

Re: [R] minimization a quadratic form with some coef fixed and some constrained

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] Comparison of linear models

be a reasonable place to start. Bon chance. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https

Re: [R] information about a function

the calculations for you; all you would be doing is reinventing the wheel. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Invoke operating system command

?system __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] modeling logit(y/n) using lrm

are called. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] if syntax error :(

readable and avoids unintended consequences. Another infelicity in your code: ``plot=F''. Use ``plot=FALSE''. (Note that the symbols ``F'' and ``T'' are assignable, *unlike* ``TRUE'' and ``FALSE''.) cheers, Rolf Turner

Re: [R] Density Estimation

(spiffy,0,3) 0.2353400 with absolute error 2e-09 cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] lm() variance covariance matrix of coefficients.

summary(object)$cov.unscaled __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] lm() variance covariance matrix of coefficients.

Peter Dalgaard wrote: Rolf Turner [EMAIL PROTECTED] writes: summary(object)$cov.unscaled You need to multiply that with sigma. However, vcov(object) is easier. Well, I thought unscaled meant unscaled --- the plain unvarnished covariance matrix! I figure that multiplying

Re: [R] Faster way to zero-pad a data frame...?

'' is to allow for 0-origin. data2 - ts(data2,start=0) ??? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing

Re: [R] Automate concatenation?

Well, I like do.call(): puddy.tat - do.call(c,lapply(paste(a,1:10,sep=),get)) cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list

Re: [R] newbie question: ROW average

Peter Dalgaard wrote: How about dim(M) - c(4,248,74) mn - apply(M, c(2,3), mean) Hey! That's sexy! Much better than my kludgy suggestion! cheers, Rolf __

Re: [R] newbie question: ROW average

be done is via m.new - t(apply(array(t(m.old),dim=c(74,4,992/4)),c(1,3),mean)) cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list

Re: [R] normality testing with nortest

, then the p-value is uniformly distributed on [0,1]. This should be kept in mind when assessing the ``instability'' of p-values. cheers, Rolf Turner

Re: [R] normality testing with nortest

. cheers, Rolf Turner [EMAIL PROTECTED] Original message: I have a question regarding normality testing with the nortest package. I have to admit, my question is so general

Re: [R] iraq statistics - OT

plausible. But still not very plausible. Or maybe major American cities are even worse than we had been led to believe. cheers, Rolf Turner __ R-help@stat.math.ethz.ch

Re: [R] How can you buy R?

Jarek wrote: At one company I was working for, I had to run all the licenses of all the software I had on my machine, through the legal department. When they read GNU Public License (GPL) their only comment was: We have no idea what that license means. Do not touch any software using it.

[R] Off topic --- help in locating a source.

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R

Re: [R] Can't there be a cd command?

Duncan Murdoch wrote (amongst other things): Statistical computing is not easy, so how could R be? Who has ever claimed it is? Any package that makes statistical computing appear to be easy is probably giving you wrong answers half the time, or is extremely limited in scope.

[R] Symbol __f95_sign not found.

for any hints. cheers, Rolf Turner [EMAIL PROTECTED] P. S. Version info: version _ platform sparc-sun-solaris2.9 arch sparc os

Re: [R] question in fitting AR model in time series?

See argument ``fixed'' in help(arima). cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] 3D pie

charts are terrible. Don't mess with them. Period. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] I am surprised (and a little irritated)

a helpful, humble, ***practical*** response, like that from Detlef Steuer, from the Microsoft(dung) Monolith. cheers, Rolf Turner [EMAIL PROTECTED

Re: [R] About poisson distribution fitting and testing

will also (almost surely!) tell you that the data are not Poisson. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

[R] Creating an environment for a function.

be missing a point or three. As I said, I don't really grok environments. Given that what I want to do makes any kind of sense at all, can someone start me off in the right direction? cheers, Rolf Turner

Re: [R] Creating an environment for a function.

Brian Ripley wrote: Can I also suggest local()? This does a similar thing in a perhaps more natural way. Sorry, I'm not with you. I'm slow, and as I said, I don't really grok environments. Let's look at a toy example. Suppose I want to create a function foo: function(x){x^n} and

Re: [R] Creating an environment for a function.

Thomas Lumley wrote: foo - local({ n - 4 function(x) {x^n} }) Bewdy! Thanks a lot. Slowly but slowly ( :-) ) I learn things! cheers, Rolf

[R] Intercepts in linear models.

. Is there a sexier way? cheers, Rolf Turner __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R

[R] smooth.spline

, it might be nice to have some sort of warning in the help file to keep young players from getting confused. cheers, Rolf Turner [EMAIL PROTECTED

[R] Derivative of a splinefun function.

? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R

[R] An optim() mystery.

. cheers, Rolf Turner [EMAIL PROTECTED] P. S. The help on optim() says the Nelder-Mead method is ``relatively slow''. However for this problem optim() seems to come in about 10 seconds when

[R] Problem with installing home-made package under Windoze.

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

[R] Puzzled over curve() syntax.

, but I would like to deepen my understanding. There's probably something about the ``expression'' concept that I'm not grokking here Thanks for any insight. cheers, Rolf Turner [EMAIL PROTECTED

[R] Getting eps into Word --- thanks.

, Rolf Turner __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Dots in function names

Peter Dalgaard wrote: Well, come the S4 revolution and dots will cause trouble no more... Yes but an infinite number of other things will cause infinitely more trouble. God save us from S4. cheers, Rolf Turner

Re: [R] User error (was arima.sim bug?)

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Getting eps into Word documents.

into a Word document. If so, would you be willing to share it with me and my student? If so, please be gentle in your explanation. I am not myself (repeat ***NOT***) a user of Word! Thanks. cheers, Rolf Turner

[R] (Off topic) Limit question --- corrected!!!

. Can anyone help me out? cheers, Rolf Turner __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R

[R] (Off topic) Limit question solved.

apologies for confusing everyone with my initial mis-statement of the problem. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https

[R] (Off topic) Limit question.

me out? cheers, Rolf Turner __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

[R] Testing.

Please ignore this message; I apologise for the annoyance. cheers, Rolf __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] monte carlo simulations/lmer

I think you want to use the function try(); see ?try. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Dating Objects

to send a separate email about this. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

[R] Bug in order with na.last=FALSE?

. It seems to be the complicated strings which form the entries of x that are messing things up; y - c(b,c,NA,a,NA) y[order(y,na.last=FALSE)] gives [1] NA NA a b c is it ought. cheers, Rolf Turner

Re: [R] interpolation function

, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Adding sum to derivatives table

deriv(expression(sum(x)), x) does not make any sense. cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing

Re: [R] Adding sum to derivatives table

Martin C. Martin wrote: Rolf Turner wrote: deriv(expression(sum(x)), x) does not make any sense. Good point. But this does: deriv(expression(sum(log(a*x))), a) where a is a scalar. Okay --- I see what you're getting at now. But I think

  1   2   3   >