Re: [R] Bug in lowess

; > > > > R version 2.2.1, 2005-12-20, i486-pc-linux-gnu > > attached base packages: > [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" > [7] "base" > > other atta

Re: [R] Extrapolated regression lines

sults plot(Y ~ X) lines(pred ~ newdat, col = "red") HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC & ENSIS, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk G

Re: [R] C code for KalmnaLike

This is not an offer (or solicitation of an offer) to buy/se...{{dropped}} > > __ > 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 &g

Re: [R] NA-handling in glm.fit?

You'll probably want to run na.omit or na.exlude on a combination of the response and predictor matrix (the arguments you are passing as x and y to glm.fit) to remove incomplete cases. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson

Re: [R] Mantel test!

s are full, square, symmetric matrices, then you can use: mantel.rtest(as.dist(your_mat1), as.dist(your_mat2)) See ?as.dist HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC & ENSIS, UCL Geography

Re: [R] nontabular logistic regression

hess Hall > Department of Biological Sciences > Auburn University > Auburn, AL 36849 > 334-329-9198 > FAX 334-844-9234 > http://www.auburn.edu/~stratja > > ______ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailm

Re: [R] concatenate 2 data.frames

_ > 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] Change Axis Size

"green"), type = "l", lty = "solid") or matplot(x, cbind(y1,y2,y3), type = "n") matlines(x, cbind(y1,y2,y3), col = c("red", "blue", "green"), lty = "solid") HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] Kmeans centers

of what you are trying to do, and how you tried to do it but failed. > and provide commented, minimal, self-contained, reproducible code. ^^^ G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin S

Re: [R] Xemacs, ESS, R config issue

ignment anymore, so the only place I can think of at the mo is that you want to type a "_" in a variable name? Also, this is best addressed on the ESS mailing list: ess-helpATstatDOTmathDOTethzDOTch (just remove the AT and DOTs) G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] Kmeans centers

t to find, why do you need a clustering algorithm to find them for you? G > > thk you > > > > 2007/3/29, Gavin Simpson <[EMAIL PROTECTED]>: > On Thu, 2007-03-29 at 15:02 +0200, Sergio Della Franca wrote: > > Dear R-Helpers, >

Re: [R] substitute NA values

gt; > > 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] Abundance data ordination in R

ailing 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] partial R

t missed out the fundamental step of loading the package from the library: install.packages('car') library(car) ?cr.plots HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC

Re: [R] division of decimal number

vide = get("/") > > > > divide(0.285,2) > [1] 0.1425 > > Is that what you want? > > Barry You can use the function directly, by quoting it: > "/"(0.285, 2) [1] 0.1425 G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] queries

rovide commented, minimal, self-contained, reproducible code. Please do. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC & ENSIS, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Street, Lo

Re: [R] Randomising matrices

0 2532 > Web: http://www.geos.ed.ac.uk/homes/s0455078 > > __ > 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 pro

Re: [R] partitioning variation using the Vegan CCA routine?

> 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] how to be clever with princomp?

ASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geogr

Re: [R] randomForest gives different results for formula call v. x, y methods. Why?

30.06 virginica 0 3470.06 HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC [f] +44 (0)20 7679 0565 UCL Department of Geography Pearson Building [e] gavin.simpso

Re: [R] Too slow to execute!

ds on my desktop - which has similar processor but a lot more RAM: > system.time({ + n.perm <- 10 + Re <- matrix(ncol = 20, nrow = n.perm) + for(i in seq_len(n.perm)) { +x <- rnorm(20) +Re[i,] <- (x-2*10)/20 + } + }) user system elapsed 3.336 0.056 3.394 HTH G -

Re: [R] xyplot with grid?

..) { panel.grid(...) panel.xyplot(x, y, ...) }) Not sure if that is the easiest way, or the best, but that's how I've learnt to use lattice recently. The v and h arguments are passed to panel.grid as part of "..." and just tell it to plot the grids at the tick m

Re: [R] Goodness of fit for hclust?

alysis. The example on that help page shows how to use the function to get the correlation between the original distances and cophenetic distances. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC

Re: [R] basic problem but can't solve it

ly for two solutions. The tapply one seems cleaner and easier to get the vector you need, the aggregate version needs an extra step: aggregate(dat$Measure, by = list(Month = dat$Month), mean)$x ^^ Note the $x at the end to subset the object

