Re: [R] Maps and plotting

2004-10-13 Thread Roger Bivand
On Wed, 13 Oct 2004, Shawn Way wrote: > At our facility we have multiple sample points that are sampled on any > given day. What I would like to do is create a map of the facility with > the sample points (and point labels) and when we have out of > specification results, place a transparent dot

Re: [R] Problems understanding qr. family of functions

2004-10-13 Thread Prof Brian Ripley
On Wed, 13 Oct 2004, Kjetil Brinchmann Halvorsen wrote: > Hola! > > By my understanding of reading ?qr > the following shold result a 3 x 3 matrix: (rw2000) > > x <- matrix( rnorm(12), 4,3) > y <- matrix( rnorm(12), 4,3) > xqr <- qr(x) > yqr <- qr(y) > > qr.qty( xqr, qr.Q(yqr) ) # dim (

Re: [R] Statistical analysis of a large database

2004-10-13 Thread Luis Torgo
On Tue, 2004-10-12 at 08:36, Prof Brian Ripley wrote: > > LuÃs Torgo, Data Mining with R. Learning by case > > studies, Maggio 2003 > > http://www.liacc.up.pt/~ltorgo/DataMiningWithR/ > > Please note that that reference is not about large datasets, nor about > `data mining' in the generally used

[R] Maps and plotting

2004-10-13 Thread Shawn Way
At our facility we have multiple sample points that are sampled on any given day. What I would like to do is create a map of the facility with the sample points (and point labels) and when we have out of specification results, place a transparent dot over the area on the map. As the number of OOS

Re: [R] Problems understanding qr. family of functions

2004-10-13 Thread Douglas Bates
Kjetil Brinchmann Halvorsen wrote: Hola! By my understanding of reading ?qr the following shold result a 3 x 3 matrix: (rw2000) x <- matrix( rnorm(12), 4,3) y <- matrix( rnorm(12), 4,3) xqr <- qr(x) yqr <- qr(y) qr.qty( xqr, qr.Q(yqr) ) # dim (3,3) [,1] [,2] [,3] [1,] 0.16

RE: [R] incomplete function output

