[R] Plot link points

Hi. This may be a straight-forward question... I have a time series from 1950 to 2000 (by year) with missing values. Sometimes only a year has value and the adjacent years have missing data, e.g. 1950 2.7 1951 1952 1953 3.4 1954 1955 1956 2.9 1957 3.6 1958 2.7 When I use plot (type='l'), the

Re: [R] Plot link points

Kitty Lee wrote: Hi. This may be a straight-forward question... I have a time series from 1950 to 2000 (by year) with missing values. Sometimes only a year has value and the adjacent years have missing data, e.g. 1950 2.7 1951 1952 1953 3.4 1954 1955 1956 2.9 1957 3.6 1958 2.7 When I use

[R] Clinical significance estimation

Hi all, On a request from a colleague, I programmed an implementation of the clinical significance estimation suggested by Jacobson and Truax some years ago, with a few generalizations and a couple of plot functions. I have wrapped this up in a package, but I am reluctant to further clutter CR

Re: [R] Plot link points

try to use xyplot of latice package. If you have a data.frame called data, whit coluns called year and x, you could do something like: require(lattice) xyplot(x~year,data, type = "a", auto.key = list(space = "right", points = FALSE, lines = TRUE)) __

[R] Help needed for Loading "tm" package

Howdy Gurus again Thanks to Tony.Breyal, I was able to writing the following script for analyzing a text document. But I got an error with "tm' package. I don't why I got the error from the R script below. I think I followed proccess of R tm manual. I use R v2.8.1. and tm_0.3-3.zip under Win XP.

Re: [R] Help needed for Loading "tm" package

On Sat, 10 Jan 2009, Kum-Hoe Hwang wrote: Howdy Gurus again Thanks to Tony.Breyal, I was able to writing the following script for analyzing a text document. But I got an error with "tm' package. I don't why I got the error from the R script below. I think I followed proccess of R tm manual.

Re: [R] rpart with interval censored data crashes R

