RE: [R] regression on a matrix

2005-03-04 Thread Martin Maechler
ReidH == Huntsinger, Reid [EMAIL PROTECTED] on Thu, 3 Mar 2005 17:24:22 -0500 writes: ReidH You might use lsfit instead and just do the whole Y ReidH matrix at once. That saves all the recalculation of ReidH things involving only X. yes, but in these cases, we have been

[R] Basic stratification calculations

2005-03-04 Thread dvrecko
Hi. I'm a student at SFU in Canada. The basic thing I want to do is calculate means of different strata. I have 2 vectors. One has the values I want to take the means from, the other is the four strata I am interested in. So I essentially want to break up the information vector into the four

Re: [R] Basic stratification calculations

2005-03-04 Thread vito muggeo
Hi, if I understand correctly, tapply() is your friend here, vito [EMAIL PROTECTED] wrote: Hi. I'm a student at SFU in Canada. The basic thing I want to do is calculate means of different strata. I have 2 vectors. One has the values I want to take the means from, the other is the four strata I am

Re: [R] Simple suggestion for improvement

2005-03-04 Thread Yan Wong
On 3 Mar 2005, at 17:17, Adaikalavan Ramasamy wrote: How will you deal with multiple word searches such as help.search(eps dev) One way to implement would be ??eps dev but this looks awkward to me. That's what you have to do with the normal help function sometimes anyway, e.g. ?+

Re: [R] Reconstructing Datasets

2005-03-04 Thread Angel Lopez
Hi Laura, You might want to have a look at function decevf in package pastecs. It uses eigenvector filtering to reconstruct a signal using only the most representative eigenvectors. It is applied for time series but you could easily modify the code to use it for spatial data also. Bests, Angel

RE: [R] ESS

2005-03-04 Thread Rau, Roland
Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mahdi Osman Sent: Thursday, March 03, 2005 11:57 PM To: r-help@stat.math.ethz.ch Subject: [R] ESS I can not start R proccess or ESS from within XEmacs. What is going wrong? Do you have a

[R] how to draw graphs within clickable hyperlink

2005-03-04 Thread Hu Chen
Hi all I want to draw a graph containing points and edges. Package graphics could do this. However I want to make each those points and edges clickable with a hyperlink on them. Anybody know how to do this? Thanks very much. __ R-help@stat.math.ethz.ch

[R] R: simulation

2005-03-04 Thread Clark Allan
hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there some arguement that goes with rnorm? please supply an

[R] bayesmix - What is or where can I find JAGS executable?

2005-03-04 Thread Karl Knoblick
Hallo! I want to use the package bayesmix. Trying the examples I had no success. The reason is, I think: haveJAGS() [1] FALSE Have read of JAGS executable (I could not find JAGS.exe in my R directory). What is JAGS? Where can I find or download it? Karl

[R] reproducibility of the loess function

2005-03-04 Thread Marie-Agnès DILLIES
Hi, We noticed a difference between R versions 1.9.1 and 2.0.1 concerning the loess function. The resulting fitted data is quite different. Furthermore, with the 2.0.1 release, if I apply the loess function several times on the same input data set, I get different results. Has anybody already

[R] Concatenate vector into string

2005-03-04 Thread Matthieu Cornec
Hello, I would like to convert c(a,b,c) into abc. Anyone could help? Thanks, Matthieu Cornec __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

RE: [R] regression on a matrix

2005-03-04 Thread Liaw, Andy
From: Martin Maechler ReidH == Huntsinger, Reid [EMAIL PROTECTED] on Thu, 3 Mar 2005 17:24:22 -0500 writes: ReidH You might use lsfit instead and just do the whole Y ReidH matrix at once. That saves all the recalculation of ReidH things involving only X. yes, but in

[R] R 2.0.1 installation

2005-03-04 Thread Ramseyer Amandine
Hello, I'm a student in biology, writting from Strasbourg, France. I need to use the logiciel R to analyse biological results. I have ADE4 since 2001, but technology's now largely evoluted ! That's why I'm just trying to get the new version of : - ADE4_1.3-3.zip ; - R 2.0.1.ter

Re: [R] Basic stratification calculations

2005-03-04 Thread Clint Harshaw
vito muggeo wrote: Hi, if I understand correctly, tapply() is your friend here, vito [EMAIL PROTECTED] wrote: Hi. I'm a student at SFU in Canada. The basic thing I want to do is calculate means of different strata. I have 2 vectors. One has the values I want to take the means from, the other is

