Re: [R] Cramer-von Mises test

Thanks Peter...I hadn't some across the sos package before, but I'm sure I'll be putting it to good use from now on! Sean. On Mon, Jul 19, 2010 at 2:55 PM, Peter Ehlers wrote: > On 2010-07-18 22:04, Sean Carmody wrote: >> >> I have found the Kolmogorov-Smirnov (ks.test) and Anderson-Darling >> (

Re: [R] simple loop(?) analysing subsets

Not to hijack the thread, but for my edification, what are the advantages/disadvantages of split() + lapply() compared to by()? Josh On Sun, Jul 18, 2010 at 9:50 PM, Dennis Murphy wrote: > Hi: > > Time to jack up your level of R knowledge, courtesy of the apply family. > > The 'R way' to do what

Re: [R] Cramer-von Mises test

On 2010-07-18 22:04, Sean Carmody wrote: I have found the Kolmogorov-Smirnov (ks.test) and Anderson-Darling (ad.test in package ADGofTest) goodness of fit tests, but was wondering whether there's a package implementing the Cramer-von Mises test for a general distribution (i.e. not just the one fo

Re: [R] simple loop(?) analysing subsets

Hi: Time to jack up your level of R knowledge, courtesy of the apply family. The 'R way' to do what you want is to split the data by species into list components, run lm() on each component and save the resulting lm objects in a list. The next trick is to figure out how to extract what you want,

Re: [R] Question about KLdiv and large datasets