Re: [R] R-About PLSR

Before you tried to use the functions you mention? HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATN

Re: [R] 'trim' must be numeric of length one？

___ > 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. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin S

Re: [R] function gdist, dist and vegdist in mvpart

et - like I do above. Readers of the list wont have your "ba12", so extra work is required to investigate the problem. > > Thanks for helping. > Sina You're welcome, hope this helps. Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simps

Re: [R] Linux editor like WinEdt?

I just haven't figured out enough of the arcane finger gestures for Xemacs yet to use it as my LaTeX environment - the one click build tools in Kile are very useful and easily customisable. G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson

Re: [R] what is the function in R for MARS and PRIM?

t suggest that this method currently isn't available in R. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. &

Re: [R] is there a formatted output in R?

> https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide! > > http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html> > > > > > > > > > > -- > > Jim Holtman

Re: [R] GAM using R tutorials?

but a starting point might be the RNews article Simon Wood wrote. That should get yourself started. But if you'd done the search yourself, you wouldn't have had to wait for someone on the list to do it for you. Finally - Please read the posting guide - it is there for a reason. HTH G

Re: [R] GAM using R tutorials?

ion of ?gam, which the posting guide does ask you to read. One can only assume you missed this reference when you looked at the page or did not realise the significance of it. And for the record, wherever possible I do try to reply constructively - if I didn't have an answer to your question I

Re: [R] multiple return values

= ran) return(retval) ## or in one step: ## return(list(mean = m, range = ran)) } foo(100) bar <- foo(100) bar HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow

Re: [R] R-squared for MARS models

], ozone[,]) (a.summ <- summary(lm(ozone[,1] ~ a$x-1))) a.summ$r.squared HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% *Note new Address and Fax and Telephone numbers from 10th April 2006* %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] message posting

u can set this up somewhere amongst the myriad of prefs in Outlook. G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% * Note new Address, Telephone & Fax numbers from 6th April 2006 * %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson

Re: [R] Problem with data frame

%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% * Note new Address, Telephone & Fax numbers from 6th April 2006 * %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson ECRC & ENSIS [t] +44 (0)20 7679 0522 UCL Department of Geography [f] +44 (0)20 7679 0565 Pearson Build

[R] Running R on Windows 2000 Terminal Services

-- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% * Note new Address, Telephone & Fax numbers from 6th April 2006 * %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson ECRC & ENSIS [t] +44 (0)20 7679 0522 UCL Department of Geography [f] +44 (0)20 7679 0565 Pearson

Re: [R] Running R on Windows 2000 Terminal Services

On Tue, 2006-04-25 at 17:32 +0100, Barry Rowlingson wrote: > Gavin Simpson wrote: > > Dear list, > > > > My employer uses a Windows 2000 Terminal Server-based system for its > > college-wide managed computer service - computers connect directly to > > the WTS se

Re: [R] www.r-project.org

op-level menu items as possible, say 5-6 max... OK the R website is not like .com type website, but you wouldn't want to flood users with too many options up front. G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% * Note new Address, Telephone & Fax

Re: [R] par(mfror=c(1,2))

nch, width = 16 * cmToInch, pointsize = 10, onefile = FALSE) opar <- par(mfrow = c(1,2), pty = "s") plot(rnorm(100), rnorm(100), asp = 1) plot(rnorm(100), rnorm(100), asp = 1) par(opar) dev.off() HTH, G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] How to use "identify"

need: plot(our.frame2$c1, our.frame2$c9) # you don't need the extra brackets identify(our.frame2$c1, our.frame2$c9) G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44

Re: [R] Mantel Ranomization test

://cc.oulu.fi/~jarioksa/softhelp/vegan.html Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] gavin.simpsonATNOSPAMucl.a

Re: [R] Creating a correlation Matrix

e amount of data21 variables for 3471 times. I want to see how each of the variables correlate to each other. Any help would be appreciated, including which package and which functions I should use to do this. Thanks, Jessica Higgs Masters Student Department of Meteorology Penn State Uni

Re: [R] strange behavior of loess() & predict()

/stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)

Re: [R] retrieving p-values in lm

