Re: [R] Help optimizing EMD::extrema()

2011-02-13 Thread Mike Lawrence
l Dunlap > Spotfire, TIBCO Software > wdunlap tibco.com > >> -Original Message- >> From: r-help-boun...@r-project.org >> [mailto:r-help-boun...@r-project.org] On Behalf Of William Dunlap >> Sent: Sunday, February 13, 2011 2:08 PM >> To: Mike Lawrence; r-h..

[R] [R-pkgs] ez version 3.0

2011-02-11 Thread Mike Lawrence
()" about which I wish I could have warned users sooner. To facilitate future rapid notification of users, I've created a discussion group (http://groups.google.com/group/ez4r) to which users can/should subscribe to keep up to date on development news. Additionally, I created a project page on

[R] Help optimizing EMD::extrema()

2011-02-11 Thread Mike Lawrence
Hi folks, I'm attempting to use the EMD package to analyze some neuroimaging data (timeseries with 64 channels sampled across 1 million time points within each of 20 people). I found that processing a single channel of data using EMD::emd() took about 8 hours. Exploration using Rprof() suggested t

[R] hierarchical mixture modelling

2010-11-22 Thread Mike Lawrence
greatly appreciated! Cheers, Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is folly... I think. ~ __ R-he

[R] [R-pkgs] ez version 2.0

2010-09-01 Thread Mike Lawrence
w of the plot facets. - see changelog for further changes Enjoy! Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is folly...

[R] compare gam fits

2010-08-05 Thread Mike Lawrence
not independent from those next to it make this an inappropriate concern? Cheers, Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calend

Re: [R] Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)

2009-11-10 Thread Mike Lawrence
nt') , variable_name='day' ) b$subject=paste(b$subject,b$treatment) #create unique sids ezANOVA( data=b , dv=.(value) , sid=.(subject) , within=.(day) , between=.(treatment) ) On Tue, Nov 10, 2009 at 6:47 AM, Mike Lawrence wrote: > Check out the reshape package for transform

Re: [R] Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)

2009-11-10 Thread Mike Lawrence
idata, idesign=~Day)) >> >> Peter Dalgaard also pointed toward an article that describes how to do the >> same thing with anova(). >> >> Regards, >>  John >> >>> -Original Message- >>> From: r-help-boun...@r-project.org [mailto:r-help-boun..

Re: [R] Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)

2009-11-09 Thread Mike Lawrence
ebraska-Lincoln > Lincoln, NE 68588-0308  USA > > > > > On Mon, Nov 9, 2009 at 2:25 PM, Mike Lawrence wrote: >> >> Have you tried ezANOVA from the ez pacakge? It attempts to provide a >> simple user interface to car's ANOVA (and when that fails, aov). &

Re: [R] Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)

2009-11-09 Thread Mike Lawrence
lp@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > -- Mike Lawrence Graduate Student Department of Psycho

Re: [R] Quickly generate all possible combinations of 2 groups

2009-11-09 Thread Mike Lawrence
In case anyone is interested in the random variant, here's some code: http://gist.github.com/229929 -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is

[R] Quickly generate all possible combinations of 2 groups

2009-11-09 Thread Mike Lawrence
the list to see if there's a computationally efficient solution to the "test all combinations" scheme. Cheers, Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calenda

[R] EM algorithm to fit circular mix of uniform+Von Mises

2009-11-07 Thread Mike Lawrence
n Mises. I couldn't find anything in the contributed packages that seemed to suit this purpose. Any pointers would be greatly appreciated! Cheers, Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calenda

Re: [R] Creating a specific skewed distribution

2009-06-10 Thread Mike Lawrence
do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a me

Re: [R] graphically representing frequency of words in a speech?

2009-06-07 Thread Mike Lawrence
my_opts+ coord_polar() ) -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is folly... I think. ~ __ R-help

Re: [R] Fitting a Weibull Distribution

2009-06-05 Thread Mike Lawrence
ect.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalho

Re: [R] What do you think about my function?

