Re: [R] extract required data from already read data

On Thu, May 6, 2010 at 10:47 PM, Mohan L l.mohan...@gmail.com wrote: snip I am using source() to run program, like source (sample.R) I want to know how to pass the argument to my sample.R, some this like source(sample.R, xxx,FeMale)  Is it possible? any link will greatly appreciated.

Re: [R] Is it possible to rearrange the facets in ggplot while keeping everything else the same?

Hi Dennis, It is similar to what I was trying to achieve. Is there a way to show the grouped time courses as separate facets. So in lattice, the time courses would be conditioned on Ab status and subject name at the same time. Thanks. George From: Dennis Murphy djmu...@gmail.com

[R] Fwd: R apply() help -urgent

-- Forwarded message -- From: Dr. Venkatesh drve...@liv.ac.uk Date: Sun, May 9, 2010 at 4:55 AM Subject: R apply() help -urgent To: r-help@r-project.org I have a file with 4873 rows of 1s or 0s and has 26 alphabets (A-Z) as columns. the 27th column also has 1s and 0s but stands

[R] reshape multiple imputed dataset from wide to long

Dear list, I did multiple imputation with mice on a wide dataset. For my repeated measures ANOVA I need to reshape the data from wide to long. In order to pool the analysis, I would still like it to be a mids-object after reshaping. Does anyone know how to do this? Thanks, Nicole

[R] How to estimate whether overfitting?

1. is there some criterion to estimate overfitting? e.g. R2 and Q2 in the training set, as well as R2 in the test set, when means overfitting. for example, in my data, I have R2=0.94 for the training set and for the test set R2=0.70, is overfitting? 2. in this scatter, can one say this

[R] Question about factor that is numeric, in aov()

I notice something curious about how aov() treats a numeric factor: score is a dependent variable and group is a factor in a one-way ANOVA. But group contains numeric codes and is not a factor (checked with is.factor). An ANOVA done using: aov(score~factor(group), data=mydata) gives the right

[R] U de Theil

Hi I need to calculate the indicator U de Theil from my forecasts, so I want to know if there is a specific command or function that provides me the indicator U de Theil Thanks.. Olga Maria Formigoni C. Walter E-mail: mailto:olgaformig...@gmail.com olgaformig...@gmail.com MSN:

Re: [R] Question about factor that is numeric, in aov()

On May 9, 2010, at 8:36 AM, Ravi Kulkarni wrote: I notice something curious about how aov() treats a numeric factor: score is a dependent variable and group is a factor in a one-way ANOVA. But group contains numeric codes and is not a factor (checked with is.factor). An ANOVA done using:

Re: [R] Fwd: R apply() help -urgent

On May 9, 2010, at 12:30 AM, Venkatesh Patel wrote: -- Forwarded message -- From: Dr. Venkatesh drve...@liv.ac.uk Date: Sun, May 9, 2010 at 4:55 AM Subject: R apply() help -urgent To: r-help@r-project.org I have a file with 4873 rows of 1s or 0s and has 26 alphabets (A-Z) as

Re: [R] U de Theil

I do not know if this is what you want, but take a look at the External links at the bottom of the Wikipedia page: http://en.wikipedia.org/wiki/Theil_index Ravi -- View this message in context: http://r.789695.n4.nabble.com/U-de-Theil-tp2164441p2164458.html Sent from the R help mailing

Re: [R] Fwd: R apply() help -urgent