Re: [R] ESS

2005-03-04 Thread Adaikalavan Ramasamy
I did not suggest the use of emacs over Xemacs (or vice versa) but was merely suggesting that I have little experience with Xemacs and thus may advice may be inaccurate to some level. I do not know why Stata does not work for you. I never used Stata and this is neither a Stata nor ESS mailing

Re: [R] R: simulation

2005-03-04 Thread Dimitris Rizopoulos
look at ?set.seed, e.g., rnorm. - function(seed, ...){ set.seed(seed) rnorm(...) } # rnorm.(100, 1, 2) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35,

Re: [R] R: simulation

2005-03-04 Thread Lukas Meier
Clark Allan writes: i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there some arguement that goes with rnorm? set.seed(79) ##

Re: [R] R: simulation

2005-03-04 Thread Clark Allan
Hi thanx for the reply. i used your code and pasted it into R and ran it a few times. the output is below. what i want is to get the same output every time the program is run. is this possible? another question? x-rnorm(100) y-rnorm(100) is x and y independent? rnorm(1) [1] 0.4251004

Re: [R] R: simulation

2005-03-04 Thread Petr Pikal
On 4 Mar 2005 at 13:07, Clark Allan wrote: hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? See ?set.seed It is

Re: [R] R: simulation

2005-03-04 Thread Detlef Steuer
On Fri, 04 Mar 2005 13:07:21 +0200 Clark Allan [EMAIL PROTECTED] wrote: hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some

Re: [R] R: simulation

2005-03-04 Thread Dimitris Rizopoulos
are you sure you tried the function I gave you? There is a . that you might have missed. When I try it I get this: rnorm. - function(seed, ...){ + set.seed(seed) + rnorm(...) + } # rnorm.(1, 1) [1] -0.6264538 rnorm.(1, 1) [1] -0.6264538 rnorm.(1, 1) [1] -0.6264538 rnorm.(1, 1)

Re: [R] bayesmix - What is or where can I find JAGS executable?

2005-03-04 Thread Bettina Gruen
Hi Karl! In the DESCRIPTION file there is the package homepage (www.ci.tuwien.ac.at/~gruen/BayesMix) indicated. You can download the windows binaries together with JAGS.exe from this page or follow the link to Martyn Plummer's JAGS homepage where you can get the sources for compiling jags for

Re: [R] Concatenate vector into string

2005-03-04 Thread Chuck Cleland
Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? ?paste paste(c(a,b,c), collapse=) -- Chuck Cleland, Ph.D. NDRI, Inc. 71 West 23rd Street, 8th floor New York, NY 10010 tel: (212) 845-4495 (Tu, Th) tel: (732) 452-1424 (M, W, F) fax: (917) 438-0894

Re: [R] R: simulation

2005-03-04 Thread Uwe Ligges
Clark Allan wrote: hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there some arguement that goes with rnorm? please

Re: [R] bayesmix - What is or where can I find JAGS executable?

2005-03-04 Thread Uwe Ligges
Karl Knoblick wrote: Hallo! I want to use the package bayesmix. Trying the examples I had no success. The reason is, I think: haveJAGS() [1] FALSE Have read of JAGS executable (I could not find JAGS.exe in my R directory). What is JAGS? Where can I find or download it? What about googling for

RE: [R] Concatenate vector into string

2005-03-04 Thread Bashir Saghir (Aztek Global)
Try: paste(c(a,b,c), collapse=) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Matthieu Cornec Sent: Friday, March 04, 2005 12:32 To: [EMAIL PROTECTED] Subject: [R] Concatenate vector into string Hello, I would like to convert c(a,b,c) into abc.

Re: [R] Concatenate vector into string

2005-03-04 Thread Uwe Ligges
Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? See ?paste Uwe Ligges Thanks, Matthieu Cornec __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Rank-based p-value on large dataset

2005-03-04 Thread Sean Davis
On 3/3/05 17:40, Deepayan Sarkar [EMAIL PROTECTED] wrote: On Thursday 03 March 2005 16:32, Deepayan Sarkar wrote: On Thursday 03 March 2005 16:22, Sean Davis wrote: I have a fairly simple problem--I have about 80,000 values (call them y) that I am using as an empirical distribution and I want

Re: [R] R: simulation

