Re: [R] Urgent help requested to modify a script

On 10/16/2009 11:01 PM, one2luv wrote: I am hoping someone will tak up this chalenge (I am new to R) I have inheritied an R script but need to change it. The script currently includes hardcoded file locations on lines 12,166 and 167. I need to modify this script to allow the folder to be passed

Re: [R] How to right-align labels in dotchart

Sean Carmody wrote: Thanks very much John, that works. I had tried using las=1 inside the dotchart function itself, but it seems to be one of those occasions where the par parameter is over-ridden by the main plotting function. You can see the results here: http://www.stubbornmule.net/2009/10/a

Re: [R] Recommendation on a probability textbook (conditional probability)

There are many examples in the book. Since I'm refreshing my memory. Is there a more concise one? On Fri, Oct 16, 2009 at 8:26 PM, Ista Zahn wrote: > I like Grinstead and Snell, not least because it's free: > http://www.dartmouth.edu/~chance/teaching_aids/books_articles/probability_book/book.html

Re: [R] Different way of scaling data

Read the whole help page for scale, then try: Newx <- scale(x, min(x), max(x)-min(x) ) -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- >

Re: [R] Generating a stochastic matrix with a specified second dominant eigenvalue

A further idea: Consider the triangular square matrices of the form with decreasing eigenvalues on the diagonal: 1 0 0 0 0 0 .3 .7 0 0 0 0 .4 .2 .4 0 0 0 .2 .4 .2 .2 0 0 .1 .3 .4 .2 .1 0 .2 .2 .2 .2 .1 .1 This would have the specified eigenvalue compo

Re: [R] length error using tapply()

George Aldridge wrote: > > I am new to R and trying to aggregate a character variable. I have > searched, and tapply() seems to hold the best hope, but I am having > trouble with it. I have a character variable 'hab' and 3 factors all of > the same length: > >> length(hab) > [1] 105219 >> leng

Re: [R] Beveridge-Nelson decomposition