2004-10-13 Thread Randy Johnson [Contr]
You can return a list of them all root <- function(var) { return(list(test1 = PP.test(var, lshort = T), test2 = ...)) } output <- root(var) Then you can print them, save them, or whatever. Or you could just print them out in the function root <- function(var) { #---Phillips

Re: [R] incomplete function output

2004-10-13 Thread Sundar Dorai-Raj
bogdan romocea wrote: Dear R users, I have a function (below) which encompasses several tests. However, when I run it, only the output of the last test is displayed. How can I ensure that the function root(var) will run and display the output from all tests, and not just the last one? Thank you, b

Re: [R] incomplete function output

2004-10-13 Thread Spencer Graves
Only the last appears, because the function returns only the last, which then gets printed. The answer to your question depends on whether you want to see the results or store them. To see them, try the following: root.print <- function(var) { #---Phillips-Perron print(PP.test(var,

RE: [R] incomplete function output

2004-10-13 Thread Austin, Matt
I would return the values from the various tests in a list. If you only want them to print and not for use in other parts of your program you could explicitly print each test using print(). root <- function(var) { #---Phillips-Perron test1 <- PP.test(var, lshort = TRUE) test2 <- PP.test(var, ls

[R] Problems with randomForest for regression

2004-10-13 Thread jmoreira
Dear list, I am trying to do a benchmark study for my case study. It is a regression problem. Among other models I use randomForest. Using the following code the result is around 0.628, and this make sense comparing with other methods. The Theil function implements Theil's U statistic. I do no

[R] incomplete function output

2004-10-13 Thread bogdan romocea
Dear R users, I have a function (below) which encompasses several tests. However, when I run it, only the output of the last test is displayed. How can I ensure that the function root(var) will run and display the output from all tests, and not just the last one? Thank you, b. root <- function(v

Re: [R] Why I can't retrieve GO identifier correctly?

2004-10-13 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hello, Uwe and Robert, Well, you are sending your message to Uwe and R-help rather than "Uwe and Robert" .. Thank you for your help! I fixed the problem as your suggestion. I run into a new issue: I can't sink the output to a file in the script within loops as bellow

[R] Problems understanding qr. family of functions

2004-10-13 Thread Kjetil Brinchmann Halvorsen
Hola! By my understanding of reading ?qr the following shold result a 3 x 3 matrix: (rw2000) x <- matrix( rnorm(12), 4,3) y <- matrix( rnorm(12), 4,3) xqr <- qr(x) yqr <- qr(y) qr.qty( xqr, qr.Q(yqr) ) # dim (3,3) [,1] [,2] [,3] [1,] 0.16815466 -0.1970936 0.2351670 [2,] 0

Re: [R] Why I can't retrieve GO identifier correctly?

2004-10-13 Thread szhan
Hello, Uwe and Robert, Thank you for your help! I fixed the problem as your suggestion. I run into a new issue: I can't sink the output to a file in the script within loops as bellows, but I can sink it in the command line( I run R1.9.1 on WinXP). The scrip: sink("C:/level3BPterm.txt") level2<-getG

RE: [R] Re:How to create a R -application

2004-10-13 Thread Kunal Shetty
Dear Ted and R- users and Helpers very sorry folks i take by my question.i have found the solution...it is Installing Packages To get a package, the easiest way is to use Rgui. In Rgui, at the top you will see the standard pull-down menus, click on Packages -> Install package from CRAN

RE: [R] Re:How to create a R -application

2004-10-13 Thread Kunal Shetty
Dear Ted and R- users and Helpers thanks ted once again for letting me know about Shafer's EM methods I am developing or rather working on windows platform. how do I use or rather install the contributed packages in my application or script. I could find documentation

RE: [R] Maximum Likelihood :- Log likehoood function

2004-10-13 Thread Ted Harding
On 13-Oct-04 Kunal Shetty wrote: > Dear R - users/Helpers > > I am dealing with bivariate Normal data with missing values. Further I > am trying to implement Expectation-Maximization (EM) algorithm to > resolve the missing data problem. > Now one of the requirements is use the Log likehood func

[R] Maximum Likelihood :- Log likehoood function

2004-10-13 Thread Kunal Shetty
Dear R - users/Helpers I am dealing with bivariate Normal data with missing values. Further I am trying to implement Expectation-Maximization (EM) algorithm to resolve the missing data problem. Now one of the requirements is use the Log likehood function i.e -2Log so as to find a reliable con

RE: [R] one more Rcmdr problem

2004-10-13 Thread John Fox
Dear Adrian, I've just been able to verify this problem. It occurs for me when I try to copy via Ctrl-C, but not when I use the Rcmdr Edit -> Copy menu, nor when I use the right-click context menu. As near as I can tell, the problem is general to all Ctrl-key combinations. I don't have time at the

Re: [R] data(eurodist) and PCA ??

2004-10-13 Thread Dan Bolser
On 13 Oct 2004, Jari Oksanen wrote: >On Wed, 2004-10-13 at 09:51, Prof Brian Ripley wrote: >> On Wed, 13 Oct 2004, Dan Bolser wrote: > >> > I have a complex distance matrix, and I am thinking about how to cluster >> > it and how to visualize the quality of the resulting clusters. >> >> Using PCA

[R] one more Rcmdr problem

2004-10-13 Thread Adrian Dusa
Hello, I'm using R 2.0.0 with the latest Rcmdr package installed from CRAN, on Windows XP Professional. When trying to copy some commands or results, either from the upper or lower text window, this causes Rcmdr to crash: "R for Windows GUI front-end has encountered a problem and needs to close"

Re: [R] data(eurodist) and PCA ??

2004-10-13 Thread Jari Oksanen
On Wed, 2004-10-13 at 09:51, Prof Brian Ripley wrote: > On Wed, 13 Oct 2004, Dan Bolser wrote: > > I have a complex distance matrix, and I am thinking about how to cluster > > it and how to visualize the quality of the resulting clusters. > > Using PCA and plotting the first two components is cl

RE: [R] Seeking advice on "introducing R"

2004-10-13 Thread John Fox
Dear Arin, I have a variety of materials at that may be helpful. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://so

[R] Seeking advice on "introducing R"

2004-10-13 Thread data Analytics
Dear List: This is in search of advices/opinions for a lecture presentation on R. I have this presentation scheduled on Saturday (16-October-2004); the topic is "introduction to R". The audience is a group of researchers and scholars who have either never used R/are beginning to use R, but othe

Re: [R] importing data

2004-10-13 Thread Prof Brian Ripley
On Wed, 13 Oct 2004, Andreas Betz wrote: > I am newcomer to R and I am loborating on this problem: > > How do I import data from a CD into R for further evaluation. > Using code "newgotcha <- read.table("e:\\asciiwin\\gotcha.dat") " > returmns a list containing an additional column with indices

[R] 'Replace' in internal editor causes R to crash

2004-10-13 Thread Luis Rideau Cruz
R-help, While using the internal editor in R (R 2.0.0) I have encountered the following problem : When I execute the facility 'Replace' (Ctrl H) R crashes (exactly when I paste to field 'look for') This is not one-off.I have faced it several days ago (and today as well). System specifications:

Re: [R] Cleaning up R-2.0.0 installation on GNU/Linux ?[SOLVED]

2004-10-13 Thread Prof Brian Ripley
On Wed, 13 Oct 2004, Thomas Schönhoff wrote: > Hello, > > > Prof Brian Ripley said the following on 13.10.2004 10:34: > > On Wed, 13 Oct 2004, Thomas Schönhoff wrote: > > > remove.packages("nlme", "/usr/local/lib/R/site-library") # a duplicate > > foo <- installed.packages(lib="/usr/local/lib/R

Re: [R] "Centered" dummy variables; non zero/one coding

2004-10-13 Thread Prof Brian Ripley
This can done by setting a contrast function or matrix on a variable. Look in e.g. chapter 6 of MASS (the only comprehensive tutorial on coding factors in R, it seems). On Tue, 12 Oct 2004, Peter Holck wrote: > I'm uncertain if this is perhaps a stupid question: > > I want to create "centered"

RE: [R] random forest -optimising mtry

2004-10-13 Thread Liaw, Andy
What I would try is repeat the RF with both mtry a couple of times and see how variable those inter-class proximities are. Even though the differences seem large, they could be within the noise. If so, your only option would be to increase the number of trees. If you are getting 0 OOB error, (an

[R] random forest -optimising mtry

2004-10-13 Thread Ute
Dear R-helpers, I'm working on mass spectra in randomForest/R, and following the recommendations for the case of noisy variables, I don't want to use the default mtry (sqrt of nvariables), but I'm not sure up to which proportion mtry/nvariables it makes sense to increase mtry without "overtunin

[R] importing data

2004-10-13 Thread Andreas Betz
Hi all! I am newcomer to R and I am loborating on this problem: How do I import data from a CD into R for further evaluation. Using code "newgotcha <- read.table("e:\\asciiwin\\gotcha.dat") " returmns a list containing an additional column with indices. The command "plot (gotcha)" gives the

Re: [R] Cleaning up R-2.0.0 installation on GNU/Linux ?[SOLVED]

2004-10-13 Thread Thomas Schönhoff
Hello, Prof Brian Ripley said the following on 13.10.2004 10:34: On Wed, 13 Oct 2004, Thomas Schönhoff wrote: remove.packages("nlme", "/usr/local/lib/R/site-library") # a duplicate foo <- installed.packages(lib="/usr/local/lib/R/site-library")[, 1] install.packages(foo, lib="/usr/local/lib/R/site-

Re: [R] Cleaning up R-2.0.0 installation on GNU/Linux ?

2004-10-13 Thread Prof Brian Ripley
On Wed, 13 Oct 2004, Thomas Schönhoff wrote: > Hello, > > first, big thank you to the developers of R, the release of recent > version is very nice.! > > Well, it seems that I've messed up my current installation by using two > update mechanism at the same time on my Debian Sid box. > > That

[R] "Centered" dummy variables; non zero/one coding

2004-10-13 Thread Peter Holck
I'm uncertain if this is perhaps a stupid question: I want to create "centered" dummy variables to use in a call to glm(), and wondering if there's some slick method in R to do so. That is, rather than have a factor, which results in a glm() fit returning coefficients specifying either absence or

[R] Cleaning up R-2.0.0 installation on GNU/Linux ?

2004-10-13 Thread Thomas Schönhoff
Hello, first, big thank you to the developers of R, the release of recent version is very nice.! Well, it seems that I've messed up my current installation by using two update mechanism at the same time on my Debian Sid box. That is : 1)Debian Way per apt-get freshing up my R installation 2) ir

Re: [R] displaying sample size in boxplots

2004-10-13 Thread Patrick Drechsler
Gabor Grothendieck wrote on 13 Oct 2004 07:16:07 MET: > Warnes, Gregory R pfizer.com> writes: > : > Since I'm still new to R: can somebody give me a pointer to the > : > docs where to find instructions on a package (not a function)? I > : > can find the man pages to specific functions with ? >

Re: [R] R: r course

2004-10-13 Thread Prof Brian Ripley
`second year statistics students' is rather imprecise, but I will guess these are undergraduate statistics majors. However, the book by Nolan & Speed listed in the R FAQ would seem to be at about that level and is rather close to your description. @Book{Nolan.Speed.00, author = {Deborah N

Re: [R] debugging non-visible functions

2004-10-13 Thread Prof Brian Ripley
On Wed, 13 Oct 2004, Murad Nayal wrote: > I would like to step-through a non-visible function. but apparently I > don't know enough about namespaces to get that to work: > > > methods(predict) > ... deleted lines ... > [27] predict.rpart* predict.smooth.spline* > [31] predict.s