2009-06-02 Thread Mike Lawrence
this message in context: > http://www.nabble.com/What-do-you-think-about-my-function--tp23831687p23831687.html > Sent from the R help mailing list archive at Nabble.com. > > ______ > R-help@r-project.org mailing list > https://stat.ethz

Re: [R] help on understanding a code

2009-06-02 Thread Mike Lawrence
> R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >

Re: [R] Most used R editors

2009-06-02 Thread Mike Lawrence
ject.org/rghuang.html > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained

Re: [R] Within Subject ANOVA question

2009-06-01 Thread Mike Lawrence
such as in t.test: >>as.numeric(t.test(serie1,serie2,na.rm=TRUE)$statistic) > > What method should I use to extract, or is it simpler to write compute > the figures that I need directly in this case? > Thanks a lot! > > - John > > On Mon, Jun 1, 2009 at 7:08 PM, Mik

Re: [R] Within Subject ANOVA question

2009-06-01 Thread Mike Lawrence
gful, something you surely do not want to do with your DV in this example. -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is folly... I think. ~ _

Re: [R] Within Subject ANOVA question

2009-06-01 Thread Mike Lawrence
subtle for me, please point out the errors that I have made if you > can. > Thank you very much! > > - John > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listin

[R] Anyone using Bespin & R?