2005-03-04 Thread Adaikalavan Ramasamy
See help(set.seed). set.seed(1) rnorm(5) [1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078 rnorm(5) [1] -0.8204684 0.4874291 0.7383247 0.5757814 -0.3053884 set.seed(1) rnorm(5) [1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078 On Fri, 2005-03-04 at 13:07 +0200, Clark

Re: [R] Concatenate vector into string

2005-03-04 Thread Adaikalavan Ramasamy
Use the collapse argument in paste. paste( c(a, b, c), collapse= ) [1] abc On Fri, 2005-03-04 at 03:32 -0800, Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? Thanks, Matthieu Cornec __

Re: [R] R: simulation

2005-03-04 Thread Christian Hennig
?set.seed On Fri, 4 Mar 2005, Clark Allan wrote: hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there

Re: [R] Concatenate vector into string

2005-03-04 Thread Pierre BADY
Hi all, ted - letters[1:3] ted paste(ted, collapse = ) cheers, P.BADY At 03:32 04/03/2005 -0800, Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? Thanks, Matthieu Cornec __ R-help@stat.math.ethz.ch

[R] lm and time series

2005-03-04 Thread Matthieu Cornec
Hello, I create a multivariate time series containing NA values (that could come directly from an imported file,) I want to compute a linear regression and obtain a time serie for both residuals and fitted values. I have tried the trick ts.intersect, without success. Could you help me out of

[R] how to draw graphs within clickable hyperlink

2005-03-04 Thread Damian Betebenner
Hu, I've found two methods of doing this---1 in R and the other in S-Plus. 1. The S-Plus solution is called graphlets which are a java-based graphics device that allows one to build various types of interactivity into their graphs. 2. The R solution uses the gridSVG package by Paul Murrell.

RE: [R] R: simulation

2005-03-04 Thread Liaw, Andy
See ?set.seed. Andy From: Clark Allan hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there some

Re: [R] Basic stratification calculations

2005-03-04 Thread Sean Davis
On Mar 4, 2005, at 6:55 AM, Clint Harshaw wrote: vito muggeo wrote: Hi, if I understand correctly, tapply() is your friend here, vito [EMAIL PROTECTED] wrote: Hi. I'm a student at SFU in Canada. The basic thing I want to do is calculate means of different strata. I have 2 vectors. One has the

RE: [R] how to draw graphs within clickable hyperlink

2005-03-04 Thread Liaw, Andy
You need to output the graph in a format that allow such interactivity. I believe svg (available through the svgDevice package on CRAN) can do that. Andy From: Hu Chen Hi all I want to draw a graph containing points and edges. Package graphics could do this. However I want to make each

RE: [R] Concatenate vector into string

2005-03-04 Thread Liaw, Andy
See ?paste, e.g., char - c(a, b, c) paste(char, collapse=) [1] abc Andy From: Matthieu Cornec Hello, I would like to convert c(a,b,c) into abc. Anyone could help? Thanks, Matthieu Cornec __ R-help@stat.math.ethz.ch mailing list

Re: [R] R: simulation

2005-03-04 Thread Spencer Graves
The help page for rnorm says See Also: ... '.Random.seed' ..., and the help page for .Random.seed says, 'set.seed' is the recommended way to specify seeds. hope this helps. spencer graves Clark Allan wrote: hi all a simple question i want to run simulations in r. i however

Re: [R] Concatenate vector into string

2005-03-04 Thread Witold Eryk Wolski
Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? Thanks, Matthieu Cornec __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] R: simulation

2005-03-04 Thread Peter Dalgaard
Clark Allan [EMAIL PROTECTED] writes: hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there some

Re: [R] bayesmix - What is or where can I find JAGS executable?

2005-03-04 Thread Spencer Graves
I did www.r-project.org - search - R site search - bayesmix. The first hit told me it was a package, which I confirmed by trying www.r-project.org - download: CRAN - (select a local mirror) - Software: Packages. There, I found bayesmix listed. Then I used

Re: [R] Concatenate vector into string

2005-03-04 Thread Spencer Graves
How about: paste( c(a,b,c), collapse=) [1] abc hope this helps. spencer graves Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? Thanks, Matthieu Cornec __ R-help@stat.math.ethz.ch

[R] Farey sequences

2005-03-04 Thread Robin Hankin
Hi has anyone coded up Farey sequences? [ The Farey sequence of order n is the set of rational numbers i/j with (i,j)=1 such that 0 = i,j = n; the sequence is ordered from lowest to highest. Thus Farey_4 = {0/1 , 1/4 , 1/3 , 1/2 , 2/3 , 3/4 , 1/1} ] My motivation is unimodular transformations:

Re: [R] Concatenate vector into string

2005-03-04 Thread Tim Cutts
On 4 Mar 2005, at 11:32 am, Matthieu Cornec wrote: Hello, I would like to convert c(a,b,c) into abc. Anyone could help? paste(c(a, b, c), sep=) Tim -- Dr Tim Cutts Informatics Systems Group, Wellcome Trust Sanger Institute GPG: 1024D/E3134233 FE3D 6C73 BBD6 726A A3F5 860B 3CDD 3F56 E313 4233

Re: [R] R 2.0.1 installation

2005-03-04 Thread Adaikalavan Ramasamy
Are you using Windows operating system ? If so, then you will need to download the executable not the source codes from http://www.cran.r-project.org/bin/windows/base/rw2001.exe Regards, Adai On Fri, 2005-03-04 at 12:54 +0100, Ramseyer Amandine wrote: Hello, I'm a student in biology,

[R] Clustering of Binary data in R

2005-03-04 Thread sima fakheran
Good afternoon! I would like to ask you about similarity measures and clustering in R for Binary data. Would you please kindly help me and let me know about that commands in R? Thanks in advance for your kind attentions. I look forward to hearing from you as soon as possible.

RE: [R] lm and time series

2005-03-04 Thread Gabor Grothendieck
From: Matthieu Cornec [EMAIL PROTECTED] I create a multivariate time series containing NA values (that could come directly from an imported file,) I want to compute a linear regression and obtain a time serie for both residuals and fitted values. I have tried the trick ts.intersect,

Re: [R] calculating of linkage-disequilibrium measures?

2005-03-04 Thread Scott D Chasalow
Taking a lead from Brian, the answer is yes. :-) But just in case you were seeking a somewhat more detailed answer: If phase is known, it is not only possible to compute all these LD measures in R, it is much simpler. For example, if you have alleles A and a at locus 1, and B and b at locus 2,

