RE: [R] How to write the output of a function into a file?

2003-03-26 Thread Gavin Simpson
does sink() do what you want? see ?sink for more #divert printed results to file sink(file = output.txt) #commands here chisq.test(x,p=probs) #turn off sink sink(file = NULL) If you are not typing the commands at the prompt but are using an external editor or text file for your commands

RE: [R] Re: Bar plot with variable width (down a drill hole) - n

2003-03-26 Thread Ted Harding
On 25-Mar-03 Phillip J. Allen wrote: Thanks Ted and Marc its works. But of course after pulling in in some real life data I discoverd one hitch. Often there are missing intervals. For example: from - c(0, 1.2, 4.0, 4.2, 5.0, 25.0, 30.1, 45) to - c(1.2, 4.0, 4.2, 5.0, 25, 30.1, 36.2, 50)

[R] ifelse()

2003-03-26 Thread Ko-Kang Kevin Wang
Hi, I'm not sure if this can be done but.. I know that with ifelse() I can do something like: ifelse(x = 3, 1, 2) to go through each element in my vector x, and if x_i = 3 substitute the number with 1 else with 2. Essentially I'll get a vector with 2 levels. Can I tweak it so I can get

Re: [R] ifelse()

2003-03-26 Thread Jason Turner
On Wed, Mar 26, 2003 at 11:20:04PM +1200, Ko-Kang Kevin Wang wrote: Hi, I'm not sure if this can be done but.. I know that with ifelse() I can do something like: ifelse(x = 3, 1, 2) to go through each element in my vector x, and if x_i = 3 substitute the number with 1 else with 2.

[R] Re: ifelse()

2003-03-26 Thread Ko-Kang Kevin Wang
Hi, I've found a much better solution than using ifelse(). I found about cut() from MASS4 Page 383, which actually does a better job . On Wed, 26 Mar 2003, Ko-Kang Kevin Wang wrote: Date: Wed, 26 Mar 2003 23:20:04 +1200 (NZST) From: Ko-Kang Kevin Wang [EMAIL PROTECTED] To: R Help [EMAIL

RE: [R] Re: Bar plot with variable width (down a drill hole) - n

2003-03-26 Thread Ted Harding
On 26-Mar-03 Ted Harding wrote: On 25-Mar-03 Phillip J. Allen wrote: However, this kind of situation needs thought about alternative ways of representing it. One possibility might be to have the vertical axis invisible, so that gaps in the data are represented by gaps in the axis (and then

Re: [R] ifelse()

2003-03-26 Thread Martin Maechler
Jason == Jason Turner [EMAIL PROTECTED] on Wed, 26 Mar 2003 22:49:27 +1200 writes: Jason On Wed, Mar 26, 2003 at 11:20:04PM +1200, Ko-Kang Jason Kevin Wang wrote: Hi, I'm not sure if this can be done but.. I know that with ifelse() I can do something like:

RE: [R] RODBC and Excel in Widows

2003-03-26 Thread r.ghezzo
HI, no sorry, so far nobody answer. So it probably does not have a solution. Excell is from you.know.who = Original Message From Meinhard Ploner [EMAIL PROTECTED] = Hello! Did you resolve the problem? I'm interested in the solution, too. Meinhard On Thursday, March 13, 2003, at 07:21

[R] CMDS

2003-03-26 Thread Heeringa W.J.
Dear statisticians, I have a small question about the Classifical Multidimensional scaling routine in R. Is this procedure always metric? Best regards, Wilbert Heeringa __ [EMAIL PROTECTED] mailing list

Re: [R] RODBC and Excel in Widows

2003-03-26 Thread David Meyer
You might look at Thomas Baier's DCOM interface as an alternative to the odbc-method for accessing EXCEL-files. -d r.ghezzo wrote: HI, no sorry, so far nobody answer. So it probably does not have a solution. Excell is from you.know.who = Original Message From Meinhard Ploner [EMAIL

RE: [R] RODBC and Excel in Widows

2003-03-26 Thread ripley
Why don't you debug the code and send the patch needed to the maintainer? R is supposed to be a collaborative project, and we need more contributors, especially on Windows. Looking forwards to your positive contribution On Wed, 26 Mar 2003, r.ghezzo wrote: HI, no sorry, so far nobody

[R] Spatial Clustering

2003-03-26 Thread Philippe Hupé
Hello, Does anyone know any package which can perform spatial clustering ? Thanks. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] Interest Calculator for Lawyers

2003-03-26 Thread David James Thorpe JD
Have you ever been in this situation? You are representing a client who is supposed to receive a monthly payment. On some months, that payment has been made on time, but on other months it has been paid late and often not at all. Your client is entitled to interest on amounts owing, but figuring

[R] Simple NLME problem (I hope)

2003-03-26 Thread Vumani Dlamini
Dear R-Users: I would like to fit a multilevel model using LME such that the parameters after fitting the multilevel in two separate groups and when I use the complete data (with interactions between the grouping variable and the other variables) set are comparable (or the same). The problem

