similar to: how to sort dataframe levels

Hi there, I have some simple frequencies I want to plot into one graph. I had it working, and now I can't figure out whats going wrong. All the data is stored in a dataframe, and i finally managed to order the factor correctly! Each column is a variable and contains integers for the same set of values in the column that contains the headers for each row (graphLabels). So, I get the data

Hi all, I've uploaded STAR (Spike Train Analysis with R) on CRAN two days ago. The package is designed to analyze neuronal spike (action potential) trains. It uses S3 classes and methods and makes heavy use of other CRAN packages like gss, R2HTML, mgcv, survival. * Analysis of both spontaneous and stimulus evoked activity is implemented for single neuron spike trains as well as for many

Hi all, I've uploaded STAR (Spike Train Analysis with R) on CRAN two days ago. The package is designed to analyze neuronal spike (action potential) trains. It uses S3 classes and methods and makes heavy use of other CRAN packages like gss, R2HTML, mgcv, survival. * Analysis of both spontaneous and stimulus evoked activity is implemented for single neuron spike trains as well as for many

Dear all, I have 3 models (from simple to complex) and I want to compare them in order to see if they fit equally well or not. From the R prompt I am not able to see where I can get this information. Let´s do an example: fit1<- lm(response ~ stimulus + condition + stimulus:condition, data=scrd) #EQUIVALE A lm(response ~ stimulus*condition, data=scrd) fit2<- lm(response ~ stimulus +

Hello everybody, I need an help because I don´t know if the command for the ANOVA analysis I am performing in R is correct. Indeed using the function aov I get the following error:"In aov (......) Error() model is singular" The structure of my table is the following: subject, stimulus, condition, sex, response Example: subject stimulus condition sex response

Hi all, I dump the result of a long simulation, which I know has run to completion and returned a fairly complex list of lists of things. However, when I try to read it back in, I get the following error: > source ('../SIMULATIONS/run.1.R') Error in structure(list(do.extinction.runs, alpha = 1.8, delta.change = 0.005, : Object "do.extinction.runs" not found I think the

Dear all, I have tried to use R for the repeated measures experiment design in a phonetic study, and in the resulting forms I met some problems which perturb the final interpretation of the results. I will explain the experiment design first: the same 7 subjects were answering a question about 25 linguistic stimuli; the stimuli were the same utterances which were processed in 3 different ways (3

Dear all, I have got two vectors coding for a stimulus presented in the current trial (mydat$Stimulus) and a prediction in the same trial (mydat$Prediciton), respectively. By applying an if-conditional I want to create a new vector that indicates if there is a match between both vectors in the same trial. That is, if the prediction equals the stimulus. When I pick out some trials randomly, I get

Dear All, I have a problem in understanding how the interactions of 2 ways ANOVA work, because I get conflicting results from a t-test and an anova. For most of you my problem is very simple I am sure. I need an help with an example, looking at one table I am analyzing. The table is in attachment and can be imported in R by means of this command: scrd<-

Dear all, I need an help because I don´t know how to perform the analysis in the right way, as I get different beheaviors using t-test and two ways ANOVA. In what follow I post the table, my goal and the strange results I got. I kindly ask you an help because I really don´t know how to solve this problem. So the table is this: number stimulus condition response 1

Dear all, I need to understand a thing in the beheaviour of the two functions aov and anova in the following case involving an analysis of ANOVA with repeated measures: If I use the folowing command I don´t get any problem: >aov1 = aov(response ~ stimulus*condition + Error(subject/(stimulus*condition)), >data=scrd) > summary(aov1) Instead if I try to fit the same model for the

Dear all, how can I perform a post hoc analysis for ANOVA with repeated measures (in presence of a balanced design)? I am not able to find a good example over internet in R...is there among you someone so kind to give me an hint with a R example please? For example, the aov result of my analysis says that there is a statistical difference between stimuli (there are 7 different stimuli). ...I

Dear List, I am familier with binary models, however i am now trying to get predictions from a ordinal model and have a question. I have a data set made up of 12 categorical predictors, the response variable is classed as 1,2,3,4,5,6, this relates to threat level of the species ( on the IUCN rating). Previously i have combined levels 1 and 2 to form = non threatened and then combined 3-6 to

Dear Colleagues, I hope to know how ordinal logistic regression with a mixed model is made in R. We (My colleague and I) are studying the behavior of a beetle. The attraction of beetles to a stimulus are recorded: the response is Slow, Mid, or Fast. They are based on the time after the presentation of the stimulus to the beetles. Because we do not observe the behavior continuously but do

Dear List, This is my first post, and I'm a relatively new R user trying to work out a mixed effects model using lme() with random effects, and a correlation structure, and have looked over the archives, & R help on lme, corClasses, & etc extensively for clues. My programming experience is minimal (1 semester of C). My mentor, who has much more programming experience, but a comparable

Hi, I am analyzing a repeated-measures Greco-Latin Square with the aov command. I am using aov to calculate the MSs and then picking by hand the appropriate neumerator and denominator terms for the F tests. The data are the following: responseFinger mapping.code Subject.n index middle ring little ---------------------------------------------------------------------------- 1 1

mgcv: Constructing probabilities for binomial GAM with crossed random intercepts and factor by variable Hello, (I'm sorry if this has been discussed elsewhere; I may not have been looking in the right places.) I ran a binomial GAM in which "Correct" is modelled in terms of the participant's age and the modality in which the stimulus is presented (written vs spoken).

Hi, I am playing with clogit(case~spontaneous+induced+strata(stratum),data=infert) from clogit help file. This line works. 1. But, why strata(stratum) doesn't have a coefficient like spontaneous and induced? 2. When I remove strata(stratum) from the command, this function seems to keep running forever. Why? 3. I think the equation for clogit looks like P=1/(1+

Dear list member, I deperately need an help in performing a MANOVA in R, but I encountered some problems both in the design and in the synthax with R. I conducted a listening experiment in which 16 participants had to rate the audio stimuli along 5 scales representing an emotion (sad, tender, neutral, happy and aggressive). Each audio stimulus was synthesized in order to represent a

Dear UseRs! We are pleased to announce the release of our new package 'PVAClone'. The 'PVAClone' package implements Population Viability Analysis (PVA) methodology using data cloning. The data cloning algorithm by Lele et al. (2007, 2010) is employed to compute maximum likelihood estimates of the state-space model parameters and the corresponding standard errors, heavily