2009-06-01 Thread Mike Lawrence
Just curious whether anyone has any experiences to relate with regards to using Bespin (https://bespin.mozilla.com/) as a text editor for your R scripting. Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public

Re: [R] ggplot2 legend

2009-05-28 Thread Mike Lawrence
1", "152", "153", "154", > "155", "156", "157", "158", "159", "160", "161", "162", "163", > "164", "165", "166

Re: [R] Labeling barplot bars by multiple factors

2009-05-27 Thread Mike Lawrence
rg mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department

Re: [R] Multiple ANOVA tests

2009-05-27 Thread Mike Lawrence
 18.8 > --- > Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > 246 observations deleted due to missingness > >> summary(q) >          Length Class       Mode > X11_20502 1      summary.aov list > X11_21067 1      summary.aov list > X11_10419

Re: [R] Multiple ANOVA tests

2009-05-27 Thread Mike Lawrence
23739615.html > Sent from the R help mailing list archive at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html &

Re: [R] Help me...!!!

2009-05-26 Thread Mike Lawrence
list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of

Re: [R] help wth boxplots

2009-05-26 Thread Mike Lawrence
; > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >

Re: [R] How to create all pairs

2009-05-25 Thread Mike Lawrence
LEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check

Re: [R] help with replacing factors

2009-05-24 Thread Mike Lawrence
dvance > > Andreas > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide htt

Re: [R] Coord_equal in ggplot2

2009-05-19 Thread Mike Lawrence
cial position of INBO, as long as the message is not confirmed by a > duly > signed document. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-proje

Re: [R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Mike Lawrence
posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is fol

Re: [R] sample variance from simulation

2009-05-18 Thread Mike Lawrence
cessary corrections, i'd be very very greatful. > > Thanks heaps guys for taking the time to look at this > >> Date: Mon, 18 May 2009 15:06:47 +0200 >> From: waclaw.marcin.kusnierc...@idi.ntnu.no >> To: konk2...@hotmail.com >> CC: mike.lawre...@dal.ca; r-help@r-

Re: [R] sample variance from simulation

2009-05-18 Thread Mike Lawrence
b.aspx?URL=http%3A%2F%2Fsecure%2Dau%2Eimrworldwide%2Ecom%2Fcgi%2Dbin%2Fa%2Fci%5F450304%2Fet%5F2%2Fcg%5F801459%2Fpi%5F1004813%2Fai%5F859641&_t=762955845&_r=tig_OCT07&_m=EXT >[[alternative HTML version deleted]] > > > __ >

Re: [R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Mike Lawrence
__ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, rep

Re: [R] Simulation

2009-05-13 Thread Mike Lawrence
ative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal,

Re: [R] read multiple large files into one dataframe

2009-05-13 Thread Mike Lawrence
* > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide h

Re: [R] Two-way Anova

2009-05-12 Thread Mike Lawrence
o read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check

Re: [R] Kumaraswamy distribution

2009-05-12 Thread Mike Lawrence
g list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Look

Re: [R] Bumps chart in R

2009-05-06 Thread Mike Lawrence
gt; On Tue, Apr 28, 2009 at 12:09 AM, Andreas Christoffersen > wrote: >> >> My legend is removed! - Couldn't find it in your ggplot2 book - but >> here it is. Brilliant - thank you very much. >> > > __ > R-help@r-proje

Re: [R] General Question about Graphics

2009-05-06 Thread Mike Lawrence
//www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certaint

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Mike Lawrence
e is great for computing summary statistics within cells of the design. Mike *note: I typically use Python for data collection (showing visual stimuli, recording responses, etc), but have it spit out raw text files of the trial-by-trial data, and thus use it for only a bare minimum of processing. -- M

Re: [R] novice question regarding R archives

2009-05-05 Thread Mike Lawrence
g-guide.html > and provide commented, minimal, self-contained, reproducible code. > > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certainty is folly... I think.

Re: [R] plot estimates and their 95% confidence intervals

2009-04-27 Thread Mike Lawrence
_ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- M

Re: [R] Bumps chart in R

2009-04-26 Thread Mike Lawrence
, y = 'Rank\n' ) #add the lines r = r + geom_line() #add the text r = r + geom_text( aes( label=countries , hjust = hjust ) ) #expand the axis to fit the text r = r + scale_x_discrete( expand=c(2,2) ) #s

Re: [R] boxplot of two variables

2009-04-23 Thread Mike Lawrence
://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tr.im/mikes_public_calendar ~ Certain

[R] OT: Strange bootstrap approach to significance testing

2009-04-21 Thread Mike Lawrence
f group. In my head this makes some sense (though I'd think a straight permutation test would be a lot simpler), but having never heard of anything like this before I thought I'd see what others on this list think of the approach. Mike -- Mike Lawrence Graduate Student Department of P

Re: [R] Sampling in R

2009-04-21 Thread Mike Lawrence
Dsk%3Atl%3Asknzsal%3Amsnnz%3A0%3Ahottag%3Aearn%5Fenough&_t=757263783&_r=Seek_NZ_tagline&_m=EXT >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help &

Re: [R] Paired test for 2-way ANOVA

2009-04-20 Thread Mike Lawrence
_ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Dep

Re: [R] Welcome to the "R-help" mailing list

2009-04-18 Thread Mike Lawrence
gt; Nothing happens. > > I appreciate your advice on what is going wrong? Thanks. > > Best, > Tony > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide

Re: [R] F test

2009-04-16 Thread Mike Lawrence
> regression linear model lm ? >> > >> > I appreciate your help >> > >> > >> > -- >> > View this message in context: >> > http://www.nabble.com/F-test-tp23078122p23078122.html >> > Sent from the R help mailing list archive at Nabble.com. >> > &

Re: [R] F test

2009-04-16 Thread Mike Lawrence
ed effect > and a random effect)? and where to find the result of the interaction? > Thanks. > > Jun > > On Thu, Apr 16, 2009 at 12:08 PM, Mike Lawrence > wrote: >> >> summary(my_lm) will give you t-values, anova(my_lm) will give you >> (equivalent) F-values. summ

Re: [R] F test

2009-04-16 Thread Mike Lawrence
com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducibl

Re: [R] performing function on data frame

2009-04-16 Thread Mike Lawrence
___ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike La

Re: [R] How to apply a function to all combinations of the values of 2 or more variables?

2009-04-15 Thread Mike Lawrence
é(s) > est interdite. > Si vous recevez ce courrier électronique par erreur, veuillez m'en aviser > immédiatement, par retour de courrier électronique ou par un autre moyen. > >        [[alternative HTML version deleted]] > > >

Re: [R] cbind

2009-04-14 Thread Mike Lawrence
4  67 >> [121,]  45  52 >> [122,]  46  86 >> [123,]  47  79 >> [124,]  48  97 >> [125,]  49  80 >> [126,]  50 101 >> [127,]  51  87 >> [128,]  52  53 >> [129,]  53  85 >> [130,]  54  79 >> [131,]  55  86 >> [132,]  56 104 >> [133,]  57 153 >> [134,]  5

Re: [R] cbind

2009-04-14 Thread Mike Lawrence
oops, of course I meant: a=c(NA,NA,2,3,NA,NA,NA) b=1:length(a) cbind(b=b-which.min(is.na(a)),a=a) On Tue, Apr 14, 2009 at 9:43 AM, Mike Lawrence wrote: >> a=c(NA,NA,2,3,NA,NA,NA) >> which.min(is.na(a)) > [1] 3 >> b=1:length(a) >> b-3 > [1] -2 -1  0  1  2  3  4 &g

Re: [R] cbind

2009-04-14 Thread Mike Lawrence
,]  72  NA > [149,]  73  NA > [150,]  74  NA > [151,]  75  NA > > could anyone help me to with a function that would be able to calculate the > sequence I require to bind to the initial sequence? > > thanks emma > -- > View this message in context: > http://www.nabble

Re: [R] weighted mean and by() with two index

2009-04-13 Thread Mike Lawrence
> Donghyun Oh > CESIS, KTH > -- > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org m

Re: [R] Group by in R

2009-04-13 Thread Mike Lawrence
_ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike

Re: [R] Quantative procedure assessing if data is normal

2009-04-12 Thread Mike Lawrence
> > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and p

Re: [R] "taking the log then later exponentiate the result" query

2009-04-12 Thread Mike Lawrence
t(temp$obs.acc.rate) > > Louisa > > > > Get the New Internet Explore 8 Optimised for MSN. Download Now > > _ > > >        [[alternative HTML version deleted]] > > _

Re: [R] Genstat into R - Randomisation test

2009-04-08 Thread Mike Lawrence
rical significance) and start using "exhaustive permutation test" (to contrast with "non-exhaustive permutation test"). On Wed, Apr 8, 2009 at 1:18 PM, Peter Dalgaard wrote: > Mike Lawrence wrote: >> >> Looks like that code implements a non-exhaustive var

