[R] Repeated Measures ANOVA and Missing Values in the data set

I am doing Repeated Measures ANOVA with missing values. When i run my model i get this error message. *aov.out = aov(values ~ time + Error(subject/time), data=mydata2)Warning message:In aov(values ~ time + Error(subject/time), data = mydata2) : Error() model is singular* The missing Values are

[R] Repeated Measures ANOVA and the Bonferroni post hoc test different results of significantly

Hi All and thanks for Help I am doing an Repeated Measures ANOVA and the Bonferroni post hoc test for my data using R project. The ANOVA gives a significantly difference between the data but not the Bonferroni post hoc test. > anova(aov2) numDF denDF F-value p-value (In

[R] the right reference for the R Stats package for a scientific journal

Dear Members list, I am writing a paper for a research where i used "the R Stats package". No one knows the right reference for this package? Thanks in Advance Gianni [[alternative HTML version deleted]] __ R-help@r-project.org mailing list ht

[R] total number of citations for R project

Dear Member list, is there a weblink or a paper where the total number of citations for R project is report? Thanks in advance Gianni [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

[R] function to calculate a SMAPE (Symmetric mean absolute percentage error) to avoid the possibility of an inflation caused by zero values in the series

Dear Reseacher, i find this modified version of SMAPE at pag 13 formula "msSMAPE" to to avoid the possibility of an inflation caused by zero values in the series using a Si component in the denominator of the symmetric MAPE http://www.stat.iastate.edu/preprint/articles/2004-10.pdf I wrote a f

[R] help to find or fix a SMAPE (Symmetric mean absolute percentage error ) with correction to avoid zero value in the series

Dear Researchers, I am looking for a library or a function to calculate SMAPE with same doneminator to avoid a problem to the presence of 0 values in the series. I find this variotion of SMAPE called mSMAPE mSMAPE page 13 http://www.stat.iastate.edu/preprint/articles/2004-10.pdf yesterday nigt

[R] some help to improve "hist to plot relative frequencies"

Dear Researches, sorry for disturb. I wish to improve my figure in R plotting the relative frequencies of my data set. library(lattice) a <- c(0,0,0,1,1,2,4,5,6,7,7,7,7,7,8,8,8,8,9,9,9,9,10,10,11) histogram(a, xlab="myData") what i wish to do is: 1) invert the order of X and Y (eg: Precent of T

[R] special simbol (±) in a legend

dear reserachers, I am looking for a expression about a special symbol (±) in a Legend and add on front of "Mean ± SD" the symbol of bracket sorry if the example is not great disptest <- matrix(rnorm(10),nrow=10) disptest.means<- rowMeans(disptest) plot(1:10,disptest.means) dispersion(1:10,disp

[R] error to convert a Compute A^-1 B from Matlab to R using solve(A, B)

Dear Researchers, I need to convert the following equation in R from Matlab a = [x y ones(size(x))]; b = [-(x.^2+y.^2)]; a\b ans = -9.9981 -16.4966 -7.6646 my solution in R is: a = cbind(x,y,rep(1,length(x))) b = cbind(-(x^2+y^2)) > head(a) xy [1,] 14.45319 5.0657

[R] Fit circle with R

Dear Researchers, I wrote two function to fit a circle using noisy data. 1- the fitCircle() is derived from MATLAB code of * zhak Bucher* from the link http://www.mathworks.com/matlabcentral/fileexchange/5557-circle-fit/content/circfit.m 2- the CircleFitByPratt() from MATLAB code of *Nikolai Cher

[R] Based-Density cluster algorithms, Help where find some algorithms

Dear Researches, Did anyone know where to find the following Based-Density cluster algorithms in R or function implemented? OPTICS (Ordering Points To Identify the Clustering Structure) DECLUN (DENsity CLUstering) Thanks in Advance Gianni [[alternative HTML version deleted]] __

[R] problem to tunning RandomForest, an unexpected result

Dear Researches, I am using RF (in regression way) for analize several metrics extract from image. I am tuning RF setting a loop using different range of mtry, tree and nodesize using the lower value of MSE-OOB mtry from 1 to 5 nodesize from1 to 10 tree from 1 to 500 using this paper as refery

[R] tuning random forest. An unexpected result

Dear Researches, I am using RF (in regression way) for analize several metrics extract from image. I am tuning RF setting a loop using different range of mtry, tree and nodesize using the lower value of MSE-OOB mtry from 1 to 5 nodesize from1 to 10 tree from 1 to 500 using this paper as refery

[R] help to setting a multiple (linear) regression model with a 5% significance level (threshold) for the inclusion of the model variables.

Dear Researchers, someone know the right syntax to chose a 5% significance level (threshold) for the inclusion of the model variables in a multiple (linear) regression in backward way? I set the formula in this way, but I don't know to choose the 5% significance? lmodelV <- step(lm(formula=MyF

[R] explanation why RandomForest don't require a transformations (e.g. logarithmic) of variables

Dear Researches, sorry for the easy and common question. I am trying to justify the idea of RandomForest don't require a transformations (e.g. logarithmic) of variables, comparing this non parametrics method with e.g. the linear regressions. In leteruature to study my phenomena i need to apply a l

Re: [R] explanation why RandomForest don't require a transformations (e.g. logarithmic) of variables

others if the situation calls for it. > > Cheers, > Andy > > > > -Original Message- > > From: r-help-boun...@r-project.org > > [mailto:r-help-boun...@r-project.org] On Behalf Of gianni lavaredo > > Sent: Monday, December 05, 2011 1:41 PM > > To: r-hel

[R] help to split a ID column in a data.frame and create a new ID column

Dear Researchers, I have a data.frame with 2 columns like this: mydf <- data.frame(value=c(1,2,3,4,5),ID=c("Area_1","Area_2","Area_3","Area_4","Area_5")) > mydf value ID 1 1 Area_1 2 2 Area_2 3 3 Area_3 4 4 Area_4 5 5 Area_5 I need to convert the *ID *in the following

[R] plot a BARPLOT with sd deviation bar up and down

dear Researchers, i am looking for a function to plot a barplot for each mean value and the related standard deviation, and i can close my week. This is an example of my data set. really Thanks in advance for any help or suggestions Gianni My.mean <- data.frame(Mean=c(0.4108926,0.3949009,0.45

[R] Different result with "kruskal.test" and post-hoc analysis with Nemenyi-Damico-Wolfe-Dunn test implemented in the help page for oneway_test in the coin package that uses multcomp

Dear Researchers, Sorry for this email but I am not a statistician, and for this I have this problem to understand. Thanks in Advance for help and suggestions. Gianni I have 21 classes (00, 01, 02, 04, ,020) with different length. I did a kruskal wall test in R with the following code krusk

[R] help to slip a file name using "strsplit" function

Dear Researchers, I have several files as this example: Myfile_MyArea1_sample1.txt i wish to split in "Myfile", "MyArea1", "sample1", and "txt", becasue i need to use "sample1" label. I try to use "strsplit" but I am able just to split as "Myfile_MyArea1_sample1" and "txt" OR "Myfile", "MyArea1",

[R] help with Box plot

Dear researchers I wish to plot a box plot without the mean line (the black line) and the i wish a full line for the standard deviation This is an example mytest <- c(2.1,2.6,2.7,3.2,4.1,4.3,5.2,5.1,4.8,1.8,1.4,2.5,2.7,3.1,2.6,2.8) boxplot(mytest) really thanks Gianni [[alternative H

[R] Help boxplot to add mean, standard error and/or stadard deviation

Dear researchers I wish to plot a box plot without the mean line (the black line) and plot only the mean (red square). Futhermore, is it possible to add standard error and/or stadard deviation? This is an example mytest <- c(2.1,2.6,2.7,3.2,4.1,4.3,5.2,5.1,4.8,1.8,1.4,2.5,2.7,3.1,2.6,2.8) boxplo

[R] Help to improve bwplot plot (lattice)

Dear Researchers, I wish to plot mean, standard deviation, and standard error and I am using bwplot(). I have the following problems and sorry if maybe there are simple questions: 1- use color black for standarddevuation line and add horizontal end bar (as the commun graphic in scientific papers)

[R] plot with ylim with regural interval

Dear Researchers, sorry for the easy question but Is it possible to plot with an interval of 1 or .5 in a plot using ylim? Thanks gianni x = 0:10; y = 0:10; plot(x~y,ylim=c(0,10),las=1) [[alternative HTML version deleted]] __ R-help@r-projec

[R] R citation for the 2012

Dear Reasearchers, I am writing a report and i need (and wish) cite R. somebody know the citation of R for the 2012? or the more actual? thanks in advance Gianni [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https:/

Re: [R] R citation for the 2012

5:38 PM, Sarah Goslee wrote: > Typing > citation() > at an R prompt will provide you with complete citation information for > the version of > R you are using. Same goes for packages, with citation("pkgname"). > > Sarah > > On Wed, Feb 15, 2012 at 11:29 AM, gianni

[R] count how many row i have in a txt file in a directory

Dear Researchers, I have a large TXT (X,Y,MyValue) file in a directory and I wish to import row by row the txt in a loop to save only the data they are inside a buffer (using inside.owin of spatstat) and delete the rest. The first step before to create a loop row-by-row is to know how many rows th

[R] Problem to resolve a step for reading a large TXT and split in several file

Dear Researchs, It's the first time I am trying to resolve this problem. I have a TXT file with 1408452 rows. I wish to split file-by-file where each file has 1,000,000 rows with the following procedure: # split in two file one with 1,000,000 of rows and one with 408,452 of rows file <- "09G001_

[R] how disable the Error massage in read.table() " no lines available in input"

Dear Researchers, I am looking a way to disable the Error massage in read.table() as warn = TRUE in readLines(), when the lines are empty Error in read.table(con, header = F, sep = " ", nrow = n) : no lines available in input thanks for all suggestions Gianni [[alternative HTML versio

[R] question how to add Standard Deviation as "Whiskers" in a simple plot

Dear Researchers, sorry for this simple question. I have a point plot with mean values and i wish to plot line with Standard Deviation as "Whiskers". I calculate the mean+sd and mean-sd, but i can not figure out the way to add the line. mydata <- data.frame(mean=c(0.42,0.41,0.41,0.43,0.45,0.43,0

Re: [R] question how to add Standard Deviation as "Whiskers" in a simple plot

The function i am looking is a bars from the mean points of the plot in boxplot style. I tryed several forum but I have no clear the way to create these bars. Gianni On Mon, May 28, 2012 at 7:13 PM, David Winsemius wrote: > > On May 28, 2012, at 9:55 AM, gianni lavaredo wrote: >