Hi, Instead of accessing the data by its name in data.frame, you may use index by [[ ]] as well. For your case, the looping may looks like for (i in 1:26){ datatable - table(data1[[i+2]], data1$pLoss) #create a new datatable2 or 3 fisher.test(datatable) } Bests, Ruihong On

Re: [R] U de Theil

On May 9, 2010, at 8:24 AM, Olga Maria Formigoni Carvalho Walter wrote: Hi I need to calculate the indicator U de Theil from my forecasts, so I want to know if there is a specific command or function that provides me the indicator U de Theil It appears to me that you have not yet

Re: [R] How to estimate whether overfitting?

On May 9, 2010, at 9:20 AM, bbslover wrote: 1. is there some criterion to estimate overfitting? e.g. R2 and Q2 in the training set, as well as R2 in the test set, when means overfitting. for example, in my data, I have R2=0.94 for the training set and for the test set R2=0.70, is

Re: [R] P values

Thank you for your replies. As I said (wrote) before, 'I am no statistician'. But I think I know what Random Variables are (not). Random variables are not random, neither are they variable. [It sounds better in french: Une variable aléatoire n'est pas variable, et n'a rien d'aléatoire.] See

[R] changing parameters of the box and whisker plot

Hello: I am plotting some data using the box and whisker plot. However, I only want to plot the median, max and min, as I only have these values and not the quartile values. It seems R arbitrarily constructs the box margins to be halfway between the median and the max/min. How do I make the

Re: [R] Question about factor that is numeric, in aov()

Dear Ravi, On Sunday 09 May 2010, Ravi Kulkarni wrote: I notice something curious about how aov() treats a numeric factor: In R, there is no such thing as a numeric factor. A numeric vector is not a factor unless declared as such. score is a dependent variable and group is a factor in a

[R] Plotting log-axis with the exponential base to a plot with the default logarithm base 10

Hello! I have a problem which I have tried to solve for several days now.. I have plottet a lineplot.CI in the library sciplot, and I am trying to plot it with a logaritmic y-axis (with exponential base). The problem is that; when I type log y, the axis transforms into the logaritmic of

Re: [R] P values

Time to rescue Random Variables before they drown! On 09-May-10 16:53:02, Bak Kuss wrote: Thank you for your replies. As I said (wrote) before, 'I am no statistician'. But I think I know what Random Variables are (not). Random variables are not random, neither are they variable. [It

[R] Error with read.ssd

Dear all, I'm trying to read and sas file hen using read.ssd, I got this error message: symbolic links are not supported I have debug read.ssd function and the function file.symlink is used. I run windows XP on virtual box, hosted by Vista Home  [[alternative HTML version

Re: [R] Plotting log-axis with the exponential base to a plot wi

On 09-May-10 18:10:27, Elisabeth Bjerke Rastad wrote: Hello! I have a problem which I have tried to solve for several days now.. I have plottet a lineplot.CI in the library sciplot, and I am trying to plot it with a logaritmic y-axis (with exponential base). The problem is that; when I type

Re: [R] changing parameters of the box and whisker plot

On May 9, 2010, at 1:33 PM, michael westphal wrote: Hello: I am plotting some data using the box and whisker plot. What function? However, I only want to plot the median, max and min, as I only have these values and not the quartile values. You could just read the help page for bxp

Re: [R] changing parameters of the box and whisker plot

Hi: Here's one way to get what you want, using the plyr and ggplot2 packages. # Fake some data dd - data.frame(g = factor(rep(LETTERS[1:10], each = 30)), y = rnorm(300)) # Summarize to get min, median and max per group # Uses function ddply() in the plyr package...

Re: [R] Question about factor that is numeric, in aov()

Hi: On Sun, May 9, 2010 at 5:36 AM, Ravi Kulkarni ravi.k...@gmail.com wrote: I notice something curious about how aov() treats a numeric factor: score is a dependent variable and group is a factor in a one-way ANOVA. But group contains numeric codes and is not a factor (checked with

[R] non-linear estimation with many firm-specific parameters

Dear R experts--- I doubt that someone has already solved my problem, but I thought I would ask quickly, just in case someone has. Let' say I start with a (flattened panel) model that says y[i] = x[i] + b*(T-x[i]) easy enough---this is just a linear model. I could also make this a

Re: [R] Plotting log-axis with the exponential base to a plot wi

Hello Ted! Thank you a lot for your reply!!! I will try to explain again; what I want is a logarithmic scaled y-axis with the base of e (not 10). And the values I would like to use in the plot are the raw values (not transformed in any way). Do you still think that the base of log does not matter

Re: [R] Plotting log-axis with the exponential base to a plot wi

On 09-May-10 20:15:56, Elisabeth Bjerke Rastad wrote: Hello Ted! Thank you a lot for your reply!!! I will try to explain again; what I want is a logarithmic scaled y-axis with the base of e (not 10). And the values I would like to use in the plot are the raw values (not transformed in any

Re: [R] P values

On Sun, May 9, 2010 at 6:53 PM, Bak Kuss bakk...@gmail.com wrote: Thank you for your replies. As I said (wrote) before, 'I am no statistician'. But I think I know what Random Variables are (not). Random variables are not random, neither are they variable. [It sounds better in french: Une

[R] Plot polygon in 3D with rgl

Dear R-helpers, an rgl-ers in particular, what is the easiest way to plot a section of a plane in 3D, that is given by the xyz coordinates of the outline? Suppose I have a polygon - which I know for sure is a set of coordinates on the same plane. One method I found is to use surf.tri from the

Re: [R] Plot polygon in 3D with rgl

On 09/05/2010 6:41 PM, Remko Duursma wrote: Dear R-helpers, an rgl-ers in particular, what is the easiest way to plot a section of a plane in 3D, that is given by the xyz coordinates of the outline? Suppose I have a polygon - which I know for sure is a set of coordinates on the same plane.

Re: [R] Plot polygon in 3D with rgl

Duncan, thanks for the tip! I actually saw this 2D function but had no idea how to use it in 3D. Works great. Remko - Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond

Re: [R] How to estimate whether overfitting?

On Sun, May 9, 2010 at 11:53 AM, David Winsemius dwinsem...@comcast.net wrote: On May 9, 2010, at 9:20 AM, bbslover wrote: 1. is there some criterion to estimate overfitting?  e.g. R2 and Q2 in the training set, as well as R2 in the test set, when means overfitting.   for example,  in my

Re: [R] How to sort a grouped barchart?

Thanks for both. The reorder function works perfectly. -- View this message in context: http://r.789695.n4.nabble.com/How-to-sort-a-grouped-barchart-tp2133579p2164506.html Sent from the R help mailing list archive at Nabble.com. __

[R] update R 2.11.0，there is error when usi ng plot(), how can I do?

a-1:5 b-2:6 plot(a,b) Error in function (width, height, pointsize, record, rescale, xpinch, : Graphics API version mismatch before, R 2.10 , plot() is ok. Now, R 2.11.0 does not work -- View this message in context:

[R] Plotting log-axis with the exponential base to a plot with the default logarithm base 10

Hello! I have a problem which I have tried to solve for several days now.. I have plottet a lineplot.CI in the library sciplot, and I am trying to plot it with a logaritmic y-axis (with exponential base). The problem is that; when I type log y, the axis transforms into the logaritmic of base

[R] Dimensions on svychisq on svydesign

Dear Forum I a running svychisq from the survey package and get errors with the number of dimensions, errors that I do not understand and do not know how to fix. I ask you kindly to help me out. The eror message follows with some information below. I hope there are enough information to help

Re: [R] How to estimate whether overfitting?

thanks for your suggestion. many I need to learn indeed. I will buy that good book. kevin -- View this message in context: http://r.789695.n4.nabble.com/How-to-estimate-whether-overfitting-tp2164417p2164847.html Sent from the R help mailing list archive at Nabble.com.

[R] Random walk

Hi everybody, I am trying to generate two random walks with an specific correlation, for example, two random walks of 200 time steps with a correlation 0.7. I built the random walks with: x-cumsum(rnorm(200, mean=0,sd=1)) y-cumsum(rnorm(200, mean=0,sd=1)) but I don't know how to fix the

Re: [R] Fwd: R apply() help -urgent

Venkatesh Is this what you are looking for? # Example data df=data.frame(A=c(1,0,0,0,1),B=c(1,1,0,0,0),C=c(1,0,0,0,0),Val=c(1,0,1,0,1)) # Variation of code of David Winsemius tbl = lapply(df[, 1:3] , function(x) table(x, df$Val)) fet = lapply(tbl, function(x) fisher.test(x)) # Identify

[R] rtracker dependencies installation problem?

Why do I get this error message? library(rtracklayer) session - browserSession() x-ucscTableQuery(session, knownGene , GenomicRanges(0,12000,chr1)) Error in solveUserSEW0(start = start, end = end, width = width) : solving row 1: range cannot be determined from the supplied arguments (too

Re: [R] How to estimate whether overfitting?

On 05/09/2010 10:53 AM, David Winsemius wrote: On May 9, 2010, at 9:20 AM, bbslover wrote: 1. is there some criterion to estimate overfitting? e.g. R2 and Q2 in the training set, as well as R2 in the test set, when means overfitting. for example, in my data, I have R2=0.94 for the training

Re: [R] rtracker dependencies installation problem?

On 05/09/2010 07:11 PM, Sigal Blay wrote: Why do I get this error message? library(rtracklayer) session - browserSession() x-ucscTableQuery(session, knownGene , GenomicRanges(0,12000,chr1)) Error in solveUserSEW0(start = start, end = end, width = width) : solving row 1: range cannot

[R] tapply function with NA

Hi R users, I have a matrix m of the type: m X4.20.2010 X4.19.2010 X4.16.2010 [1,] 0.008319468 0. -0.008250825 [2,] 0.005574136 0.01816118 0.073081608 [3,] -0.047830688 0.01612903 -0.030239833 [4,] NA NA NA [5,] 0.008746356 0.02848576

Re: [R] tapply function with NA

It is exactly the same tmp - matrix(1:24,6,4) tmp[4,] - NA tmp apply(tmp, 2, sum, na.rm=TRUE) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] About counting letters.-

Hi All, I'm working with text processing and I have a problem. The problem is the following: I have a text, e.g.: 'R is a free software environment for statistical computing and graphics.' and I need create a vector of characters such that each position in the vector will be each letter in the

Re: [R] How to estimate whether overfitting?

thank you, I have downloaded it. studying -- View this message in context: http://r.789695.n4.nabble.com/How-to-estimate-whether-overfitting-tp2164417p2164932.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] About counting letters.-

tmp - 'R is a free software environment for statistical computing and graphics' strsplit(tmp, ) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] Warning message

Hello, I want to draw a histogram of the mean of sample observations drawn from multivariate t distribution. I am getting the following error corresponding to the code I used. Though I am getting the graph, but I am curious to know the warning message. Warning messages: 1: In if (freq)

Re: [R] How to estimate whether overfitting?

many thanks . I can try to use test set with 100 samples. anther question is that how can I rationally split my data to training set and test set? (training set with 108 samples, and test set with 100 samples) as I know, the test set should the same distribute to the training set. and what