[[???]][[]] > > -- > Patrick Kuss > PhD-student > Institute of Botany > University of Basel > Schönbeinstr. 6 > CH-4056 Basel > +41 61 267 2976 > > ______ > R-help@stat.math.ethz.ch mailing list > https://stat.eth

Re: [R] tapply and weighted means

y(dat, dat$x, function(z) weighted.mean(z$y, z$w))) but if you want to easily access the numbers you need to do a little work, e.g. as.vector(res) Also, I don't see a function wtd.mean in standard R and weighted.mean() doesn't have a weights argument, so I guess you are using a function f

Re: [R] In which application areas is R used?

sage- > >> From: [EMAIL PROTECTED] > >> [mailto:[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/p

Re: [R] Timeliness of precompiled binaries--R-2.2.1 still not available as RPM

as in all the EL3/4 and FC3/4 i386 and x86_64 directories I could be bothered to look at on the UK Bristol mirror for example. G > Regards, > Scott Waichler > Pacific Northwest National Laboratory > scott.waichler _at_ pnl.gov -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] Timeliness of precompiled binaries--R-2.2.1 still not available as RPM

On Mon, 2006-01-30 at 18:53 +0100, Peter Dalgaard wrote: > Gavin Simpson <[EMAIL PROTECTED]> writes: > > ...Mores the point, though, which CRAN mirror were you looking at? R > > 2.2.1 was in all the EL3/4 and FC3/4 i386 and x86_64 directories I could > > be bothered t

Re: [R] Timeliness of precompiled binaries--R-2.2.1 still notavailable as RPM

edhat you were enquiring about. All the best, G > > Scott Waichler > Pacific Northwest National Laboratory > scott.waichler _at_ pnl.gov -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Rese

Re: [R] Warning message when returning multiple items

1] 3 $squaredx [1] 81 HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] gavin.simpsonATNOSPAMucl.ac.uk UCL Depa

Re: [R] Prediction method for lowess,loess,lokerns,lpepa,ksmooth

bsite for something that meets your needs? HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] gavin.simpsonATNOSPAMucl.ac.uk UCL Department of Geography

Re: [R] symbol() function in 3d

.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 -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow

[R] Error when running configure on xeon with r-1-9-patched

cifics of what has gone wrong here could someone enlighten me as to my mistake? Many thanks in advance, Gavin -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7

Re: [R] Error when running configure on xeon with r-1-9-patched

x27;t have LD_LIBRARY_PATH set. libg2c.so is in /usr/lib/gcc-lib/i386-redhat-linux/3.3.3 and libg2c.so.0 and libg2c.so.0.0.0 both link from /usr/lib to libg2c.so in /usr/lib/gcc-lib/i386-redhat-linux/3.3.3. This is the same setup as the other Fedora boxes I have compiled R on. On Mon, 6 Sep 2004, Gav

Re: [R] Error when running configure on xeon with r-1-9-patched

Peter Dalgaard wrote: Gavin Simpson <[EMAIL PROTECTED]> writes: configure:25282: gcc -o conftest -D__NO_MATH_INLINES -mieee-fp -O2 -g -pipe -march=i386 -mcpu=i686 -I/usr/local/include -L/usr/local/lib conftest.c -lreadline -ldl -lncurses -lm -L/usr/local/lib -L/usr/lib/gcc-lib/i386-redhat

[R] generalized eigenvalues - porting code using eig() from Matlab to R

,] 0.000 0.000 0.000 0.2307692 How would I go about doing this in R? Thanks in advance, Gavin [1] http://www.mathworks.com/access/helpdesk/help/techdoc/ref/eig.html?cmdname=eig -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T

Re: [R] biplot.princomp with loadings only

All the best, Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] [EMAIL PROTECTED] UCL Department of Geography

Re: [R] Strange Matrix Multiplication Behaviour

expect both vectors to be equivalent. Hopefully my longwinded explanation helps. Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC

Re: [R] Performing Analysis on Subset of External data

that can connect to your database and only extract the rows you need. See the R Data Import/Export manual for more on scan and using relational databases with R. Hope this helps, Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T

Re: [R] Equivalents of Matlab's 'find' and 'end'

Y [1] 2 1 1 2 3 Not quite the same I'll grant you, but results in the same thing. If you wanted Y = X(4:end) you could use: > Y <- x[-c(1:3)] > Y [1] 1 2 3 HTH Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson

