search for: bmdp

I used both BMDP and SAS in my earlier years, side by side. At that time the BMDP statistical methods were much more mature and comprehensive: we treated them as the standard when the two packages disagreed. (It was a BMDP manual that clearly explained to me what the hypothesis of "Yate's weighted mean t...

[Sorry, In case this is repeating a message already sent to the list]. I am trying to move a project to R (base or nlme), for which I have a partial solution in BMDP 8V. Here is the 8V control language: /input title='Augenbewegungen'. variables=4. file='latm.dat'. format='11x,f7.0,f7.0,f12.4,f10.0'. /variables names= llat,rlat,vg,diff. /design ndep = 4. level=21,2,2,9,3. names = person,retest,direct,flicker,grou...

I am trying to translate a given BMDP 8V problem treatment to something approximately equivalent in nlme package (or base R, if sufficient). (Of course, it is not my goal to end there. I want to get access to the advanced flexibility of a more modern treatment.) Here is the problem statement in BMDP control language: /input title=...

Hi All, I've been fiddling around with various ways to estimate the popularity of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS, R-PLUS and Revolution R. It's not an easy task. You can see what I've come up with so far at http://r4stats.com/popularity . I'm sure people will have plenty of ideas on how to improve this, so please let me know what you think. Cheers, Bob ==================================...

I have design with two repeated-measures factor, and no grouping factor. I can analyze the dataset successfully in other software, including my legacy DOS version BMDP, and R's 'aov' function. I would like to use 'Anova' in 'car' in order to obtain the sphericity tests and the H-F corrected p-values. I do not believe the data are truly deficient in rank. I followed the methods for this kind of analysis outlined in Bennett'...

I wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

...having to suspend the > metaphor of a "libname as > folder" and think of the libname check as pointing > to a file: > "testing.xpt". It works that way in trying to read > spss files with the > spss engine too. I would assume it's also true for > using bmdp and osiris > engines, but I haven't had occasion to use them. > > > Hope this helps, > > Scot > > > > -- > Scot W. McNary email:smcnary at charm.net > >

...>--- >Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 > >Error: Within > Df Sum Sq Mean Sq F value Pr(>F) >Residuals 48 318 7 > This is mostly identical the analysis by BMDP 4V, except for the Greenhouse-Geisser epsilons, which are not estimated this way. I have to analyse a similar dataset, which is not balanced. So I need to change the method. Following Pinheiro/Bates p.90f, I tried > hf2.lme <- lme(ampl~gapf*bl,hframe2,random=list(VP=pdDiag(~gapf*bl),bl=pdDia...

Data from Fisher's paper: Confidence Limits for a Cross-Product Ratio. > y col1 col2 [1,] 10 3 [2,] 2 15 fisher.test(y) Fisher's Exact Test for Count Data data: y p-value = 0.0005367 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 2.753438 300.682787 sample estimates: odds ratio 21.30533 The crude odds

Hi there, I need to do the same thing as cumsum but with the variance and skewness. I have tried to do a loop for like this: var.value <- vector(mode = "numeric", length = length(daily)) for (i in (1:length(daily))) { var.value[i] <- var(daily[1:i]) } But because my dataset is so huge, I run out of memory..... Any ideas?!?! Much appreciate

...TIA, Chris <H>Long question for those with time and wanting more detail</H>: I'm not a statistician though I've spent a lot of the last thirty years of my life working with statisticians and doing my own stats part time. I've worked from roll my own in dBaseII, through BMDP, loads of SPSS, SAS & SAS/IML, S+ and now, over the last two to three years, from S+ to R and loving it. However, at times I really hit the limits of my amateur skills. .... and this is one of them! I've put together some scripts to deal with Jacobson et al.'s ideas of "relia...

SAS has DUD (Does not Use Derivatives)/Secant Method for nonlinear regression, does R offer this option for nonlinear regression? I have read the helpfile for nls() and could not find such option, any suggestion? Thanks, Derek [[alternative HTML version deleted]]

Dear All, The institute I work for is organizing an internal workshop for High Performance Computing (HPC). I am planning to attend it and talk a bit about fluid dynamics, but there is also quite a lot of interest devoted to data post-processing and management of huge data sets. A lot of people are interested in image processing/pattern recognition and statistic applied to geography/ecology, but I