search for: thatmike

...#make sure id, iv1, and iv2 are factorized a$id=factor(a$id) a$iv1=factor(a$iv1) a$iv2=factor(a$iv2) #run ANCOVA covariate_aov = aov(dv~covariate+group*iv1*iv2+Error(id/(iv1*iv2)),data=a) summary(covariate_aov) -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Check my public calendar: http://www.thatmike.com/mikes-public-calendar ~ Certainty is folly... I think. ~

On my machine (Mac OS 10.5.6, R 2.8.1) the following plot is drawn such that a substantial portion of the descender of the title is covered by the plot: library(ggplot2) qplot(x=1:10,y=1:10,main='p q j g') -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Do so at: http://www.timetomeet.info/with/mike/ ~ Certainty is folly... I think. ~ [[alternative HTML version deleted]]

...ether to collapse across repetition before running the stats, particularly since errors will leave unequal numbers of observations per cell if it's left in. Any advice anyone can provide would be great! Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Check my public calendar: http://www.thatmike.com/mikes-public-calendar ~ Certainty is folly... I think. ~

...p; google). Can anyone direct me to resources on using R for computational neuroscience? Input on my possibly naive assumption that R is a sufficient tool for this field would also be appreciated. Cheers, Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Check my public calendar: http://www.thatmike.com/mikes-public-calendar ~ Certainty is folly... I think. ~

...is. I am not yet savvy enough with 'R' to create a package, but if anyone out there is curious, let me know & I can supply the code. Thanks! Mike On Thu, Dec 18, 2008 at 4:21 AM, Mike Lawrence <mike@thatmike.com> wrote: > On Mon, Dec 15, 2008 at 9:48 PM, Mike Williamson <this.is.mvw@gmail.com>wrote: > >> I want to make some simple CDF (cumulative distribution function) plots >> to check whether distributions are Gaussian / normal. > > > qqnorm() > > >&gt...

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, 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:

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

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

...,col='red',lty=2) # It's clear that it's inappropriate to use the standard "glm(y~x,family='binomial')" # method to fit the modified logistic data, but what is the alternative? -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Check my public calendar and find a time that is mutually free: http://tinyurl.com/mikescalendar (tiny url redirects to google calendar) ~ Certainty is folly... I think. ~ [[alternative HTML version deleted]]

Dear List, One persistent feedback I am getting to people who are newly introduced to R ( especially in this cost cutting recession) is - 1) The website looks a bit old. While the current website does have a lot of hard work behind it, should n't a world class statistics package have a better website instead. You can check out www.knime.org which is an open source software , and free, and

Hello everyone I'm working with R 2.8.1 on a windows machine I have a question regarding time series analysis The first question is how does R expect the input file to be structured? I'm working with a *.txt file similar to the abbreviated one here: Date,stage 4/2/1953,7.56 4/3/1953,7.56 4/4/1953,7.54 4/5/1953,7.53 4/6/1953,7.5 4/7/1953,7.47 4/8/1953,7.44 4/9/1953,7.41 4/10/1953,7.37

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

...a, b=b, sub=sub), mean) ss.aBYb.cells = levs.d*sum((aBYb.cells$x-grand.mean)^2) ss.aBYb.error = ss.aBYb.cells - ss.aBYb - ss.subs round(ss.aBYb.error, 10)==round(summary(this_aov)[[5]][[1]]$Sum[2], 10) #not ok :( -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com