similar to: (robust) mixed-effects model with covariate

Dear all, I am having quite a hard time in trying to figure out how to correctly spell out a model in R (a repeated-measures anova with a within-subject covariate, I guess). Even though I have read in the posting guide that statistical advice may or may not get an answer on this list, I decided to try it anyway, hoping not to incur in somebody's ire for misusing the tool. For the sake of

Dear all I would like to have some data in the form of a 2-dimensional array (matrix) of lists, so that I can easily find the desired list object by indexing the structure by rows and columns. In matlab there exists a data type called "cell array": a matrix of "cells", which are composite objects very similar to R lists. I know that in R you can create 1-dimensional

Hi all I am running a 2-way repeated measure anova with 1 between-subjects factor (Group=treatment, control), and 1 within-subject factor (Time of measurement: time1, time2). I extract the results of the anova with: summary(aov(effect ~ Group*Time + Error=Subj/Time, data=mydata)) Now, this must be clearly a dumb question, but how can I quickly extract in R all the post-hoc t-tests for the

Dear all, I am estimating a mixed-model in Ubuntu Raring (13.04¸ amd64), with the code: fm0 <- lme(rt ~ run + group * stim * cond, random=list( subj=pdSymm(~ 1 + run), subj=pdSymm(~ 0 + stim)), data=mydat1) When I check the approximate variance-covariance matrix, I get: > fm0$apVar [1] "Non-positive definite

Hi all I am running a 2-way repeated measure anova with 1 between-subjects factor (Group=treatment, control), and 1 within-subject factor (Time of measurement: time1, time2). I extract the results of the anova with: summary(aov(effect ~ Group*Time + Error=Subj/Time, data=mydata)) Now, this must be clearly a dumb question, but how can I quickly extract in R all the post-hoc t-tests for the

Dear R list I am looking for a function in R that computes the integration of a discrete curve, such as a power spectrum, in a specified interval (in my case, that would be 'power in a certain frequency band'). I found only functions, such as 'integrate', that perform adaptive quadrature on analytic functions, and not on a curve specified as a set of (x,y) pairs. I have the

Dear R users I have a really simple question (hoping for a really simple answer :-): Having estimated the spectral density of a time series "x" (heart rate data) with: x.pgram <- spectrum(x,method="pgram") I would like to compute the power in a specific energy band. Assuming that frequency(x)=4 (Hz), and that I am interested in the band between f1 and f2, is the

Hello, I'm new in this R world and I don't know much about statistics, but now I have to analize some data and I've got some first queries yet: I have 5 sets of area mesures and each set has 5 repetitions. My first step is to check data looking for outliers. I've used the outliers package. I have to use the cochran test and the grubbs test in case I find any outlier. The problem

I have a statistical question. The data sets I am working with are right-skewed so I have been plotting the log transformations of my data. I am using a Grubbs Test to detect outliers in the data, but I get different outcomes depending on whether I run the test on the original data or the log(data). Here is one of the problematic sets: fgf2p50=c(1.563,2.161,2.529,2.726,2.442,5.047)

Hi! I'm trying to implement an outlier test once/row in a large dataframe. Ideally, I'd do this then add the Pvalue results and the number flagged as an outlier as two new separate columns to the dataframe. Grubbs outlier test requires a vector and I'm confused how to make each row of my dataframe a vector, followed by doing a Grubbs test for each row containing the vector of numbers

Hello, I am trying to plot a few correlograms on the same figure, with the function corrgram() from the package corrgram. However, the function does not seem to use the base graphic system, as setting out the multiple figure layout with, e.g., par(mfrow=c(2, 2,)) does not work. Does anybody know a workaround for this? Many thanks in advance for any advice best giuseppe -- Giuseppe Pagnoni,

Dear All, I have small samples of data (between 6 and 15) for numerious time series points. I am assuming the data for each time 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

Happy holidays all! I know it's very subjective to determine whether some data is outlier or not... But are there reasonally good and realistic methods of identifying outliers in R? Thanks a lot! [[alternative HTML version deleted]]

Hello all, I would like to rigorously test whether observations in my dataset are outliers. I guess all the main tests in R (Grubbs) impose the assumption of normality. My data is surely not normal, so I would like to use something else. As far as I can tell from wikipedia, Peirce's criterion is just that. The data I am interested in testing is: 1) Continuous on the unit interval 2)

Hi, this is both a statistical and a R question... what would the best way / test to detect an outlier value among a series of 10 to 30 values ? for instance if we have the following dataset: 10,11,12,15,20,22,25,30,500 I d like to have a way to identify the last data as an outlier (only one direction). One way would be to calculate abs(mean - median) and if elevated (to what extent ?) delete the

Hi Dear All, Can someone give me a suggestion about which robust statistics are most appropriate for outlier detection in linear models, and is available with R ? Thanks for any idea. -- View this message in context: http://r.789695.n4.nabble.com/Robust-Statistics-for-Outlier-Detection-tp3438493p3438493.html Sent from the R help mailing list archive at Nabble.com.

I have been learning about some outlier tests -- Dixon and Grubb, specifically -- for small data sets. When I try help.start() and search for outlier tests, the only response I manage to find is the Bonferroni test avaiable from the CAR package... are there any other packages the offer outlier tests? Are the Dixon and Grubb tests "good" for small samples or are others more

Hi list, I want to write a general function so that it would take an lm object, extract its data element, then use the data at another R function (eg, glm). I searched R-help list, and found this would do the trick of the first part: a.lm$call$data this would return a name object but could not be recognized as a data.frameby glm. I also tried call(as.character(a.lm$call$data)) or

Dear R-help, I am trying to set up linear mixed effects models in R using the (recommended) nlme package (R version 2.3.1 on a Linux platform). When trying to reproduce an example from Jose Pinheiro & Douglas Bates (2000, p 210) I get the following error message (code to produce message pasted as well): library("nlme") data("Orthodont") vf1Ident <- varIdent(

Hi, there: As you all know, correlation is not a very robust procedure. Sometimes correlation could be driven by a few outliers. There are a few ways to improve the robustness of correlation (pearson correlation), either by outlier removal procedure, or resampling technique. I am wondering if there is any R package or R code that have incorporated outlier removal or resampling procedure in