[R] 2 scatter plots in the same graphic window....

2003-03-26 Thread Vasudevan, Geetha
Hello, I am trying to do a scatter plot of x1,y1 and x2,y2 in the same graphics window with diff col/pch values. i cannot get it to do them in the same window. what is the cmd to do this? thanks. __ [EMAIL PROTECTED] mailing list

RE: [R] 2 scatter plots in the same graphic window....

2003-03-26 Thread Henrik Bengtsson
plot() and points(), e.g. x - matrix(rnorm(300), ncol=3) plot(x[,1],x[,2], pch=1, col=red) points(x[,1],x[,3], pch=2, col=blue) Cheers Henrik Bengtsson Lund University -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Vasudevan, Geetha Sent: den

Re: [R] 2 scatter plots in the same graphic window....

2003-03-26 Thread Ben Bolker
Either n1 - length(x1) n2 - length(x2) plot(c(x1,x2),c(y1,y2),col=rep(c(1,2),c(n1,n2)), pch=rep(c(1,2),c(n1,n2)) OR plot(x1,y1,xlim=range(c(x1,x2)),ylim=range(c(y1,y2)),col=1,pch=1) points(x2,y2,col=2,pch=2) will work. Ben Bolker On Wed, 26 Mar 2003, Vasudevan, Geetha wrote: Hello,

[R] hist overlay...

2003-03-26 Thread Vasudevan, Geetha
thanks to all for the 2d scatter plot. i have one more. how do i plot 'hist(y1, col=red) and hist(y2,col=blue) in the same window? thanks again. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] a statistic question about chisq.test()

2003-03-26 Thread aprilsun
Hi, In the chisq.test(), if the expected frequency for some categories is 5, there will be a warning message which says Warning message: Chi-squared approximation may be incorrect in: chisq.test(x, p = probs) I am wondering whether there are some methods to get rid of this mistake... Seems the

Re: [R] hist overlay...

2003-03-26 Thread Jerome Asselin
Do you mean two separate plots in the same window? If so, read the documentation on multiple figure environment in the chapter on Graphical Procedures in the manual An Introduction to R. Alternately, reading plot help file will lead you to the par help file, which will enable you to do

RE: [R] hist overlay...

2003-03-26 Thread Vasudevan, Geetha
Cool, thanks! -Original Message- From: Ross Ihaka [mailto:[EMAIL PROTECTED] Sent: Wed 3/26/2003 12:26 PM To: Vasudevan, Geetha Cc: [EMAIL PROTECTED] Subject:Re: [R] hist overlay... Vasudevan, Geetha wrote: thanks to all for the 2d scatter plot. i have one more.

[R] Categorical data (Monte Carlo exact p-value)

2003-03-26 Thread Jorge Magalhes
In R program, I can perform categorical data test analysis like Odds ratio test in stratified 2x2 contingency tables? I do that in statistical package StatXact, but i would like perform the same test in R environment. Thanks very much Jorge Magalhães

[R] R TclTk iwidgets::comboboc

2003-03-26 Thread James Wettenhall
Hi, I am trying to create a drop-down combobox in R TclTk. The following works fine for a ListBox but fails for a combobox: # THIS WORKS FINE - CREATES AN EMPTY LISTBOX ## tt-tktoplevel() win - .Tk.subwin(tt) .Tcl(paste(listbox,.Tk.ID(win),.Tcl.args())) tkpack(win)

Re: [R] Categorical data (Monte Carlo exact p-value)

2003-03-26 Thread Spencer Graves
A search in the standard R documentation for Fisher's exact test reveals a command fisher.test, and a search for Chi-square reveals a chisq.test. Spencer Graves Jorge Magalhães wrote: In R program, I can perform categorical data test analysis like Odds ratio test in stratified 2x2

Re: [R] R TclTk iwidgets::combobox

2003-03-26 Thread James Wettenhall
Thanks Peter, I've started a new R session (in Windows) and managed to get both ways working now : ### THIS WORKS !!! ### library(tcltk) tclRequire(IWidgets) tt-tktoplevel() combo - tkwidget(tt,iwidgets::combobox) tkpack(combo) ### AND THIS WORKS TOO !!! ### tt-tktoplevel() win -

RE: [R] a statistic question about chisq.test()

2003-03-26 Thread Wiener, Matthew
One other option. I usually find that when I do the chisq.test with exact p-value calculation, I find the p-values are nearly identical to the results when I use the approximation and get the warnings (I'm usually dealing with just a few bins with less than 5, and many bins with more than 5). So

Re: [R] Re: Bar plot with variable width (down a drill hole) - nowmissing intervals?

2003-03-26 Thread Marc Schwartz
Phillip J. Allen wrote: Hi again, Thanks Ted and Marc its works. But of course after pulling in in some real life data I discoverd one hitch. Often there are missing intervals. For example: from - c(0, 1.2, 4.0, 4.2, 5.0, 25.0, 30.1, 45) to - c(1.2, 4.0, 4.2, 5.0, 25, 30.1, 36.2, 50)