Re: [R] how to draw graphs within clickable hyperlink

2005-03-04 Thread Barry Rowlingson
Damian Betebenner wrote: Hu, I've found two methods of doing this---1 in R and the other in S-Plus. Here's something I did a few years ago that lets you create PNG files and an HTML imagemap http://www.maths.lancs.ac.uk/Software/Imagemap/ It might be tricky to add clickable areas to plots drawn

Re: [R] Clustering of Binary data in R

2005-03-04 Thread Christian Hennig
Hi, as a starting point, see ?dist (method=binary) library(prabclus) ?jaccard ?kulczynski For distance based clustering methods see library(cluster) ?agnes ?pam Best, Christian On Fri, 4 Mar 2005, sima fakheran wrote: Good afternoon! I would like to ask you about similarity

[R] gbm

2005-03-04 Thread WeiWei Shi
Hi, there: Is there anyone who read the codes for gbm package before? Before i sent this email, I also sent an email to ask for help from the author, Greg. But still I am wondering if someone here can share some understanding like the roadmap or document on the implementation too. Thanks, Ed.

[R] test of significance for nlme coefficients

2005-03-04 Thread Lindsey Root
Hello, I am using the nlme package to fit an exponential decay model to a longitudinal data set: ##Model a variable called RASCH16A as an exponential function ## of an initial value A and a decay rate K fm2 - nlme(model = RASCH16A ~ A* (exp(-K*DAYS)), + ###A has nonzero fixed and random

[R] Listas de E-mails para Mala Direta , Marketing Direto

2005-03-04 Thread marcosjbd_g5a
cadastros de e-mails divididos por profissão cadastros de e-mails divididos por segmento Cadastros para mala direta livre de spam: http://www.gueb.de/divulgamail Listas de E-mails para Mala Direta , Marketing Direto cadastros de e-mails divididos por profissão cadastros de e-mails divididos

[R] Need suggestions for finding dose response using nls

2005-03-04 Thread Ken
I am relatively new to R and am looking for advice, ideas or both... I have a data set that consists of pathogen population sizes on individual plant units in an experimental field plot. However, in order to estimate the pathogen population sizes I had to destroy the plant unit and could not