Re: [R] Number of characters per line

(The valid values are in file 'Print.h' and can be changed by re-compiling R.) Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 E

Re: [R] Selecting from a character vector with logical operators

ot meaningful for factors in: Ops.factor(letfact, "j") You need an ordered factor for this to work: > letfact <- as.ordered(letters) > letfact[letfact > "j"] [1] k l m n o p q r s t u v w x y z 26 Levels: a < b < c < d < e < f < g < h < i

[R] sub- and superscript in plot labels

u, "eq cm"^{-2}, " yr"^{-1}))) Strangely, I can do this: plot(1:10, xlab = expression(paste("nm SO"[4]^2, " ", mu, "eq cm"^{-2}, " yr"^{-1}))) With almost the desired effect (except I need to add two characters to the superscript). Can

Re: [R] sub- and superscript in plot labels

Gavin Simpson wrote: Dear List, I need to add a subscript and a superscript to some of the ions in the labels on some plots. I have got to here but now I'm stuck: plot(1:10, xlab = expression(paste("nm SO"[4], " ", mu, "eq cm"^{-2}, " yr"^{-1})))

Re: [R] sub- and superscript in plot labels

Prof Brian Ripley wrote: On Thu, 4 Nov 2004, Gavin Simpson wrote: Dear List, I need to add a subscript and a superscript to some of the ions in the labels on some plots. I have got to here but now I'm stuck: plot(1:10, xlab = expression(paste("nm SO"[4], " ", mu,

Re: [R] R works on Fedora Core 3

that command a few times, closing the opened faulty device window in between to get the title bar to appear. Same result running R from within xemacs/ess. Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5

Re: [R] time series plots: labels language & tickmarks

you have told > us your platform we could have told you how. If you had shown us some > reproducible code we would not have had to guess as to what you want to > do.) > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson

[R] confidence interval for StructTS?

time-series [e.g. from tsSmooth]). predict.StructTS() has the option of getting the standard errors of forecasts. I'd like to get the same but for observed series. Thanks in advance, Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Si