Re: [R] Genstat into R - Randomisation test

2009-04-08 Thread Mike Lawrence
*sigh* and the elapsed time line should be: cat('Elapsed time:',proc.time()[1]-start.time) On Wed, Apr 8, 2009 at 11:35 AM, Mike Lawrence wrote: > *Hits head* > Of course, the approach taken by your Genstat code of only shuffling > one variable is sufficient and faster: &

Re: [R] Genstat into R - Randomisation test

2009-04-08 Thread Mike Lawrence
1:num.permutations){ IDorder = sample(index) perm.cor[i] = .Internal(cor(ID[IDorder], CD, 4, FALSE)) } cat('Elapsed time:',start.time-proc.time(1)) sum(perm.cor>obs.cor)/num.permutations On Wed, Apr 8, 2009 at 11:33 AM, Mike Lawrence wrote: > "Assuming you h

Re: [R] Genstat into R - Randomisation test

2009-04-08 Thread Mike Lawrence
"Assuming you have a data frame with columns ID & CD, this should do it" Oops, the code I sent doesn't assume this. It assumes that you have two vectors, ID & CD, as generated in the first 3 lines. On Wed, Apr 8, 2009 at 11:31 AM, Mike Lawrence wrote: > Looks like

Re: [R] Genstat into R - Randomisation test

2009-04-08 Thread Mike Lawrence
n > Altenbergrain 21 > CH-3103 Bern > kem...@ips.unibe.ch > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html &g

Re: [R] threshold distribution

2009-04-05 Thread Mike Lawrence
x27;lognormal') > pdf("adj.pdf") > plot(density(x)) > lines(dlnorm(x, -.1348, .1911), col='red') > points(x,dlnorm(x, -.1348, .1911), col='blue') > curve(dlnorm(x, -.1348, .1911), col='green',add=T) > dev.off() > > > > _

Re: [R] repeated measures ANOVA - among group differences

2009-04-01 Thread Mike Lawrence
Hm, it seems I possibly used the technical term "nested" inappropriately in my response. I meant: "If Month is a repeated measure within each Quadrat..." and "If Treatment is also a repeated measure within each Quadrat..." On Wed, Apr 1, 2009 at 8:21 AM, Mike

