similar to: ggplot2's qplot() not rendering title descender

Hi all, I have data from an experiment with 3 independent variables, 2 are within and 1 is between. In addition to the dependent variable, I have a covariate that is a single measure per subject. Below I provide an example generated data set and my approach to implementing the ANCOVA. However the output confuses me; why does the covariate only appear in the first strata? Presumably it should

All, Thanks for all of your help & advice! I created a plot starting with qqnorm, if I remember right, but I had to add a LOT of extras. I was surprised that it doesn't include the probability (e.g. 50% at mean) and instead provides the standard deviations (e.g., 0 at mean) on the axis. Plus, there is no line through Q1 & Q3 to help give a sense of the normality of the graph.

Hi guRus, I'm looking for advice on a good way to approach analysis of some multi-level data I've obtained. I had humans classify words and measured response time. Words were 10 positive words ("happy", "joy", etc) and 10 negative words ("sad","grumpy", etc). Words were also presented in either white, red or green color. All variables were

Hi all, I've noticed that many computational neuroscience research groups use MATLAB. While it's possible that MATLAB may have some features unavailable in R, I suspect that this may instead simply be a case of costly tradition, where researchers were taught MATLAB as students and pay for it as researchers because it's all they know. I'd like to attempt to break the cycle by

This question might seem silly, because I felt that it MUST be in the mailing list archives or help files somewhere, but I simply couldn't find it. I want to make some simple CDF (cumulative distribution function) plots to check whether distributions are Gaussian / normal. But in order to check how "normal" the distribution is, I really need the y-axis to be Gaussian as well

I all, is there a way with R to perform an exact permutation test to replace the wilcoxon test to compare paired series and/or to perform pairwise multiple comparisons for related series after a Friedman test ? Thanks Gilles

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Hi all, Where f(x) is a logistic function, I have data that follow: g(x) = f(x)*.5 + .5 How would you suggest I modify the standard glm(..., family='binomial') function to fit this? Here's an example of a clearly ill-advised attempt to simply use the standard glm(..., family='binomial') approach: ######## # First generate some data ######## #define the scale and location of

Dear R users, I have a matrix "final" which looks like this: final oSO4 oNO3 mSO4 mNO3 [1,] 3.3728 0.2110 1.9517421 1.01883602 [2,] 0.8249 0.0697 1.5970292 0.11368781 [3,] 0.2636 0.1004 0.6012445 0.24356332 [4,] 8.0072 0.3443 6.1016998 3.63207149 [5,] 13.5079 0.6593 12.4011068 1.55323386 [6,] 6.1293 0.1989 5.7620926 0.12884845 [7,] 0.6004 0.0661

Dear all R users, I am going to use R to process some of my physiological data about eye. The problem is the recording machine does not sample in a reliably constant rate: the time intervals between data sampled can vary from 9msec to ~120msec, while most around in the 15-30msec range. The below is a fraction of a single data file of a trial: Time CursorX CursorY Pupilsize 1811543 -1 -1 -1

Hi, when performing a spearman_test stratified by a given factor in package "coin", how is it possible to obtain the value of Rho, the Spearman correlation coefficient ? Thanks in advance Gilles (F) [[alternative HTML version deleted]]

<> Hello One sometimes (quite often really ) marvel at the choice of defaults in R's graphic engine. For some obscure reason, mtext uses the typographical descender (bottom of letters) to align text. That is: "gG" will end up slightly higher that "GG". Depending on the font, "Q" might end up higher than "O". (for explanation of baseline

Hello, I have a very simple data set i imported from excel including 96 averages in a column along with 96 standard errors associated with those averages (calculated in excel). I plotted the 95 averages using r and I am wondering if it is possible to plot the second column of standard errors while applying error bars to each value so they represent the error corresponding to each average? thanks,

Hi!!! First time 'R' user looking for a little assistance. Here is what I have so far: practice1 = matrix ((runif(5000, min =0, max = 12)), 100) which is creating 50 samples, for 100 cases, distributed between 0-12. What I would like is to be able to set the mean and SD so that the data is normally distributed around lets say 7. Any help I can get with achieving that goal would be

To understand some functions about spatial point patterns in "spatstat" ,I should know some background about it, and the best way is to read the monograph, and "Statistical Analysis of Spatial Point Patterns" (2nd edt.) is a better choise. But I can not find it anywhere I can. Who can help me? Thank you! ----- unangu at gmail.com -- View this message in context:

The aggregate function does "almost" all that I need to summarize a datasets, except that I can't specify exclusion of NAs without a little bit of hassle. > set.seed(143) > m <- data.frame(A=sample(LETTERS[1:5], 20, T), B=sample(LETTERS[1:10], 20, T), C=sample(c(NA, 1:4), 20, T), D=sample(c(NA,1:4), 20, T)) > m A B C D 1 E I 1 NA 2 A C NA NA 3 D I NA 3 4 C I

I'm trying to create a colored map that would show the number of students per state. My data frame consists of two columns - state and count. I'm using the following code library(maps) map("usa") library(plotrix) state.col<-color.scale(gre$count,0,0,c(0,1)) map("state",fill=TRUE,col=state.col) I'm getting a map, but the values are not being mapped to correct

Hopefully this is a pretty simple question:   Is there a function in R that calculates the "mode" of a sample?   That is, I would like to be able to determine the value that occurs the most frequently in a data set.   I tried the default R "mode" function, but it appears to provide a storage type or something else.    I tried the RSeek and some R documentation that I

Hi all, I am trying to run a two factor anova, but one of the factors is a random factor, now I am also running in SPSS and it seems its dividing by the wrong term to get the appropriate F term. here is my data. In SPSS the F scores about double the ones in R, how can I specify one of my factors as a random factor or change it to where it does the right model fitting? I am using the lm command

Hello, here I have some code, which behaves quite strange, IMHO. There is a "res.txt" which will be collected before printing it. There is a paste-statement, which has a comment at the end of the line, which is this one: "# !!!!!HERE!!!!!". If you throw out the first hash-mark "#" on that line, the printout behaves like if it is done for each of the for-loop