Re: [R] Warning message

Hi: A univariate t-distribution with 1 df is equivalent to the standard Cauchy distribution, for which it is well established that the population mean does not exist. You're basically simulating a vector version of the Cauchy distribution (an IID 8-dimensional version), so the same problem is like

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 met

[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) x$count

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] 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 te

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 mailin

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 gu

[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 -0.0255661

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

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 se

[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 argume

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 intern

[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 c

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] 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 me

[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 ba

[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: http://r.789695.n4.nabble.com/update-R-2-11-0-there-is-error-wh

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-help@r-project.

Re: [R] How to estimate whether overfitting?

On Sun, May 9, 2010 at 11: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 R

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 N

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. One

[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 ge

Re: [R] P values

On Sun, May 9, 2010 at 6:53 PM, 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 sounds better in french: Une variable al

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] 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 i

[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 fixed-effe

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

Hi: On Sun, May 9, 2010 at 5: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

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... library(ggpl

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 and

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

[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 dele

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. > [

[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 > l

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 i

[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] 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 t

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 o

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 tr

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 05/09/

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 li

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

On May 9, 2010, at 12:30 AM, Venkatesh Patel wrote: -- Forwarded message -- From: Dr. Venkatesh 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

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 us

[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: olgaformig...@gmail.com MSN: o

[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

[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 overfi

[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] Fwd: R apply() help -urgent

-- Forwarded message -- From: Dr. Venkatesh 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 for a different var

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" To: "George Chen"