Re: [R] chisq.,test`

ing guide! http://www.R-project.org/posting-guide.html -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] gavin.sim

Re: [R] How do I get the row indices?

___ > 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 -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson

Re: [R] indicator value in labdsv

help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & E

Re: [R] (no subject)

ailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fel

[R] merging dataframes

order of rows the same, so that rows from mat1 come before the rows of mat2? Many thanks, Gavin -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC [E] gavin.simpsonAT

Re: [R] merging dataframes

ust pay more attention in class. Roger's and Andy's approaches seem like the most fool-proof way of canning this into a function that does all the I asked and a bit of tidying up of NA's etc. Cheers, G > > > From: Gavin Simpson > > > > Dear List, > >

[R] plots in a matrix

I'm sure there is something simple I'm missing, but I've tried various combinations of split.screen(), layout(), changing oma, mar but I still can't quite get this right. Any suggestions? Thanks in advance, Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] garch function in R

ckage: zoo > ?garch #displays help file as expected > > Which package(s) do I really need to run garch models? Looks like you need quadprog and zoo to be installed and available as well. > Thanks, > Shelton HTH, Gav -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] loadings matrices in plsr vs pcr in pls pacakage

; https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. & ECRC

Re: [R] Why is the R mailing list so hard to figure out?

est posts, for example, is at (via the "threaded HTTP" link in the Groups section of the above mentioned link): http://news.gmane.org/gmane.comp.lang.r.general I personally let my emailer manage the posting for me as I find that the easiest way to work... HTH G -- %~%~%~%~%~%~%~%~

Re: [R] Spectral analysis

ow and there isn't a ready built package on CRAN so you'll have to compile it yourself. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC [f] +44 (0)20 7

Re: [R] Formating the data

1 Is that what you want? It is just altering how the data are printed. You still get the [1] at the start though. G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography,

Re: [R] How do you do an e-mail post that is within an ongoing thread?

, but many people on the list do and they end up in the correct place in the thread. Note that in Gmail to see the headers for a message you can select "show original" from the little drop down menu (down triangle) next to the reply button. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%

[R] Efficiently calculate sd on an array?

m = TRUE, dims = 2)) instead of temp <- apply(arr, c(2,1), mean, na.rm = TRUE) but I am having difficulty seeing how to calculate the standard deviations efficiently. Any idea how I might go about this? All the best, G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] BIC and Hosmer-Lemeshow statistic for logistic regression

og(n), where n is number of observations, will get BIC. See ?AIC. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Buildin

Re: [R] How to extract diagonals

g) > all.equal(diags, diag(as.matrix(dis.bc))) [1] TRUE So you'll have to reformulate your question if this is not what you wanted. A word of warning, do not do diag(dis.bc)) on the above as it brought my Linux box to it's knees trying to do something silly - easily recoverable, but be

Re: [R] Dissimilarity

rning over symbol.For is an issue with the labdsv > > library itself > > - nothing you are doing wrong. The results will still be valid - > > but the use of > > symbol.For is something that will eventually need to be changed in > > the labdsv > > li

Re: [R] How to extract diagonals

is a learning experience. You only need to grovel and apologise if you have not done your homework before posting and not read the FAQ, the documentation or searched the archives, or followed the posting guide. Which is not the case here. HTH G > > Am 20.06.2007 um 14:11 schrieb Gavin Simpson

Re: [R] How to extract diagonals

my turn to learn about it. > > > Once again: thanks > > > Greetings > > > Birgit > > > > Am 20.06.2007 um 18:02 schrieb Gavin Simpson: > > > On Wed, 2007-06-20 at 15:09 +0200, Birgit Lemcke wrote: > > > Hello Gavin and thank

Re: [R] Distance function

deleted]] > > __ > 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, s

[R] Names of objects passed as ... to a function?

y returns a single "name": > foo <- function(...) + deparse(substitute(...)) > dat1 <- rnorm(10) > dat2 <- runif(10) > foo(dat1, dat2) [1] "dat1" Can anyone suggest to me a way to get the names of objects passed as the ... argument of a function? TIA G --

Re: [R] Names of objects passed as ... to a function?

On Sat, 2007-06-23 at 16:52 +0100, Prof Brian Ripley wrote: > On Sat, 23 Jun 2007, Gavin Simpson wrote: > > > Dear list, > > > > I have a function whose first argument is '...'. Each element of '...' > > is a data frame, and there will be

Re: [R] aggregating daily values

4, , 5556, 5557, > 5558, 5559, 5560, 5561, 5562, 5563, 5564, 5565, 5566, 5567, 5568, > 5569, 5570, 5571, 5572, 5573, 5574, 5575, 5576, 5577, 5578, 5579 > ), class = "Date"), class = "zoo") > > TIA > > Antonio > > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] using self-written functions

_ > 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] Fees to use R

On Fri, 2007-07-06 at 09:31 +0200, [EMAIL PROTECTED] wrote: > Good morning to all, > > I work for a bank in Italy, I want to know if i can install R and > relative add on like Rbloomberg for free or my company has to pay some > fee. > tanks to all. > Stefano Colucci R is released under the GNU GP

Re: [R] random sampling with some limitive conditions?

. > > > How to get the new random sampling data? I have no idea. > > > Thanks. > > > > __ > > R-help@stat.math.ethz.ch mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting gui

Re: [R] why function 'sum' within 'aggregate' need an 'index'?

2007 month 07 day05 svn rev42131 language R version.string R version 2.5.1 Patched (2007-07-05 r42131) -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson

Re: [R] distance function (analogue)

, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.500) Dist.Gower <- distance(Table1, Table0, method ="mixed", weights = my.weights) -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] Alternative to xyplot()?

On Tue, 2007-07-17 at 16:39 -0400, Manuel Morales wrote: > Hello list, > > Is anyone aware of a non-lattice-based alternative to xyplot()? x <- rnorm(20) y <- rnorm(20) plot(x, y) ? If you mean some specific aspect of xyplot(), you'll have to tell us what this is. HTH G __

Re: [R] Subsetting Enigma: More rows after dataframe[-list,]?

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. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL G

Re: [R] dates() is a great date function in R

4425 34547 34759 34773 > > The dates() function makes calendar dates from Julian dates, > > >dmp<-dates(ss$MPdate,origin=c(month = 1, day = 1, year = 1900)) > > > dmp[1:5] > [1] 10/12/93 04/03/94 08/03/94 03/03/95 03/17/95 > > I would appreciate the comments of

<    1   2   3   4   >