I just answered this but realize that I did so off-list. So, for completeness, here's what I said: I think I see the problem. From ?KLdiv, you're getting the modeltools help page. What you need is the flexmix help page for KLdiv. Just get to the flexmix index page (you can do ?flexmix and then cl

Re: [R] Connecting to MS Access database

Hey Xin this is a piece I copied from my blog. library(RODBC) mdbConnect<-odbcConnectAccess("C:\\temp\\demo.mdb") sqlTables(mdbConnect) demo<-sqlFetch(mdbConnect, "tblDemo") odbcClose(mdbConnect) rm(demo) On Mon, Jul 19, 2010 at 12:05 AM, Xin Ge wrote: > Hi All, > > Can anyone please suggest me

[R] Connecting to MS Access database

Hi All, Can anyone please suggest me from where should I start to learn about 'how to connect to access db' ? How if someone has some written code and I can go over that to understand and make necessary changes... any help would be highly appreciated, -- Xin Ge. [[alternative HTML vers

[R] Cramer-von Mises test

I have found the Kolmogorov-Smirnov (ks.test) and Anderson-Darling (ad.test in package ADGofTest) goodness of fit tests, but was wondering whether there's a package implementing the Cramer-von Mises test for a general distribution (i.e. not just the one for the normal distribution in the nortest pa

Re: [R] Question about KLdiv and large datasets

Is the 'eps' argument part of KLdiv (was not able to find that in the help pages) or part of a general environment (such as the graphics parameters 'par' ) ? I am asking so that I can read about it what it actually does to resolve the question you already raised about its reliability... Ralf On F

Re: [R] specifying column names in a vector of characters and the use?

On Jul 18, 2010, at 10:09 PM, Seth wrote: Hi, What I would like to do is have a data.frame with column names and have these column names stored as strings in another vector. Then I would like to be able to access the data.fram columns via referencing the vector of names. The code belo

Re: [R] specifying column names in a vector of characters and the use?

Try: table.1[[hold[1]]] Cheers, Simon. On 19/07/10 12:09, Seth wrote: Hi, What I would like to do is have a data.frame with column names and have these column names stored as strings in another vector. Then I would like to be able to access the data.fram columns via referencing the vector

Re: [R] specifying column names in a vector of characters and the use?

Hello, What I would like to do is have a data.frame with column names and have these column names stored as strings in another vector. Then I would like to be able to access the data.fram columns via referencing the vector of names. The code below shows the last few executions that failed to r

[R] specifying column names in a vector of characters and the use?

Hi, What I would like to do is have a data.frame with column names and have these column names stored as strings in another vector. Then I would like to be able to access the data.fram columns via referencing the vector of names. The code below shows the last few executions that failed to retri

Re: [R] Help with Reshaping from Wide to Long

Hi Phil and Jeff, Thanks so much for taking the time to help me solve this issue! Both approaches work perfectly. Each of your approaches helped me learn more about what R can do. I really appreciate your help! Very best regards, John -- View this message in context: http://r.789695.n4

Re: [R] re. Mathematica and R

The wiki on rcom.univie.ac.at in section R(D)COM, rcom, and other software systems has an example how to use R from within Mathematica on Windows. You will need to install statconnDCOM and rcom available from the same site. On Jul 17, 2010, at 1:06 PM, Alan Kelly wrote: > David - information on ca

Re: [R] Inserting an image into a PDF file

Hi On 12/07/2010 4:33 a.m., Marc Schwartz wrote: On Jul 11, 2010, at 10:51 AM, David Winsemius wrote: On Jul 11, 2010, at 11:06 AM, Dennis Fisher wrote: Colleagues, I am creating a PDF document from R using pdf(), the document contains text / images created with R. I also have an image

Re: [R] any one knowing how to install library (equate)?

Hi Grace, I'm using the UCLA mirror and this code installs equate fine for me. #install the package (note that the quoted name) install.packages("equate") #load the package library(equate) It also looks like the packages has passed the CRAN checks for all OSes, so availability is probably not a

Re: [R] simple loop(?) analysing subsets

Hello, Take a look at ?by Something like: by(data, data[ , "species id"], function(x) {lm(o2con ~ bm, data = x)}) As an aside, it can be a bit cumbersome to deal with spaces in a variables name (because it has to be quoted then). Also since data is a function, something like mydata might make

Re: [R] any one knowing how to install library (equate)?

On 18/07/2010 2:57 PM, ying_chen wang wrote: I was told to go to packages on the tool bar and select 'install packages'. Then choose library (equate). My R is version 2.11. I tried multiple times with different locations of R within US. Yet, I didn't see library(equate) being offered in the pack

[R] simple loop(?) analysing subsets

Hi All, I have a large data set with many columns of data. One of these columns is a species identifier and the remainder are variables such as temperature or mass. Currently I am carrying out a single regression on subsets of the data set, e.g. separated data sets with only the data from one spe

Re: [R] Neural Network

I'd start with the nnet library type: ?nnet CS - Corey Sparks, PhD Assistant Professor Department of Demography and Organization Studies University of Texas at San Antonio 501 West Durango Blvd Monterey Building 2.270C San Antonio, TX 78207 210-458-3166 corey.sparks 'at' utsa.edu https://row

[R] Neural Network

Hi all, I am working for my master's thesis and I need to do a neural network to forecast stock market price, with also external inputs like technical indicators. I would like to know which function and package of R are more suitable for this study. Thanks a lot for your response, Arnaud TREB

Re: [R] Reading files with varying number of columns

Thanks for the suggestion.  I had read "fill" but didn't pick up on the details.  It only uses the first five columns to determine column dimensions, so it ignored my last two columns.  Specifying "col.names" with fill=TRUE did the trick! Aloha, Tim - Original Mess

[R] Question from day2 beginner

Hello, I just started learning R and have a very basic question: When I try ar model and extract residuals it returns Null. Could someone explain what's going on here? Does that mean the model is null? But it seems residual has a length=# observation of the original series according to the mod

Re: [R] R: package "plotrix"

On Jul 18, 2010, at 1:32 PM, wrote: -Messaggio originale- Da: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Inviato: dom 18/07/2010 18.58 A: mau...@alice.it Cc: j...@bitwrit.com.au; r-h...@stat.math.ethz.ch Oggetto: Re: [R] package "plotrix" On 18.07.2010 17:59, mau...@alice.

Re: [R] Continuing on with a loop when there's a failure

Thanks very much again David for your helpful answers. However, the code STILL does not appear to be working properly! Even though the third time through the loop *should* work, it appears that R has given up after the second time through the loop. What I mean is: although y2 causes the lrm fu

[R] any one knowing how to install library (equate)?

I was told to go to packages on the tool bar and select 'install packages'. Then choose library (equate). My R is version 2.11. I tried multiple times with different locations of R within US. Yet, I didn't see library(equate) being offered in the packages. Thanks. Grace [[alternative HT

Re: [R] loop troubles

On Jul 18, 2010, at 4:21 PM, Joe P King wrote: This is the latest code I have been trying, but it when I try to turn the vectors of SD into a vector of variances, it turns into a scalar. n=c(NULL,13,89,50) ybar=c(NULL,1.58,1.26,0.80) sd=c(NULL,2.19,4.18,5.47) mu=1;sigma=1;p0=0;var1=NULL;p=NU

Re: [R] loop troubles

This is the latest code I have been trying, but it when I try to turn the vectors of SD into a vector of variances, it turns into a scalar. n=c(NULL,13,89,50) ybar=c(NULL,1.58,1.26,0.80) sd=c(NULL,2.19,4.18,5.47) mu=1;sigma=1;p0=0;var1=NULL;p=NULL;pn=NULL for (i in 2:length(n)){ var1[i] = sd[i]^

Re: [R] Import of specific column of many space-delimited text files

That worked perfectly well. Thanks a lot! -- View this message in context: http://r.789695.n4.nabble.com/Import-of-specific-column-of-many-space-delimited-text-files-tp2293273p2293342.html Sent from the R help mailing list archive at Nabble.com. __ R

Re: [R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

On Jul 18, 2010, at 2:52 PM, Hadley Wickham wrote: The problem is in data.frame[ and any NA in a logical vector will return a row of NA's. This can be avoid by wrapping which() around the logical vector which seems entirely wasteful or using subset(). The basic philosophy that causes this

Re: [R] Import of specific column of many space-delimited text files

Hi, There might be better ways, but here's something that might help you get started. Obviously you'll need to tweak a lot of things. If your files do not have a uniform number of columns and the column you want is not always in the same position, it may be easier to just read it all in and then

Re: [R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

> The problem is in data.frame[ and any NA in a logical vector will return a > row of NA's. This can be avoid by wrapping which() around the logical vector > which seems entirely wasteful or using subset(). The basic philosophy that causes this behaviour is sensible in my opinion: missing values m

[R] Import of specific column of many space-delimited text files

Hi, I have about 300 space-delimited text files and from each file I want to import one specific column into R to create a data frame where all imported columns are included. Is there a smart way to do so? Thanks! -- View this message in context: http://r.789695.n4.nabble.com/Import-of-specifi

Re: [R] R: package "plotrix"

On 18.07.2010 19:57, David Winsemius wrote: On Jul 18, 2010, at 1:32 PM, wrote: -Messaggio originale- Da: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Inviato: dom 18/07/2010 18.58 A: mau...@alice.it Cc: j...@bitwrit.com.au; r-h...@stat.math.ethz.ch Oggetto: Re: [R] package

[R] R: package "plotrix"

-Messaggio originale- Da: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Inviato: dom 18/07/2010 18.58 A: mau...@alice.it Cc: j...@bitwrit.com.au; r-h...@stat.math.ethz.ch Oggetto: Re: [R] package "plotrix" On 18.07.2010 17:59, mau...@alice.it wrote: > I installed package plotrix

Re: [R] log transformation

jarry hotmail.com> writes: > > > Hi, > > I am reading an article in which, first, some given "p-results" are > obtained, and, second, these are afterwards transformed and expressed as the > "NEGATIVE base-10 logarithm of the p-value". > > My question is if anyone could indicate how is such t

Re: [R] latex table question

I' ve seen latex questions being solved here that's why I sent my question to this list but someone already told me about the latex forum so the same question has already been solved there, Sorry about that.   - Original Message > From: Uwe Ligges > To: Felipe Carrillo > Cc: r-h...@

Re: [R] Correct function name for text display of S4 object

On 13.07.2010 03:35, Andrew Liu wrote: Hello, I am working on an R package for storing order book data. I currently have a display method that has the following output (ob is an S4 object): display(ob) Current time is 09:35:02 Price Ask Size -- 11.42

Re: [R] latex table question

Why do you think R-help is a mailing list about LaTeX? Best, Uwe Ligges On 13.07.2010 23:05, Felipe Carrillo wrote: Hi: My head is spinning with this latex doc so hopefully after I align my tables to the left of the page my headache are going to be over. I always use: \hspace*{-0.1in} to move

Re: [R] Error using the mi package

Since we do not have the data and cannot reproduce: What about sending a reproducible example to the "mi" maintainer in the first step? Best, Uwe Ligges On 15.07.2010 16:48, Andrew Miles wrote: I'm trying to impute data using the mi package, but after running through almost the entire first

Re: [R] Toggle between the various pages for multi-page figures

On 16.07.2010 18:38, mahesh samtani wrote: Hello, I am a new R user having transitioned over from S-plus recently. I have a question that is probably very trivial but I am having trouble finding a solution. In S-plus, graphic pages are created as tabs when multi-page figures are created. I hav

Re: [R] Continuing on with a loop when there's a failure

On Jul 18, 2010, at 12:28 PM, Josh B wrote: Thanks very much again David for your helpful answers. However, the code STILL does not appear to be working properly! Even though the third time through the loop *should* work, it appears that R has given up after the second time through the loo

Re: [R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

On Jul 18, 2010, at 12:02 PM, stephen sefick wrote: I am confused by the behavior of the below piece of code. The NAs are making it past the logical call ==0. I am sure that I am missing something. I just don't understand this behavior. Thanks for your help in advance. code###

Re: [R] Eliminating case numbers in a dendrogram

On 18.07.2010 19:02, Holger Steinmetz wrote: Hi folks, I conducted a hierarchical cluster analysis. As I wanted to illustrate the result, I created a dendrogram with the code plot(as.dendrogram(fit),sub="",xlab="",ylab="Heterogeneity",nodePar = list(lab.cex=.5,pch=NA,xlab="")) However, the

Re: [R] Reading files with varying number of columns

See argument "fill" in ?read.table Uwe Ligges On 18.07.2010 12:14, Tim Clark wrote: Dear R list, I am trying to read files with a varying number of columns and can't figure out how to get them in the proper format. Some rows have five value and some have seven. Is there a way to read them in

Re: [R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

On 18.07.2010 18:02, stephen sefick wrote: I am confused by the behavior of the below piece of code. The NAs are making it past the logical call ==0. I am sure that I am missing something. I just don't understand this behavior. Thanks for your help in advance. code

[R] Eliminating case numbers in a dendrogram

Hi folks, I conducted a hierarchical cluster analysis. As I wanted to illustrate the result, I created a dendrogram with the code plot(as.dendrogram(fit),sub="",xlab="",ylab="Heterogeneity",nodePar = list(lab.cex=.5,pch=NA,xlab="")) However, the dendrogram contains the case numbers and, as I ha

Re: [R] package "plotrix"

On 18.07.2010 17:59, mau...@alice.it wrote: I installed package plotrix because reading its vignette it looks like it can help me solve a "legend" problem. The package instaleed correctly on my Mac OS/X 10.5.8 But I cannot reproduce the examples centered on function "lgendg". You mean "lege

[R] Reading files with varying number of columns

Dear R list, I am trying to read files with a varying number of columns and can't figure out how to get them in the proper format.  Some rows have five value and some have seven.  Is there a way to read them in so that two empty columns are added to every row that has only five values?  An examp

[R] log transformation

Hi, I am reading an article in which, first, some given "p-results" are obtained, and, second, these are afterwards transformed and expressed as the "NEGATIVE base-10 logarithm of the p-value". My question is if anyone could indicate how is such transformation achieved with R and how are they

[R] NA preserved in logical call - I don't understand this behavior because NA is not equal to 0

I am confused by the behavior of the below piece of code. The NAs are making it past the logical call ==0. I am sure that I am missing something. I just don't understand this behavior. Thanks for your help in advance. code### left <-

[R] package "plotrix"

I installed package plotrix because reading its vignette it looks like it can help me solve a "legend" problem. The package instaleed correctly on my Mac OS/X 10.5.8 But I cannot reproduce the examples centered on function "lgendg". > library(plotrix) > plot(0.5,0.5,xlim=c(0,1),ylim=c(0,1),type="n

[R] heatmap.2 - change column & row locations; angle / rotate

Esteemed R user's, I'm struggling to achieve some details of a heatmap using heatmap.2(): 1. Change label locations, for both rows & columns from the default right & bottom, to left and top. Can this be done within heatmap.2()? Or do i need to suppress this default behavior (how) and call a ne

Re: [R] re. Mathematica and R

The wiki on rcom.univie.ac.at in section R(D)COM, rcom, and other software systems has an example how to use R from within Mathematica on Windows. You will need to install statconnDCOM and rcom available from the same site. On Jul 17, 2010, at 1:06 PM, Alan Kelly wrote: > David - information on

Re: [R] Identify points (was Plot error)

You can always use the identify() function which will leave the id of the point on the plot; see ?identify. But you do eventually have to sort out your Mac problem. Sorry, I can't help with that. -Peter Ehlers On 2010-07-18 7:43, James Platt wrote: This is exactly what I want as I will have

[R] producing biplot of canonical correlation analysis

Dear List, I would like to obtain biplot containing arrows and variable labels. However, after checking the previous message and package CCA and function anacor, my attempt remains an open question. Please kindly advise code or package I might miss to achieve it. Thank you. Elaine [[al

Re: [R] Identify points (was Plot error)

On Jul 18, 2010, at 9:43 AM, James Platt wrote: This is exactly what I want as I will have several thousand data points on the final graph i make, so the scroll over option is ideal. I've read the TeachingDemos pdf, I'm working on a Mac so Cannot use HWidentify. Not true. I have install

[R] Giovanna Jonalasinio è fuori ufficio, Out of Office

Risposta automatica dal 18/7/10 fino al 27/7/10 I'm not in Rome and I won't read my email regularly untill the 27th of July Non sono a Roma e leggerò la posta sporadicamente fino al 27 luglio [[alternative HTML version deleted]] __ R-help@r-p

Re: [R] Help with a problem

You can also use 'embed' to create a list of indices you can use to do the test: > dat dsc1 c2 1 2010-04-03 100 0 2 2010-04-30 11141 15 3 2010-05-01 3 16 4 2010-05-02 7615 14 5 2010-05-03 6910 17 6 2010-05-04 5035 3 7 2010-05-05 3007 15 8 2010-05-06 4 14 9 2

Re: [R] Identify points (was Plot error)

This is exactly what I want as I will have several thousand data points on the final graph i make, so the scroll over option is ideal. I've read the TeachingDemos pdf, I'm working on a Mac so Cannot use HWidentify. I have installed and loaded the TeachingDemos package but im having trouble

Re: [R] loop troubles

On Jul 18, 2010, at 9:12 AM, Joe P King wrote: I tried that, this is what I tried, but I only get it to do one iteration and then it wont cycle back. I am not sure how to tell it how to look at the previous precision to add to the current new sample precision to get the new precision. A

Re: [R] slplot issue

It's truly amazing how many people asking for help can't be bothered to report what add-on packages they're using. Shawn's query is a particularly egregious case: slplot is mentioned in pkg:ccmn, which can be found on CRAN, but which depends on pkg:pcaMethods which is on Bioconductor. And pcaMetho

Re: [R] loop troubles

I tried that, this is what I tried, but I only get it to do one iteration and then it wont cycle back. I am not sure how to tell it how to look at the previous precision to add to the current new sample precision to get the new precision. This is my first try at a loop so I am not sure if I am doin

Re: [R] loop troubles

On Jul 17, 2010, at 9:09 PM, Joe P King wrote: Hi all, I appreciate the help this list has given me before. I have a question which has been perplexing me. I have been working on doing a Bayesian calculating inserting studies sequentially after using a non-informative prior to get a meta-analys

Re: [R] sort file names in numerical order

Another option: require(gtools) ?mixedsort > mixedsort(fileNames) [1] "A1" "A2" "A10" "B1" "B2" "B10" -- David On Jul 18, 2010, at 5:16 AM, Duncan Mackay wrote: Hi Yes it is possible- one way is: fileNames[order(sprintf("%02s", sub("[[:upper:]]","", fileNames)))] [1] "A1" "B1" "A

Re: [R] Continuing on with a loop when there's a failure

On Jul 18, 2010, at 4:25 AM, Josh B wrote: Hello Peter, I tried your suggestion, but I was still not able to get it to work. Would you mind looking at my code again? Here's what I'm trying: x <- read.table(textConnection("y1 y2 y3 x1 x2 indv.1 bagels donuts bagels 4 6 indv.2 donuts donuts

Re: [R] R graphic help

On 07/18/2010 06:58 PM, Timothy O'Brien wrote: Dear All, I've done some searching, but to no avail -- I'm plotting x-y data via the "plot" command, and the "log=x" command. The graphed x values are in scientific notation (1e-02 1e-01 1e+00 etc). Might you have some idea on how I can get the pl

Re: [R] R graphic help

Hi, There may be a simpler way but try this, plot(10^jitter(seq(-2,4,length=10)), 1:10, log="x", xaxt="n") axis(1, at = axTicks(1),labels = format(axTicks(1),scientific=FALSE)) HTH, baptiste On 18 July 2010 10:58, Timothy O'Brien wrote: > Dear All, > > I've done some searching, but to no avai

[R] R graphic help

Dear All, I've done some searching, but to no avail -- I'm plotting x-y data via the "plot" command, and the "log=x" command. The graphed x values are in scientific notation (1e-02 1e-01 1e+00 etc). Might you have some idea on how I can get the plot to uses the x values (0.01 0.10 1.0 10.0 etc)

Re: [R] Namibia becoming NA

Berwin A Turlach wrote: > Did you think of trying some variations of "na.strings"? ;-) > > IMO, the simplest way of coding missing values in CSV files is to have > two consecutive commas; not some code (whether NA, 99, 999, -1, ...) > between them. Yes. Arguably, na.strings=NULL should be the d

Re: [R] sort file names in numerical order

Hi Yes it is possible- one way is: fileNames[order(sprintf("%02s", sub("[[:upper:]]","", fileNames)))] [1] "A1" "B1" "A2" "B2" "A10" "B10" Regards Duncan Duncan Mackay Department of Agronomy and Soil Science University of New England ARMIDALE NSW 2351 Email home: mac...@northnet.com.au A

Re: [R] Continuing on with a loop when there's a failure

Oops, forgot to add the line: library(Design) #the lrm function is in the Design library ...but even when I load the Design library, the loop still doesn't work. It stops after failing on the second run of the loop: > results y1 y2 y3 [1,] 0.6976063 NA NA (The third run through th

Re: [R] Continuing on with a loop when there's a failure

Hello Peter, I tried your suggestion, but I was still not able to get it to work. Would you mind looking at my code again? Here's what I'm trying: x <- read.table(textConnection("y1 y2 y3 x1 x2 indv.1 bagels donuts bagels 4 6 indv.2 donuts donuts donuts 5 1 indv.3 donuts donuts donuts 1 10 indv.

Re: [R] Namibia becoming NA

G'day Ted, On Sun, 18 Jul 2010 09:25:09 +0100 (BST) (Ted Harding) wrote: > On 18-Jul-10 05:47:03, Suresh Singh wrote: > > I have a data file in which one of the columns is country code and > > NA is the > > code for Namibia. > > When I read the data file using read.csv, NA for Namibia is being >

Re: [R] Namibia becoming NA

On 18-Jul-10 05:47:03, Suresh Singh wrote: > I have a data file in which one of the columns is country code and NA > is the > code for Namibia. > When I read the data file using read.csv, NA for Namibia is being > treated as > null or "NA" > > How can I prevent this from happening? > > I tried th

Re: [R] Help with a problem

Hi all, zoo::rollmean() is a nice idea. But if I understand Mike correctly, he wants 5 out of any 7 consecutive logicals to be TRUE, where these 5 do not necessarily need to be consecutive themselves. (remaining open question: could, e.g., the condition on c1 be TRUE for rows 1,2,3,4,5 and on