[R] Exact wilcoxon may differ in R and SPSS/StatXact (due to round off in the latter pair)

2008-01-16 Thread Nilsson Fredrik X
Dear R-users, If you use the exact Wilcoxon test in the coin package, I would like make you aware of that SPSS/StatXact MAY perform a round-off before doing their exact Wilcoxon-Mann-Whitney test (if you ever are unlucky enough not to use R). I have data from two treatments and was surprised to f

[R] How to do a time-stratified case-crossover analysis for air pollution data?

2008-03-07 Thread Nilsson Fredrik X
Dear Experts, I am trying to do a time-stratified case-crossover analysis on air pollution data and number of myocardial infarctions. In order to avoid model selection bias, I started with a simple simulation. I'm still not sure if my simulation is right. But the results I get from the "ts

[R] How to do a time-stratified case-crossover analysis for air pollution data? Unformatted text-version, with an additional note

2008-03-07 Thread Nilsson Fredrik X
Dear Experts, I am trying to do a time-stratified case-crossover analysis on air pollution data and number of myocardial infarctions. In order to avoid model selection bias, I started with a simple simulation. I'm still not sure if my simulation is right. But the results I get from the "ts-ca

Re: [R] [R-sig-ME] mixed model testing

2007-11-08 Thread Nilsson Fredrik X
Dear John, Forgive me for putting my nose out, I hope that I'm not rude, but I am a bit bewildered by your mail (and by statistical modelling). I agree that if your model is: Lawndepression~lawn.roller.weight +(1|lawn.id), When, in fact, it *should be* (because you simulated the data or you're

Re: [R] question about non-linear least squares in R

2007-09-11 Thread Nilsson Fredrik X
Dear Warren, I had similar problems, this is (roughly) how I solved it translated to your problem: x <- c(-0.06,-0.04,-0.025,-0.015,-0.005,0.005,0.015,0.025,0.04,0.06) y <- c(1866760,1457870,1314960,1250560,1184850,1144920, 1158850,1199910,1263850,1452520) dafa<-data.frame(x,y) plot(x,y)