Re: [R] repeated measures ANOVA - among group differences

2009-04-01 Thread Mike Lawrence
    0.461538462 0.746898594 3     4 > 4     0     0     4     4 > 4     0.125 0.361367124 5     4 >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/li

Re: [R] help about: anova and population no normal

2009-03-31 Thread Mike Lawrence
simple to instead compute, say, the best fit weibull to each condition and simultaneously test shift, shape, and scale. On Tue, Mar 31, 2009 at 12:18 PM, Mike Lawrence wrote: > Those with more formal statistical backgrounds may provide better > advice, but in my own informal training I&#x

Re: [R] help about: anova and population no normal

2009-03-31 Thread Mike Lawrence
> Many thanks for your help with these questions. > > Maurizio > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and

Re: [R] R Language Bloggers and Web Sites Owners

2009-03-31 Thread Mike Lawrence
Hm, this feels reminiscent of the Wellsphere free blog content scam: http://neurocritic.blogspot.com/2009/02/wealthcentral.html -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tinyurl.com/mikes

Re: [R] Calculating First Occurance by a factor

2009-03-30 Thread Mike Lawrence
e in context: >> http://www.nabble.com/Calculating-First-Occurance-by-a-factor-tp22789964p22789964.html >> Sent from the R help mailing list archive at Nabble.com. >> >> __ >> R-help@r-project.org mailing list >> https

Re: [R] how to input multiple .txt files

