search for: oulier

...e point is normally distributed. The problem is that the data arrvies sporadically and I would like to detect the number of outliers after I have six data points for any time period. Essentially, I would like to detect the number of outliers when I have 6 data points then test whether there are any ouliers. If so, remove the outliers, and wait until I have at least 6 data points or when the sample size increases and test again whether there are any outliers. This process is repeated until there are no more data points to add to the sample. Is it valid to use the grubbs.test in this way? If not, ar...

...tipo. Tienes bibiligrafía de artículos cientificos en las que se quitan valores siguiendo los criterios que te he dicho solo has de buscar bibliografía del campo en el que te mueces. Osea que nos quitarlos por que sí, aunque en cada area de ciencia tienes una serie de criterios para acabar con los ouliers. Yo personalmente me fio más de mi buena intuición(y validación) que de un test. De todas maneras seguro que tienes un informe de como recogiste los datos y puedes trazar el origen de este dato en concreto. Creo que eso es mejor que cualquier test. De todas maneras a mi el test me sale: > setw...

...hen there is an outlier IO in the model. Hi R users I have a problem when a use the predict() method in an object of class arimax ( These objects are the results of the implementation of the function arimax() from the TSA library) . The object is a model of a time series in which I identified an IO oulier at the element 33 of the serie (using the function detectIO() from the package TSA in an object arimax that corresponds a model ARIMA(1,1,0) that not include the outlier) , so I incorporated this in the model ARIMA(1,1,0): mz3B<-arimax(z3B,io=c(33),order=c(1,1,0)) where the object z3B is a time...

Dear List, I just googled to find out if notched box plots are possible with ggplot2, but couldn't find a answer to it. boxplot() has the option: notch = TRUE, e.g.: boxplot(mpg$hwy, notch=TRUE) My example code (taken from the net) is: require(ggplot2) qplot(class, hwy, fill=factor(year), data=mpg, geom="boxplot", position="dodge")+theme_bw() Thank you for you help!

To enhance my understanding, and that of my students, I have a question about cox.zph in the survival package. If I have correctly gleaned the high-level point from the 1994 Biometrika paper of Grambsch and Therneau, it looks to me like cox.zph provides a mechanism to test for a simple trend in plots of a function of time, g(t) versus the scaled schoenfeld residuals and it also provides some

Tu fichero tiene los decimales como puntos y no como comas como tu le indicas. Te dejo un ejemplo #--------------------------------------------------------------------------------------------------------------------- setwd(dir="c:/Users/usuario/Desktop/") library(outliers) filename<-"timediff.csv" time<-read.csv(filename, sep=";",header=TRUE,dec=".")