sdlywjl666 wrote: > > Dear all, > Is there some function in R to carry out Beveridge-Nelson > decomposition? > Is there some function in R to carry out seasonal adjustments of time > series? > Is there some function in R to carry out x11 or x12 seasonal adjustments ? > > RSiteSearch("Bev

[R] avoiding loops in equation

To illustrate my problem, I have a complex code of several matrices and a vector. To simplify I only used two matrices and a vector as an example: pex<-function(acc, pga, std){ (acc-pga)/std } acc<-seq(.1,1,.1) pga<-matrix(rnorm(9,4,.1),3,3) std<-matrix(rnorm(9,4,.5),3,3) I tried calculating

Re: [R] Recommendation on a probability textbook (conditional probability)

On Fri, Oct 16, 2009 at 9:37 PM, Peng Yu wrote: > What's the title? Introduction to Probability. > > On Fri, Oct 16, 2009 at 8:16 PM, Yi Du wrote: >> Hogg's book is enough for you considering your problems. >> >> Yi >> >> On Fri, Oct 16, 2009 at 7:12 PM, Peng Yu wrote: >>> >>> I need to refres

Re: [R] Recommendation on a probability textbook (conditional probability)

What's the title? On Fri, Oct 16, 2009 at 8:16 PM, Yi Du wrote: > Hogg's book is enough for you considering your problems. > > Yi > > On Fri, Oct 16, 2009 at 7:12 PM, Peng Yu wrote: >> >> I need to refresh my memory on Probability Theory, especially on >> conditional probability. In particular,

[R] Recommendation on a probability textbook (conditional probability)

I need to refresh my memory on Probability Theory, especially on conditional probability. In particular, I want to solve the following two problems. Can somebody point me some good books on Probability Theory? Thank you! 1. Z=X+Y, where X and Y are independent random variables and their distributi

Re: [R] Odp: : Question about correlation between data.

Petr PIKAL wrote: Hi r-help-boun...@r-project.org napsal dne 16.10.2009 15:24:05: hi everybody, I'm a student, and I'm new using R! I'm looking for statistical help hoping somebody can answer me! This is my problem: I have 2 temporal series. The firstone is a series of mesured data (heigh

Re: [R] Frequencies, proportions & cumulative proportions

Muenchen, Robert A (Bob) wrote: Ted, I know how to do that. It's just such a standard display in SAS, SPSS and Stata that I figured someone had done it and I had just overlooked it. Thanks! Bob I don't think there is a ready-made one, but it is very little effort to make your own: mkMyTable

[R] custom selfStart model works with getInitial but not nls

Hello, I'm having problems creating and using a selfStart model with nlme. Briefly, I've defined the model, a selfStart object, and then combined them to make a selfStart.default model. If I apply getInitial to the selfStart model, I get results. However, if I try usint it with nls or nlsL

[R] ozone data

Hi just a quick inquiry Sorry if this is too trivial but I am just a beginner in R. I am looking at the plyr package and I am intrigued at how data(i.e. ozone, baseball) is loaded without having to type data(ozone). Are they automatically loaded when i call library(plyr)? I want to do the same

[R] RODBC sqlSave does not append the records to a DB2 table

I am running R version 2.9.2 on Windows XP OS with RODBC version Version: 1.3-0. Has anyone out there in the R user community successfully appended records to a DB2 table on a remote database using the sqlSave function in the RODBC package? (or by any other means from R?) I posed a similar que

Re: [R] Breusch-pagan and white test - check homoscedasticity

On Fri, 16 Oct 2009, Gautier RENAULT wrote: Hi r-programmers, I performe Breusch-Pagan tests (bptest in package lmtest) to check the homoscedasticity of the residuals from a linear model and I carry out carry out White's test via bptest (formula, ~ x * z + I(x^2) + I(z^2)) include all regressor

Re: [R] Frequencies, proportions & cumulative proportions

Ted, I know how to do that. It's just such a standard display in SAS, SPSS and Stata that I figured someone had done it and I had just overlooked it. Thanks! Bob I don't think there is a ready-made one, but it is very little effort to make your own: mkMyTable <- function(X){ Table <- data.f

Re: [R] Matrixes as data

Create a matrix out of a list. In this example column 1 contains matrices, column 2 contains the number 1 and 2 and column 3 contains letters: > L <- matrix(list(m, m+10, 1, 2, "a", "b"), 2); L [,1] [,2] [,3] [1,] Integer,4 1"a" [2,] Numeric,4 2"b" > L[[2,1]] [,1] [,2] [1,]

Re: [R] Proper syntax for using varConstPower in nlme

Hi Dieter, Thanks for the reply. I had played with the initial conditions, but apparently not enough. I finally found some that avoided the singularity issue. Many thanks. More generally, I went over the documentation yet again in P&B and I still find it a bit confusing. They talk about

Re: [R] Frequencies, proportions & cumulative proportions

On 16-Oct-09 20:51:06, Muenchen, Robert A (Bob) wrote: > Dear R-Helpers, > > I've looked high and low for a function that provides frequencies, > proportions and cumulative proportions side-by-side. Below is the table > I need. Is there a function that already does it? > > Thanks, > Bob > >> # G

Re: [R] Matrixes as data

Magnus Torfason-2 wrote: > > > Hola! > > > > I am working on a problem where data points are (square) matrices. Is > > there a way to make a "vector" of matrices, such that it can be stored > > in a data.frame? > > I agree with previous posters that in most cases, you would want to > sto

Re: [R] [stats-rosuda-devel] how to install JGR manually?

On Oct 16, 2009, at 16:49 , Liviu Andronic wrote: On 10/16/09, c...@witthoft.com wrote: It doesn't get that far. I launch jgr.exe and it announces that it can't connect to the internet, and shuts itself down. In R, does the following work? require(JGR) JGR() No. JGR() creates a start

[R] Modify the plot from countCDFxt

Hello, I am running the plot from countCDFxt (popbio). I would like to report the y-axis as a percent instead of the log scale (e^01...). I can add an axis with axis(2, 0:1, line =2) but I'm having trouble understanding how to assign the tic marks (with 'at ='). I'd like to tell it to make tics wi

[R] Frequencies, proportions & cumulative proportions

Dear R-Helpers, I've looked high and low for a function that provides frequencies, proportions and cumulative proportions side-by-side. Below is the table I need. Is there a function that already does it? Thanks, Bob > # Generate some test scores > myValues <- c(70:95) > Score <- ( sample( myVal

Re: [R] how to install JGR manually?

On 10/16/09, c...@witthoft.com wrote: > It doesn't get that far. I launch jgr.exe and it announces that it can't > connect to the internet, and shuts itself down. > In R, does the following work? require(JGR) JGR() Liviu __ R-help@r-project.org maili

[R] what's the R code for wavelet decomposition (Haar transformation)?

Dear all, Using R function "dwt", it seems that I cannot specify the wavelet transformation like Haar. What's the R code for wavelet decomposition which allows me to specify Haar wavelet transformation? Of course, if it can include "db2", that is even better. In general, I want an R function

Re: [R] How to right-align labels in dotchart

Thanks very much John, that works. I had tried using las=1 inside the dotchart function itself, but it seems to be one of those occasions where the par parameter is over-ridden by the main plotting function. You can see the results here: http://www.stubbornmule.net/2009/10/asylum-seekers/ Now if on

[R] Breusch-pagan and white test - check homoscedasticity

Hi r-programmers, I performe Breusch-Pagan tests (bptest in package lmtest) to check the homoscedasticity of the residuals from a linear model and I carry out carry out White's test via bptest (formula, ~ x * z + I(x^2) + I(z^2)) include all regressors and the squares/cross-products in the auxilia

Re: [R] How can I run a function to a piece of text?

An alternative could be: do.call(mean, as.list(parse(text = unlist(strsplit(script, ", ") On Fri, Oct 16, 2009 at 3:39 PM, Javier PB wrote: > > Dear users, > > I got really stuck trying to apply a function to a piece of code that I > created using different string functions. > > To make thi

Re: [R] how to install JGR manually?

It doesn't get that far. I launch jgr.exe and it announces that it can't connect to the internet, and shuts itself down. I've installed every dependent package listed in the jgr documentation. Carl Oct 16, 2009 09:14:38 AM, landronim...@gmail.com wrote: ==

Re: [R] How can I run a function to a piece of text?

Javier PB wrote: > > I got really stuck trying to apply a function to a piece of code > that I created using different string functions. > > To make things really easy, this is a wee example: > > x<-c(1:10) > script<-"x, trim = 0, na.rm = FALSE"##script created by a > number of paste()

Re: [R] How odds ratio is computed in fisher.test()?

From ?fisher.test: estimate: an estimate of the odds ratio. Note that the _conditional_ Maximum Likelihood Estimate (MLE) rather than the unconditional MLE (the sample odds ratio) is used. Only present in the 2 by 2 case. from fisher.test: mle <- function(

Re: [R] How can I run a function to a piece of text?

Based solely on what you told us, this can be done using eval(parse(text=...)) cmd <- sprintf("mean(%s)", script) eval(parse(text = cmd)) However, with more context, there may be a better solution. See, for example, install.packages("fortunes") library(fortunes) fortune("parse()") HTH, --sunda

[R] How can I run a function to a piece of text?

Dear users, I got really stuck trying to apply a function to a piece of code that I created using different string functions. To make things really easy, this is a wee example: x<-c(1:10) script<-"x, trim = 0, na.rm = FALSE"##script created by a number of paste() and rep() steps mean(

Re: [R] Cannot calculate mean() for double vector

It looks like I should read the documentation a bit more carefully. Simply ignoring the NA values when calculating the mean *is* a much better solution. (Although having 10 values out of a vector of 4 or so set to zero is *not* going to bias the mean toward zero in any significant way.) Thanks

Re: [R] Matrixes as data

> Hola! > > I am working on a problem where data points are (square) matrices. Is > there a way to make a "vector" of matrices, such that it can be stored > in a data.frame? I agree with previous posters that in most cases, you would want to store matrices in a list. However, if you already have

Re: [R] Cannot calculate mean() for double vector

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Reuben Bellika > Sent: Friday, October 16, 2009 11:02 AM > To: r-help@r-project.org > Subject: Re: [R] Cannot calculate mean() for double vector > > OK. It looks like I just have

[R] How odds ratio is computed in fisher.test()?

I'm wondering how odds ratio is computed. I thought that it is (n11/n12)/(n21/n22), but it is not what fisher.test() computes. Could somebody let me know? > n11=3 > n12=1 > n21=1 > n22=3 > > n1_=n11+n12 > n2_=n21+n22 > > n_1=n11+n21 > n_2=n12+n22 > > x=rbind(c(n11,n12),c(n21,n22)) > > threshold=dh

Re: [R] Cannot calculate mean() for double vector

Are you sure that's the right solution? If you set them to 0, those values are averaged in with the others, and that could make a substantial difference. A much better approach: mean(x_ema, na.rm=TRUE) - see the help for mean for more information. Sarah On Fri, Oct 16, 2009 at 2:01 PM, Reuben

Re: [R] Cannot calculate mean() for double vector

On Fri, Oct 16, 2009 at 2:01 PM, Reuben Bellika wrote: > OK. It looks like I just have several NA values at the start of my array: > >> which (is.na(x_ema)) > [1] 1 2 3 4 5 6 7 8 9 > > That make sense, because the moving average is not defined for those > positions. I'll just have to set those val

Re: [R] How to right-align labels in dotchart

Oops replying to the wrong post but anyway does this do what you want? aa <- c(3,6,3,5,8) lbs <- c('cat','goat', 'elephant', 'horse', 'whale') dotchart(aa, pch=(16), col = 1:5, main="A Dotchart") axis(side = 2, seq_along(aa), lbs, las=1) --- On Thu, 10/15/09, Sean Carmody wrote: > From: Sean

Re: [R] Different way of scaling data

library(reshape) ?rescaler I think something along the lines of rescaler(data.frame, type="range") should do what you want. --- On Fri, 10/16/09, Noah Silverman wrote: > From: Noah Silverman > Subject: [R] Different way of scaling data > To: "r help" > Received: Friday, October 16, 2009,

Re: [R] Cannot calculate mean() for double vector

OK. It looks like I just have several NA values at the start of my array: > which (is.na(x_ema)) [1] 1 2 3 4 5 6 7 8 9 That make sense, because the moving average is not defined for those positions. I'll just have to set those values to zero: > x_ema = replace(x_ema, which(is.na(x_ema)), 0) > wh

Re: [R] "negative length vectors are not allowed" in wilcox.exact() and perm.test()

Number of politicians' promises kept? --- On Fri, 10/16/09, Patrick Burns wrote: > From: Patrick Burns > Subject: Re: [R] "negative length vectors are not allowed" in wilcox.exact() > and perm.test() > To: "Ben Bolker" > Cc: r-help@r-project.org > Received: Friday, October 16, 2009, 1:45 PM >

Re: [R] "negative length vectors are not allowed" in wilcox.exact() and perm.test()

Ben Bolker wrote: David Croll wrote: I want to compare two datasets and I get the message Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col), : negative length vectors are not allowed after specifying the exact test. I'm using the exactRankTests package. Do you suggest

[R] Different way of scaling data

Hi, I have a data.frame that I need to scale. I've been using the scale function and it works nicely. Some of the libraries I'm testing won't accept negative values for data, so I need to find a way to scale the data from 0 to 1 Any ideas? Thans! ___

Re: [R] Cannot calculate mean() for double vector

Reuben Bellika wrote: > > I've been using R recently to analyze some data, but I'm having a > problem using the mean() function. > > I imported the original data set as a vector of integers, x and then > calculated a exponential moving average of the data, x_ema. This part > worked fine. > > The

[R] DCC Garch model estimation

Hi all: I am a master student writing my master thesis about using the DCC garch model analysing the correlation between the insurance stock market and the local market. Right now I have a question about how to estimate the DCC garch, I know all the other meanings of the arguments in the

Re: [R] Cannot calculate mean() for double vector

Have you tried: mean(x) mean(as.numeric(as.character(x))) mean(x_ema) mean(as.numeric(as.character(x_ema))) What is the result of the following: which(is.na(as.numeric(as.character(x_ema Abit hard since you don't provide the data, but there may be an NA or character value that is causing the

[R] Problem with the stl function

Hi there, My name is Renan X. Cortes, student of Statistics, from south of Brazil, and I'd like to ask you a few questions about decomposition of time series. In R, when I fit the decomposition using the "stl" function, an object is returned when ask the summary of the fit, called STL.

Re: [R] Division of data frame and deletion of values from column

?subset ?split --- On Fri, 10/16/09, Joel Fürstenberg-Hägg wrote: > From: Joel Fürstenberg-Hägg > Subject: [R] Division of data frame and deletion of values from column > To: r-help@r-project.org > Received: Friday, October 16, 2009, 12:16 PM > > Hi all, > > I guess this might be an easy que

Re: [R] Generating a stochastic matrix with a specified second dominant eigenvalue

ooops! I tested it on a 3x3, and a 4x4, and it worked both times. Then I ran it once more and pasted it in without looking carefully... Here's another half-baked idea. Generate a random graph with given degree distribution, turn the incidence matrix into a stochastic matrix (ie the transition pr

Re: [R] Subset returning unexpected result

Steve Murray wrote: Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable MIME-Version: 1.0 Bill=2C It seems to be 'character' - odd...! Not really very "odd". I'll wager that your data was originally recorded in Excel or equivalent and that the column

Re: [R] Test model for singular gradient matrix

dmhultst wrote: > > > Hello, > > I am working with a real-time hydrologic modeling system, and I am using R > (R batch script on Linux) to create a non-linear relationship > (exponential) > between observed data. I want to extract the non-linear coefficients (“b” > and “m”) if the relationsh

[R] Cannot calculate mean() for double vector

I've been using R recently to analyze some data, but I'm having a problem using the mean() function. I imported the original data set as a vector of integers, x and then calculated a exponential moving average of the data, x_ema. This part worked fine. Then, I tried to find the mean squared error

Re: [R] generalization of tabulate()

If I have correctly understood what Robin wants, I don't think Gabor's elaborate solution is necessary (though I grant that it may be more general). It's also slow due to the "apply's". A more straightforward and faster approach is to convert the rows to individual character strings via paste() and

Re: [R] Division of data frame and deletion of values from column

On Oct 16, 2009, at 12:16 PM, Joel Fürstenberg-Hägg wrote: Hi all, I guess this might be an easy question, but I've searched multiple help pages without finding any answear... so now I put my trust in you! I have a data frame (36 variables and 556 observations). One column contains t

Re: [R] Matrixes as data

Thanks. The points of having the column of matrices (all the same dimension) in a data.frame, is that there are also other data, each matrix is at a location, so there are geographical coordinates and possibly other measurements at the same location. Kjetil On Fri, Oct 16, 2009 at 12:46 PM, Barry

[R] Division of data frame and deletion of values from column

Hi all, I guess this might be an easy question, but I've searched multiple help pages without finding any answear... so now I put my trust in you! I have a data frame (36 variables and 556 observations). One column contains three factors, and I would like to divide the data frame into three n

Re: [R] tapply() and using factor() on a factor

On Oct 16, 2009, at 11:33 AM, Alexander Peterhansl wrote: Thank you Mohamed and Bill for your replies. (I did not send the data because it is unwieldy.) Yes Bill, the issue arises directly from what you had guessed. I was working with a subset of the data (which implicitly had factors for

Re: [R] Removing Embedded Null characters from text/html

[David contacted me directly, so I am sending my off-line reply to the list just for the record in case others encounter a simple problem.] Hi David. No problem contacting me at all. I saw your mail at one point on the mailing list, but didn't have a chance to respond. Indeed, it seems like th

[R] Test model for singular gradient matrix

Hello, I am working with a real-time hydrologic modeling system, and I am using R (R batch script on Linux) to create a non-linear relationship (exponential) between observed data. I want to extract the non-linear coefficients (“b” and “m”) if the relationship can be created, if the relationship

Re: [R] Matrixes as data

Bonjour, On Oct 16, 2009, at 11:36 AM, Kjetil Halvorsen wrote: Hola! I am working on a problem where data points are (square) matrices. Is there a way to make a "vector" of matrices, such that it can be stored in a data.frame? Can that be done with S3? or do I have to learn S4 objects & method

Re: [R] populating an array

R doesn't access arrays like C, use [i,j] to access a 2-d array, e.g.: my_array <- array(0,dim=c(2,2)) for(i in seq(1,2,by=1)){ + for(j in seq(1,2,by=1)){ + my_array[i,j] = i+j + } + } my_array [,1] [,2] [1,]23 [2,]34 tdm wrote: Hi, Can someone please give m

Re: [R] Matrixes as data

On Fri, Oct 16, 2009 at 4:36 PM, Kjetil Halvorsen wrote: > Hola! > > I am working on a problem where data points are (square) matrices. Is > there a way to make a > "vector" of matrices, such that it can be stored in a data.frame? Can > that be done with S3? > or do I have to learn S4 objects & me

Re: [R] Proper syntax for using varConstPower in nlme

Michael A. Gilchrist wrote: > > - >> nlme(Count ~ quad.PBMC.model(aL, aN, T0), > + data = tissueData, > + weights = varConstPower(form =~ Count), > + start = list( fixed = c(rep(1000, 8), -2, -2) ), > + fixed = list(T0 ~ TypeTissue-1, aL ~ 1,

[R] Matrixes as data

Hola! I am working on a problem where data points are (square) matrices. Is there a way to make a "vector" of matrices, such that it can be stored in a data.frame? Can that be done with S3? or do I have to learn S4 objects & methods? Kjetil __ R-help@r

[R] Determining a linear model based on a factor

I guess I should disclose up front that am not a statistician by schooling, but I am intersted in getting the terminology correct so please correct it if I butcher it too badly. I have been able to very easily build a linear model showing the correlation between two variables, e.g. year built

Re: [R] tapply() and using factor() on a factor

Thank you Mohamed and Bill for your replies. (I did not send the data because it is unwieldy.) Yes Bill, the issue arises directly from what you had guessed. I was working with a subset of the data (which implicitly had factors for the complete data set). On this, what is the best way take a su

Re: [R] "negative length vectors are not allowed" in wilcox.exact() and perm.test()

David Croll wrote: > > I want to compare two datasets and I get the message > > Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col), : > negative length vectors are not allowed > > after specifying the exact test. I'm using the exactRankTests package. Do > you suggest me

Re: [R] "negative length vectors are not allowed" in wilcox.exact() and perm.test()

David Croll wrote: > > I want to compare two datasets and I get the message > > Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col), : > negative length vectors are not allowed > > after specifying the exact test. I'm using the exactRankTests package. Do > you suggest me

Re: [R] Generating a stochastic matrix with a specified second dominant eigenvalue

On Oct 15, 2009, at 6:24 PM, Ravi Varadhan wrote: Hi, Given a positive integer N, and a real number \lambda such that 0 < \lambda < 1, I would like to generate an N by N stochastic matrix (a matrix with all the rows summing to 1), such that it has the second largest eigenvalue equal to

Re: [R] package installation from source

On 10/16/2009 9:31 AM, Petr PIKAL wrote: Dear all I noticed from NEWS 2.11.0,dev SIGNIFICANT USER-VISIBLE CHANGES o Packages must have been installed under R 2.10.0 or later, as the current help system is the only one now supported. So I tried to follow instructions in manual, Du

Re: [R] using a custom color sequence for image()

Rajarshi Guha asked: > > Hi, I'd like to use a custom color sequence (black - low values, > green - high values) in am image() plot. While I can specify colors > (say a sequence of grays) to the col argument, the ordering is > getting messed up. I have two questions: > > 1. How can I get a seq

Re: [R] Generating a stochastic matrix with a specified second dominant eigenvalue

A valiant attempt, Albyn! Unfortunately, the matrix B is not guaranteed to be a stochastic matrix. In fact, it is not even guaranteed to be a real matrix. Your procedure can generate a B that contains negative elements or even complex elements. > M = matrix(runif(9),nrow=3) > M = M/apply(

[R] "negative length vectors are not allowed" in wilcox.exact() and perm.test()

Dear R friends, I want to compare two datasets and I get the message Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col), : negative length vectors are not allowed after specifying the exact test. I'm using the exactRankTests package. Do you suggest me using the coin lib

Re: [R] [stats-rosuda-devel] how to install JGR manually?

On Oct 16, 2009, at 9:14 , Liviu Andronic wrote: (cc'ing JGR specific list) Yes, thanks, stats-rosuda-devel is where this post should go. Carl, On 10/15/09, Carl Witthoft wrote: Here's the problem: on Windows, the 'jgr.exe' tool starts up by checking for a connecting to the 'net in order

[R] using a custom color sequence for image()

Hi, I'd like to use a custom color sequence (black - low values, green - high values) in am image() plot. While I can specify colors (say a sequence of grays) to the col argument, the ordering is getting messed up. I have two questions: 1. How can I get a sequence of say 256 colors starting from

Re: [R] two graphs 1 x-axis

You can use ggplot2. library(ggplot2) a<-c(1,2,3,4,5,6) b<-c(3,5,4,6,1,1) c<-c(1,1,1,1,1,1) dframe = data.frame(a,b,c);dframe melt.dframe <- melt(dframe, id= "a");melt.dframe qplot(a,value,data=melt.dframe) + facet_grid(variable~.,scales="free") Felipe D. Carrillo Supervisory Fishery Biologis

Re: [R] currency conversion function?

On Fri, Oct 16, 2009 at 8:04 AM, Liviu Andronic wrote: > On 10/15/09, Jeff Ryan wrote: >>  > getFX("EUR/USD",from="2009-04-01") >> > Indeed, with the date correctly specified, the function no longer > generates errors. There is one issue though (similar to the one in the > code posted by Henrique

Re: [R] reference on fisher.test()

And indeed Agresti's book (first edition on the help page, also has a 2002 second edition) is a good source for the 'minutiae'. I find the description in his 1992 survey article even better, especially with regard to the noncentral hypergeometric distribution: @article{Agresti:92,

[R] Odp: : Question about correlation between data.

Hi r-help-boun...@r-project.org napsal dne 16.10.2009 15:24:05: > hi everybody, I'm a student, and I'm new using R! > I'm looking for statistical > help hoping somebody can answer me! > > This is my problem: > I have 2 temporal > series. The firstone is a series of mesured data (height of m

[R] package installation from source

Dear all I noticed from NEWS 2.11.0,dev SIGNIFICANT USER-VISIBLE CHANGES o Packages must have been installed under R 2.10.0 or later, as the current help system is the only one now supported. So I tried to follow instructions in manual, Duncan Murdoch presentation and help pages t

[R] Confidence intervals and p-values of ARMAX model

Is there any way to get confidence intervals and/or p-values in the ARMAX model? I tried the search using: ?ar ?arma ?arima But, AFAIK, "ar" and "arma" don't accept exogenous (independent) variables, and "arima(x, order, xreg)" does not return statistics on the xreg terms. BTW, just in case so

[R] : Question about correlation between data.

hi everybody, I'm a student, and I'm new using R! I'm looking for statistical help hoping somebody can answer me! This is my problem: I have 2 temporal series. The firstone is a series of mesured data (height of monitorated points), the second is a series of temperature (in Celsius degree).

Re: [R] how to install JGR manually?

(cc'ing JGR specific list) Hello On 10/15/09, Carl Witthoft wrote: > Here's the problem: on Windows, the 'jgr.exe' tool starts up by checking > for a connecting to the 'net in order to grab the support packages. Well, > we have machines at work that are not and never will be connected to the > In

Re: [R] currency conversion function?

On 10/15/09, Jeff Ryan wrote: > > getFX("EUR/USD",from="2009-04-01") > Indeed, with the date correctly specified, the function no longer generates errors. There is one issue though (similar to the one in the code posted by Henrique): > getFX("EUR/PEN",from="2009-10-16") [1] "EURPEN" > EURPEN

Re: [R] reference on fisher.test()

For me, the classical reference for FET is: @BOOK{Siegel56, title = {Nonparametric Statistics for the Behavioral Sciences}, publisher = {McGraw-Hill}, year = {1956}, author = {Siegel, Sidney}, address = {New York} } Tom Prof Brian Ripley wrote: On Fri, 16 Oct 2009, Robin Hankin wrote: H

Re: [R] two graphs 1 x-axis

Hi Naomi, Take a look at the lattice package for plotting. An example using your data: library(lattice) library(reshape) a<-c(1,2,3,4,5,6) b<-c(3,5,4,6,1,1) c<-c(1,1,1,1,1,1) bla = data.frame(a,b,c) # melt is from reshape bla2 = melt(bla, id.vars = "a") xyplot(value~a | variable, bla2, layout

Re: [R] reference on fisher.test()

For some alternative views (and references) for FET see: http://www.stat.columbia.edu/~cook/movabletype/archives/2009/10/what_is_the_bay.html#comments kjetil On Fri, Oct 16, 2009 at 8:38 AM, Prof Brian Ripley wrote: > On Fri, 16 Oct 2009, Robin Hankin wrote: > >> Hi >> >> fexact.c points you to

Re: [R] Estimation in a changepoint regression with R

Package bcp does Bayesian changepoint analysis, though not in the general regression framework.  The most recent reference is Bioinformatics 24(19) 2143-2148; doi:   10.1093/bioinformatics/btn404; slightly older is JSS 23(3). Both reference some alternatives you might want to consider (including s

Re: [R] currency conversion function?

On 10/15/09, Henrique Dallazuanna wrote: > foo <- function(from, to, date){ >url <- > "http://www.oanda.com/convert/classic?script=..%2Fconvert%2Fclassic&language=en&value=1"; > params <- > > sprintf("%s&date=%s&exch=%s&exch2=&margin_fixed=0&expr=%s&expr2=&SUBMIT=Convert+Now&lang=

Re: [R] Urgent help requested to modify a script

On 10/16/2009 02:01 PM, one2luv wrote: I am hoping someone will tak up this chalenge (I am new to R) I have inheritied an R script but need to change it. The script currently includes hardcoded file locations on lines 12,166 and 167. I need to modify this script to allow the folder to be passed

[R] Urgent help requested to modify a script

I am hoping someone will tak up this chalenge (I am new to R) I have inheritied an R script but need to change it. The script currently includes hardcoded file locations on lines 12,166 and 167. I need to modify this script to allow the folder to be passed as a command line argument to Rscript.ex

Re: [R] reference on fisher.test()

On Fri, 16 Oct 2009, Robin Hankin wrote: Hi fexact.c points you to the original ACM paper: Well, you'll get a better idea from the help page as to the real 'original' source reference: the reference below is to a revised version in a remark. And indeed Agresti's book (first edition on the

Re: [R] generalization of tabulate()

Using the generalized inner product defined in this post: https://www.stat.math.ethz.ch/pipermail/r-help/2006-July/109311.html try this: cbind(S, d = rowSums(inner(S, obs, identical))) On Fri, Oct 16, 2009 at 4:29 AM, Robin Hankin wrote: > Hi > > I want a generalization of tabulate() wh

Re: [R] two graphs 1 x-axis

On Fri, Oct 16, 2009 at 12:22:06PM +0200, Duijvesteijn, Naomi wrote: > I have a question concerning plotting graphs. > Here an example dataset > > > a<-c(1,2,3,4,5,6) > b<-c(3,5,4,6,1,1) > c<-c(1,1,1,1,1,1) > d<-as.data.frame(cbind(a,b,c)) > plot.new() > plot(d$a, d$b, col="red") > par(new=TRUE)

Re: [R] two graphs 1 x-axis

On 10/16/2009 09:22 PM, Duijvesteijn, Naomi wrote: Dear R-people I have a question concerning plotting graphs. Here an example dataset a<-c(1,2,3,4,5,6) b<-c(3,5,4,6,1,1) c<-c(1,1,1,1,1,1) d<-as.data.frame(cbind(a,b,c)) plot.new() plot(d$a, d$b, col="red") par(new=TRUE) plot(d$a,d$c, col="red"

[R] two graphs 1 x-axis

Dear R-people I have a question concerning plotting graphs. Here an example dataset a<-c(1,2,3,4,5,6) b<-c(3,5,4,6,1,1) c<-c(1,1,1,1,1,1) d<-as.data.frame(cbind(a,b,c)) plot.new() plot(d$a, d$b, col="red") par(new=TRUE) plot(d$a,d$c, col="red", pch="|") What I would want is to plot de second p

  1   2   >