On a Leopard Mac with the Urbanek compiled 64 bit R, one sees this: > library(rpart) > library(survival) Loading required package: splines > fit<-rpart(Surv(N,Y,type="interval2")~Salt+pH+Temp, data=myD) *** caught segfault *** address 0x0, cause 'memory not mapped' Traceback: 1: .C(C_rpartexp

[R] AT&T Researchers and the New York Times

I've submitted a posting to the NY Times blog mentioned by David Smith, specifically mentioning John, Rick and Allan. (Now awaiting moderation.) Of course there are many others deserving credit but far to numerous to list so I didn't get started. On Thu, Jan 8, 2009 at 5:58 PM, Robert Wilkins

Re: [R] The R Inferno

Extremely nice. Tons of useful information presented in a very entertaining way. Thanks! Tom Patrick Burns wrote: "The R Inferno" is now on the Burns Statistics website at http://www.burns-stat.com/pages/Tutor/R_inferno.pdf Abstract: If you are using R and you think you're in hell, this is

Re: [R] The R Inferno

Excellent read, Patrick. A very useful and clear guide. On 2009.01.09 16:14:49, Patrick Burns wrote: > "The R Inferno" is now on the Burns Statistics website at > > http://www.burns-stat.com/pages/Tutor/R_inferno.pdf > > Abstract: If you are using R and you think you're in hell, > this is a map

Re: [R] Ashlee Vance's article on R in the New York Times

On Thu, Jan 8, 2009 at 1:26 AM, Robert Wilkins wrote: > ...The user interface for R, otherwise known as the S programming language > has the same origins as C and Unix We could take this one step further, and note that C's design (its "user interface"?) was based on BCPL, which was developed

[R] Print specific matrix value???

Hello All, I'm trying to print specific row and column for Observed_Scores matrix, however, when I execute the command "Observed_Scores[1,1]", I get the message "Error in Observed_Scores[1, 1] : incorrect number of dimensions". Could you please help and let me know where is the mistake? Here is

[R] Fitting pareto to some data

Dear R-users, I am trying to fit pareto distribution to some data but i've one problem. Using optim to calculate the maximum of the likelihood function of the pareto I use as start parameters the moments method(using the distribution function in the package actuar): media=mean(x) var=mean(x^2)-me

[R] Print specific matrix value???

Hello All, I'm trying to print specific row and column for Observed_Scores matrix, however, when I execute the command "Observed_Scores[1,1]", I get the message "Error in Observed_Scores[1, 1] : incorrect number of dimensions". Could you please help and let me know where is the mistake? Here is

Re: [R] The R Inferno

Hi , I have been trying to convince Pat to come up with a book on this. He can add in the chapters on *Purgatorio* (Purgatory), and *Paradiso* (Paradise ) Regards, Ajay On Sat, Jan 10, 2009 at 9:33 PM, Jason Morgan w

Re: [R] Plot link points

Look at: http://geography.uoregon.edu/geogr/topics/maps.htm -Roy On Jan 10, 2009, at 1:48 AM, Jim Lemon wrote: Kitty Lee wrote: Hi. This may be a straight-forward question... I have a time series from 1950 to 2000 (by year) with missing values. Sometimes only a year has value and the adjace

Re: [R] Print specific matrix value???

Hi Nidhi, in your very last line, you set Observed_Scores = GenData[1] After that, Observed_Scores is not a matrix any more, but a list with a single entry Observed_Scores$model, which is in fact a data.frame, but with different dimensions than you initially set Observed_Scores to. HTH, Step

Re: [R] Fw: Re: R Stacked Histogram

Hi Jason, You can do: qplot(mpg, wt, data = mtcars, colour = factor(cyl)) + opts(legend.position = c(0.5, 0.5)) Hadley On Thu, Jan 8, 2009 at 3:04 PM, Jason Rupert wrote: > Thank you again for your response. > > This worked great. > > Quick question about the legend for qplot. Instead of bein

Re: [R] R in the NY Times

“We have customers who build engines for aircraft. I am happy they are not using freeware when I get on a jet.” The lady who made this comment, Anne H. Milley, director of technology product marketing at SAS, has written a response to try and clarify what she meant (funilly enough, i got this link

[R] R version requires Bioconductor 2.3 version

Dear support, I am new to R. I installed R 2.8.1 under Windows, and upon starting it I get the message that it requires Bioconductor version 2.3. I thought this was kind of installed in parallel or so. Where can I upgrade it ? Best regards, Andreas -- Andreas Scherer, PhD,APM CEO Spheromics

[R] Arguments for Rcmd BATCH

Sorry to bother you and I hope this will be easy to solve. I am trying to run R scripts in batch, called from another programme, under Windows XP. I have R 2.4.1 I need to be able to pass a list of numeric arguments to the script as well. The Rscript and its input data file are attached

[R] Applying 'lm' in each group

Hi,  I want to do regression in each group. I made the group the following way. Now I want to run regression in each of the three groups.. Suppose , x<-c(0.5578196,6.5411662,13.2728619,2.0217271,6.7216176,3.37220617,2.5773252,7.2600583,15.3731026,3.4140288,8.1335874,     15..6476637,4.3014084,

[R] Calling C from FORTRAN and vice versa

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Dear R-Help, as I am not sure where to put this suggestion for improving the documentation in "Writing R extensions" Section 6.6 (Calling C from FORTRAN and vice versa) I send it to R-Help as suggested by the R Bug Report page. In Sect 6.6 an example

[R] JGR editor setting problem

Hi all, I am using the JGR (Java GUI for R ) and I am facing some problem. I want to use the help agent in the editor. I try to change the setting in the preference to use the help agent in editor, but it does not work. Though I have apply and save this seting ,it will return to the default setin

Re: [R] Arguments for Rcmd BATCH

Your R is far too old, and in particular too old for the example you are trying. Please do you as asked in the posting guide to do before posting, and upgrade. Then you will have Rscript, as described in the current 'An Introduction to R' manual for this purpose (with examples). On Sat, 10

[R] Setting a contingency table

Hi, I want to set a make a contingency table which will look like this.. The problem is that I can't  set the table like the following. col1 col2   Total     row1  a b n10 rp1   rp2    100    

Re: [R] R in the NY Times

On Sat, Jan 10, 2009 at 8:11 AM, Tony Breyal wrote: > "We have customers who build engines for aircraft. I am happy they are > not using freeware when I get on a jet." > > The lady who made this comment, Anne H. Milley, director of technology > product marketing at SAS, has written a response to t

Re: [R] R in the NY Times

2009/1/10 Tony Breyal : > [SAS marketroid quote] > "In fact, SAS values open-source software." But clearly not enough to open-source SAS itself. It would seem that SAS values _other_people's_ open source. If SAS was open source and free, then SAS would collect on all the other things "Customer

Re: [R] R in the NY Times

more on the reasons R is bad for you http://www.decisionstats.com/2009/01/top-ten-rrreasons-r-is-bad-for-you/ On Sun, Jan 11, 2009 at 12:01 AM, Barry Rowlingson < b.rowling...@lancaster.ac.uk> wrote: > 2009/1/10 Tony Breyal : > > > [SAS marketroid quote] > > "In fact, SAS values open-source soft

Re: [R] R version requires Bioconductor 2.3 version

Andreas Scherer wrote: Dear support, This is not "support", but the R-help mailing list. Please do read docuemntation on how to upgrade R and contributed packages such as those from the Bioconductor project. I am new to R. I installed R 2.8.1 under Windows, and upon starting it I get the

Re: [R] Print specific matrix value???

Nidhi Kohli wrote: Hello All, I'm trying to print specific row and column for Observed_Scores matrix, however, when I execute the command "Observed_Scores[1,1]", I get the message "Error in Observed_Scores[1, 1] : incorrect number of dimensions". Could you please help and let me know where i

Re: [R] bug in R2WinBUGS

John Smith wrote: In newest version of R2WinBUGS, the default directory is changed to working.directory, but never changed back once finished bugs call. Thank you for the report! I will see what happened and release a fixed version tomorrow. Best wishes, Uwe [[alternative HTML

Re: [R] R in the NY Times

I think the substance of the issue is that the more eyes on code, the fewer the bugs (assuming a well-designed examination and debugging process is in place, as is typical for large open source projects like R). By this (obvious?)criterion, both the remarks about the dangers of proprietary code an

Re: [R] Applying 'lm' in each group

Bhargab Chattopadhyay yahoo.com> writes: > > Hi, > >  I want to do regression in each group. I made the group the following way. Now I want to run regression > in each of the three groups.. > Suppose , > w > > Now if I write lm(yg~xg) then it won't work. Please advice how to proceed. > > Tha

Re: [R] Applying 'lm' in each group

Did you read the error message? Generally saying "doesn't work" will elicit undesirable responses. Now would be a good time to read the Posting Guide: http://www.R-project.org/posting-guide.html Including: "Say exactly what happened, including any error messages you received. " > lm(yg~xg) E

Re: [R] Applying 'lm' in each group

On 1/10/2009 12:55 PM, Bhargab Chattopadhyay wrote: > Hi, > > > I want to do regression in each group. I made the group the following way. > Now I want to run regression in each of the three groups.. > Suppose , > x<-c(0.5578196,6.5411662,13.2728619,2.0217271,6.7216176,3.37220617,2.5773252,7.2

Re: [R] R in the NY Times

The reactions to the NYT article have certainly made for some interesting reading. Here are some of the links: http://overdetermined.net/site/content/new-york-times-article-r http://jackman.stanford.edu/blog/?p=1053 http://ggorjan.blogspot.com/2009/01/new-york-times-on-r.html several posts on

[R] install.views()

Dear Rxperts, Using R 2.8.1 and trying install.views("Cluster") getting error Error: could not find function "install.views" Please help:-( -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted]] __

Re: [R] install.views()

did you install and load the 'ctv' package? install.packages('ctv') library('ctv') then try your code... Kingsford Jones On Sat, Jan 10, 2009 at 1:00 PM, oscar linares wrote: > Dear Rxperts, > > Using R 2.8.1 and trying > > install.views("Cluster") > > getting error > > Error: could not fin

Re: [R] Setting a contingency table

library(gmodels) ?CrossTable# by Marc Schwartz a<-67; b<-10; c<-79; d<-67; > dft<- with(as.data.frame.table(as.table(c(a=77,b=10,c=79,d=67))), data.frame(Var1 = rep(Var1,Freq)) ) There are probably more efficient ways to replicate entries. I am just drawing a blank right now.

Re: [R] install.views()

Does your session info show that you have ctv installed? > sessionInfo() R version 2.8.0 Patched (2008-11-14 r46932) i386-apple-darwin9.5.0 locale: en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base o

[R] Hmisc-xtable label

Dear all: Does anybody know about label conflicts between xtable and Hmisc? I found a couple of e-mails similar to this problem but is not clear to me how to get around the label problem. The first table(longtable below) is generated with the latex function from Hmisc but for some reason when I

Re: [R] R in the NY Times

Bert Gunter [Sat, Jan 10, 2009 at 08:31:03PM CET]: [...] > Perhaps an argument is that certain code might not get written at all if it > were not proprietary. Device drivers might be an example. Device drivers are not an example. Linux is ubiquitous _despite_ device manufacturers being secretiv

Re: [R] R in the NY Times

There do exist device manufacturers who GPL their device drivers, e.g. http://freshmeat.net/projects/wanpipe/?branch_id=73783&release_id=290741 On Sat, Jan 10, 2009 at 2:31 PM, Bert Gunter wrote: > > I think the substance of the issue is that the more eyes on code, the fewer > the bugs (assuming

Re: [R] Hmisc-xtable label (Solved)

As soon as I sent the e-mail I found a solution to the hyperlink problem. I removed the squared brackets and the number 1 from the line below. I don't know why I had the number there in the first place. \pdfbookmark[1]{Contents}{table} --- On Sat, 1/10/09, Felipe Carrillo wrote: > From: Felipe

[R] how to get the signed distance in SVM?

Dear all, In the svm() function in the package e1071, is there anyway to get the signed distance of a sample point to the separating hyperplane in the feature space? Namely, if the separating hyperplane is given by f(x) = h(x)^T * w - rho, is there any way to get f(x)? Also, in the returned

[R] Problem with compiling shared C/C++ library for loading into R (Linux)

I am using the .Call interface to call c++ code from R. For that, I am trying to create a dynamic library (mylib.so) using "R CMD SHLIB" by linking my own c++ code and an external c++ library (blitz++). The makefile works fine on my Mac, produces mylib.so and I am able to call .Call() from R, but

[R] Problem with compiling shared C/C++ library for loading into R (Linux)

Dear all, Sorry about posting this to R-help by mistake. I just realized it belongs to the r-devel list. So I am reposting it there. - I am using the .Call interface to call c++ code from R. For that, I am trying to create a dynamic library (mylib.so) using "R CMD SHLIB" by linking my o

[R] Rserve/RandomForest does not work with a CSV?

Hi all, We're using Rserve and RandomForest to do classification from within a Java program. The total is about 4 lines of R code: library('randomForest') x y future fit<-randomForest(x,y,no.action=na.roughfix,importance=T,proximity=T) p<-predict(fit, future) What is very frustrating is that we

Re: [R] The R Inferno

Hi Pat & R-list, Congratulation for this very informative and funny reading! I enjoyed so much. Best wishes, miltinho astronauta brazil On Sat, Jan 10, 2009 at 1:31 PM, Ajay ohri wrote: > Hi , > I have been trying to convince Pat to come up with a book on this. > > He can add in the chapters o

Re: [R] Arguments for Rcmd BATCH

I suggested you use Rscript, please do. R CMD BATCH is the Unix-alike notation. Rcmd BATCH is the recommended way on Windows, but R CMD BATCH is also accepted. There is nothing to reproduce here, so here is a simple example. tystie% cat foo.R args <- commandArgs(TRUE) print(args) tystie% Rs

Re: [R] Zipf fitting using R

Hi there, I haven't seen any answers to your question on the list, so here's a rather late response. I'm not familiar with the VGAM package, but if you just want to fit Zipf's law (did you mean Zipf's law or an actual probability distribution?) to some data, why not do this directly with nl

Re: [R] Arguments for Rcmd BATCH

Dear Professor Ripley, Firstly, thank you very much for your advice. I was not aware my R was so old. Sorry I must have missed the instruction on reading the posting guidelines. I have updated it and it now runs again. This is a big step forward. However it is not running properly (or I am doing

[R] temporal join

I have dataframe a: sym date val1 === foo 20090101 a1 foo 20090102 a2 foo 20090103 a3 and dataframe b: sym date val2 === foo 20090104 b1 I would like to join/merge them to generate the following: sym date val2 val1 === foo 20090104 b1 a3 i.e. an equijoin on

[R] problems installing package XML to a computer without an internet connection

Hello, I am hoping for some advice regarding how I can install the XML package which I require to run package tm. Normally I would use the install package option, however, I have to install the packages to a laptop running XP. The laptop does not have an internet connection. Firstly I trie

Re: [R] problems installing package XML to a computer without an internet connection

On a computer with an internet connection use your browser to find and download the XML package zip file. (Note that this step does not involve using R.) Transfer zip file to XP. Start up R on the XP & use the menus: Packages | Install package(s) from local zip file to install it. On Sat, Ja

[R] I'm looking for a book about spatial point patterns (Diggle, 2003)

To understand some functions about spatial point patterns in "spatstat" ,I should know some background about it, and the best way is to read the monograph, and "Statistical Analysis of Spatial Point Patterns" (2nd edt.) is a better choise. But I can not find it anywhere I can. Who can help me? Th

Re: [R] install.views()

Thanks! Your suggestion did it:-) OAL On Sat, Jan 10, 2009 at 3:00 PM, oscar linares wrote: > > Dear Rxperts, > > Using R 2.8.1 and trying > > install.views("Cluster") > > getting error > > Error: could not find function "install.views" > > > Please help:-( > > > -- > Oscar > Oscar A. Linares >

[R] Convert date/time string to date

I am new to R and am trying to import a text file that contains date/time and various fields. I want to sort this matrix by date and also perform calculations on the date field (difference between two dates etc). The format of the string looks as follows: x<-c("25/2/2003 0:00:00") I tried: as.

Re: [R] problems installing package XML to a computer without an internet connection

Hello Gabor, I believe I used the method you described to download zip files for rJava, tm, matrix, filehash etc etc. There is no windows zip file for XML, only a gz file which I can't get to load through install local zip file. The XML zip file I located on http://www.stats.ox.ac.uk/pub/RWin

[R] How to get solution of following polynomial?

Hi, I want find all roots for the following polynomial : a <- c(-0.07, 0.17); b <- c(1, -4); cc <- matrix(c(0.24, 0.00, -0.08, -0.31), 2); d <- matrix(c(0, 0, -0.13, -0.37), 2); e <- matrix(c(0.2, 0, -0.06, -0.34), 2) A1 <- diag(2) + a %*% t(b) + cc; A2 <- -cc + d; A3 <- -d + e; A4 <- -e fn <- fu