2009-03-30 Thread Mike Lawrence
, > > > a <- do.call(rbind, lapply(cust1_files, read.table)) > > (i believe expanding objects in a for loop belong to the R Inferno) > > baptiste > > On 30 Mar 2009, at 12:58, Mike Lawrence wrote: > >> >> cust1_files = >> list.files(path=pa

Re: [R] how to input multiple .txt files

2009-03-30 Thread Mike Lawrence
cust2_files = list.files(path=path_to_my_files,pattern='cust2',full.names=TRUE) a=NULL for(this_file in cust2_files){ a=rbind(a,read.table(this_file)) } write.table(a,'cust2.master.txt') On Mon, Mar 30, 2009 at 8:55 AM, Mike Lawrence wrote: > my_files = list.files

Re: [R] how to input multiple .txt files

2009-03-30 Thread Mike Lawrence
elp > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Ch

Re: [R] loading and manipulating 10 data frames-simplified

2009-03-28 Thread Mike Lawrence
-- >> >> >> I have tried : >> >> Bin <- lapply(1:10, function(.file){ >> >> #Loads bin data frame from csv files with acres and TAZ data >> fileName >> <-paste("I:/Research/Samba/urb_transport_modeling/LUSDR/Workspace/BizLand

Re: [R] Way to rotate a histogram?

2009-03-26 Thread Mike Lawrence
,axis.ticks = theme_blank() ,axis.text.y = theme_blank() ,axis.title.y = theme_blank() ) print(p1,vp = viewport(width = 1/3,height = 1,x = 1/3*.5,y = .5)) print(p2,vp = viewport(width = 1/3,height = 1,x = 1/3+1/3*.5,y = .5)) print(p3,vp = viewport(width = 1/3,height = 1,x = 2/3+1/3*.5,y = .5)

Re: [R] geometric mean of probability density functions

2009-03-18 Thread Mike Lawrence
starship(sim.x) cat('Sim convergence:',sim.fit$optim.results$convergence,'\n') sim.gld.par = sim.fit$lambda #assess the fit using sum squared scaled error gld.par.sum.sq.sc.err = sum(((emp.gld.par.means$x-sim.gld.par)/sim.gld.par)^2) print(gld.par.sum.sq.sc.err) -- Mike Law

Re: [R] Way to rotate a histogram?

2009-03-17 Thread Mike Lawrence
t; > Perhaps with a red line through the hinge points. And probably need to get > the axes aligned proerly but it's very close. > > -- > david Winsemius > > > > On Mar 17, 2009, at 5:14 PM, Mike Lawrence wrote: > >> library(ggplot2) >> test_d

Re: [R] Way to rotate a histogram?

2009-03-17 Thread Mike Lawrence
posting guide >> > http://www.R-project.org/posting-guide.html >> > and provide commented, minimal, self-contained, >> > reproducible code. >> >> __ >> R-help@r-project.org mailing list >> https://stat.ethz.ch/mailman/

Re: [R] How to find maximum values on the density function of a random variable

2009-03-12 Thread Mike Lawrence
ml > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tinyurl.com/mikes-public-calendar ~ Certainty is folly... I think. ~ __

Re: [R] Centering multi-line strip text in lattice

2009-03-11 Thread Mike Lawrence
Perfect, thanks! On Wed, Mar 11, 2009 at 7:41 PM, Paul Murrell wrote: > Hi > > > Mike Lawrence wrote: >> I'm having trouble centering multi-line strip text in lattice. As the >> code below demonstrates bounding box of the text is centered within >> the strip, b

[R] OT: Likelihood ratio for the randomization/permutation test?

2009-03-11 Thread Mike Lawrence
observed difference score. Any guesses on how to translate this into a quantification of evidence? Mike -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tinyurl.com/mikes-public-calendar ~ Certainty is f

[R] Centering multi-line strip text in lattice

2009-03-10 Thread Mike Lawrence
m to do much. Suggestions? a=data.frame( x=rep(1:10,2) ,y=rep(1:10,2) ,z=rep(c('First Line\nLonger Second Line (1)','First Line\nLonger Second Line (2)'),each=10) ) xyplot( y~x|z ,data=a ,par.strip.text = list(cex = .75, lineheight=

Re: [R] Changing factor to numeric

2009-03-10 Thread Mike Lawrence
at.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange

Re: [R] simple question beginner

2009-03-10 Thread Mike Lawrence
ive HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, s

Re: [R] find the important inputs to the neural network model in "nnet" package

2009-03-10 Thread Mike Lawrence
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://ti

Re: [R] info

2009-03-09 Thread Mike Lawrence
__ > > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented

Re: [R] how to estimate distribution?

2009-03-09 Thread Mike Lawrence
> Simone > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, repr

Re: [R] Computing sd across an array with missing values

2009-02-26 Thread Mike Lawrence
)) >>>>> >>>>> the apply function takes a long time to run. >>>>> >>>>> Is there a faster, more efficient way to do this? >>>>> >>>>> Thanks in advance >>>>> >>>>> Matt >>&

Re: [R] colored maps again

2009-02-19 Thread Mike Lawrence
g mapped to correct states. > What do I need to do to fix that? > >[[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting

Re: [R] error bars

2009-02-18 Thread Mike Lawrence
osting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Check my

Re: [R] anova help

2009-02-15 Thread Mike Lawrence
-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Mike Lawrence Graduate Student Department of Psyc

Re: [R] Mixed ANCOVA with between-Ss covariate?

2009-02-12 Thread Mike Lawrence
I'm surprised there were no takers on this query; I thought it would be an easy answer, particularly where I provided example data set and code. Did my request run afoul of the list etiquette? On Tue, Feb 10, 2009 at 5:12 PM, Mike Lawrence wrote: > Hi all, > > I have data from an e

[R] Mixed ANCOVA with between-Ss covariate?

2009-02-10 Thread Mike Lawrence
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-pub

Re: [R] lme() direction

2009-02-07 Thread Mike Lawrence
y(lme( fixed = rt~valence*luminance , data = a ,random = ~1|id )) On Sat, Feb 7, 2009 at 12:34 PM, Dieter Menne wrote: > Mike Lawrence thatmike.com> writes: > >> >>Would it improve things if "type" were a continuous variable rather >>th

Re: [R] lme() direction

2009-02-07 Thread Mike Lawrence
Would it improve things if "type" were a continuous variable rather than categorical? I chose words at the extreme ends of a valence rating scale but I still have the raw valence ratings for each word. On Sat, Feb 7, 2009 at 12:02 PM, Dieter Menne wrote: > Mike Lawrence thatmik

  1   2   >