[R] MGCV-confidence interval

2005-03-04 Thread Aloisio Joaquim Freitas Ribeiro
I am using MGCv to smoothing mortality rates. I need to= construct simultaneous confidence intervals for the mortality function. Is= It a possible using MGCV? Thanks Aloisio -- Open WebMail Project ([1]http://openwebmail.org) -- Esta mensagem foi verificada

[R] R-2.01 and RSPerl-0.6.2

2005-03-04 Thread Zachariah Zachariah
Hi, I am somewhat new to R and RSPerl, but I think this particular problem has to do with RSPerl and so I am not sure if this is the right forum to ask for help. Nevertheless I am quite sure that many of you would have used RSPerl with R. My hardware platform is a Sun/Solaris V440 server running

[R] R: Moran's I

2005-03-04 Thread Theresa Sinicrope Talley
Hi- I would like to assign distance classes (lags) to the moran's I test in R. I don¹t have equally spaced points since I was mapping shrubs across a landscape (vs. an experiment or survey in a grid). I somehow need to select points based on the distances in the distance matrix (i.e., I need

Re: [R] R-2.01 and RSPerl-0.6.2

2005-03-04 Thread Duncan Temple Lang
The omega-help mailing list is more appropriate. (See http://www.omegahat.org/mailman/listinfo) One of the things that comes to mind is whether you built R via the --enable-R-shlib entirely from scratch or did you reconfigure and rebuild from an existing compiled source. If it is the latter,

[R] here's why it make s sense

2005-03-04 Thread Paul Johnson
if you go x[i] you are giving x an index vector, which we had mistakenly thought was an integer. Rather, it is a vector of indices for observations. Here's data x - c(1 , 4, 3, 2, 5) x[1] would be 1 x[2] would bd 4 but if you put an index vector in the brackets, you have x [ c(1,2,1,2] ] it

[R] vector memory allocation?

2005-03-04 Thread Sam Yeaman
Hi all, I have a vector size allocation problem with R 2.0.1 (script and output shown): var1 - sum (input1 * input2, na = TRUE) gc() used (Mb) gc trigger (Mb) Ncells 199327 5.4 785113 21.0 Vcells 71039552 542.0 206003790 1571.7 var2 - sum (input1 * input2 / input2,

[R] Emacs keystroke to toggle T/F for setting logical values

2005-03-04 Thread Waichler, Scott R
I'd like to have an Emacs keystroke that would let me toggle between T and F when editing logical settings in R code. I looked in the ESS documentation and in my O'Reilly emacs books but found nothing. Any ideas? Scott Waichler Pacific Northwest National Laboratory [EMAIL PROTECTED]

[R] Gap plots

2005-03-04 Thread Liu, Jane
Hi, Does anyone know how to do gap plots for k-mean clustering in R? Thanks a lot! Jane __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: leap years (Was: [R] lm and time series)

2005-03-04 Thread Kjetil Brinchmann Halvorsen
A followup: How do people treat dsaily time series, when there is a yearly cycle? For the moment I just left out the 29 of February's, but that dsoes'nt look really good. Is the only way to include them to use irregular time series? Kjetil Gabor Grothendieck wrote: From: Matthieu Cornec [EMAIL

Re: leap years (Was: [R] lm and time series)

2005-03-04 Thread Gabor Grothendieck
The other approach is the pastecs::daystoyears approach where each year consists of 365.25 days. Kjetil Brinchmann Halvorsen kjetil at acelerate.com writes: : : A followup: : : How do people treat dsaily time series, when there is a yearly cycle? : : For the moment I just left out the 29 of

[R] S-code for piecewise regression

2005-03-04 Thread Simon Chamaillé
Dear R-helpers, a S-code for piecewise regressions was provided by Toms Lesperance (2003) Ecology, 84, 2034-2041 (paper can be found on the web).The code is quite complete with different types of transitions around breakpoints and model selection fonctions. It doesn't work directly under R due

Re: [R] vector memory allocation?

2005-03-04 Thread Prof Brian Ripley
On Fri, 4 Mar 2005, Sam Yeaman wrote: I have a vector size allocation problem with R 2.0.1 (script and output shown): And your computer's OS and RAM size are? var1 - sum (input1 * input2, na = TRUE) What do you think na=TRUE does? There is an na.rm argument, but sum(1, na=TRUE) may surprise you.