similar to: unhappy with aov performance

Subj. says it. The patch contains several lines complaining about "No newline". After removing these the patch was applicable. -- Dipl.-Math. Wilhelm Bernhard Kloke Institut fuer Arbeitsphysiologie an der Universitaet Dortmund Ardeystrasse 67, D-44139 Dortmund, Tel. 0231-1084-257 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

[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

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='Augenbewegungen'. variables=4.

I tried to reproduce a result from a former colleague which he got with S-plus bootstrap method. I don't have S-plus at hand. In R, there are 2 packages related to bootstrap method, bootstrap and boot. The former has a function called 'bootstrap' but this does not seem to conform either to the function used in S-plus nor to that described in MASS, 3d ed., p.144. The latter seems to be

I just detected, that det() is not working on complex matrices any more, due to the fix to the bug reports noted above. I am not happy with this, as determinants are perfectly usable on complex matrices. AFAIUI the bugs resulted from less than optimal behaviour of qr() in certain cases. IMHO this is due to the unhappy decision to use a default for parameter tol to decide whether the the

I have a problem moving from multistratum aov analysis to lmer. My dataset has observations of ampl at 4 levels of gapf and 2 levels of bl on 6 subjects levels VP, with 2 replicates wg each, and is balanced. Here is the summary of this set with aov: >> summary(aov(ampl~gapf*bl+Error(VP/(bl*gapf)),hframe2)) > >Error: VP > Df Sum Sq Mean Sq F value Pr(>F) >Residuals

Given a matrix B, where B=A'A, how can I find A? In other words, if I have a matrix B which I know is another matrix A times its transpose, can I find matrix A? Thanks, Mike

We get caught occasionally, and we all receive spam; sorry for the last one (on R-help). Of course, there will never be a 100% sure prevention... Question is what should happen when mail from R-help to an adressee is refused with the following message : At least I won't be able to easily send e-mail to them {and am not willing to use a non-easy way}.... ----- Transcript of session

Dear Rxperts.. what settings in barplot and histogram do I use, to show bars in an monotonously increasing or decreasing order of the frequency of a categorical variable? an example is provided below.. histogram(~factor(sample(letters,200,rep=T))) I was able to get it in 3 - 4 steps.. l a1 <- sample(letter,200,rep=T) a2 <- table(a1) barplot(a2); barplot(sort(a2)) # all variables are not

Full_Name: W.B.Kloke Version: 1.9.1 OS: FreeBSD-5.2.1 Submission from: (NULL) (195.253.16.182) Subject line says it. I had problems installing lme4. 1. The dependency on package Matrix was not resolved (I am not sure that this is really a bug; but it is annoying, anyway). 2. Installing Matrix failed with a message saying something like "no rule for %_D.o" after compiling a lot of

On 5/8/2018 12:26 PM, varin sacha via R-help wrote: > > Dear R-experts, > > Here below the reproducible example. I am trying to get the average of the 100 results coming from the "lst" function. I have tried lst$mean and mean(lst). It does not work. > Any help would be highly appreciated > > #################### > > ?## R script for getting MedAe and

Full_Name: W.B.Kloke Version: 1.8.0 OS: FreeBSD-4.7 Submission from: (NULL) (195.253.22.63) I found a parameter evaluation inconsistency in NLME package. I tried to use gls() inside a function, and I wanted use gls() for different subsets of a data frame: prgls <- function(name){ gls( log10(Y)~(cond-1)+(cond-1):t ,pr,subset=subject==name)} Applying this function with a string as parameter

Dear R-experts, Here below the reproducible example. I am trying to get the average of the 100 results coming from the "lst" function. I have tried lst$mean and mean(lst). It does not work. Any help would be highly appreciated. #################### ?## R script for getting MedAe and MedAeSQ from HBR model on Testing data install.packages("robustbase") install.packages(

Dear R-experts, I am trying to do cross-validation for different models using the cvTools package. I can't get the CV for the "FastTau" and "hbrfit". I guess I have to write my own functions at least for hbrfit. What is going wrong with FastTau ? Here below the reproducible example. It is a simple toy example (not my real dataset) with many warnings, what is important to

Dear R-experts, Doing cross-validation for 2 robust regressions (HBR and fast Tau). I can't get the 2 errors rates (RMSE and MAPE). The problem is to predict the response on the testing data. I get 2 error messages. Here below the reproducible (fictional example) R code. #install.packages("MLmetrics") # install.packages( "robustbase" ) # install.packages(

Dear R-experts, I guess I have a problem with my fast function (fast tau estimator) here below. Indeed, zero errors look highly suspicious. I guess there is a bug in my R code. How could I correct my R code ? # install.packages( "robustbase" ) # install.packages( "MASS" ) # install.packages( "quantreg" ) # install.packages( "RobPer" ) #

Chapman & Hall/CRC: The R Series We are delighted to announce that our new series of books on R is up and running, with two books already published and another nine forthcoming (including three set to publish in 2012). We are keen to receive proposals for books covering all aspects of the development and application of R software. If you have an idea for a book, please contact one of the

now with the proper flags passed and the jobs.c fix dash allmost compiles in DEBUG mode, the last two failures happen in show.c, which assumes some buffering functions that klibc doesn't provide: usr/dash/show.o: In function `opentrace': show.c:(.text+0x36): undefined reference to `freopen' show.c:(.text+0x86): undefined reference to `setlinebuf' afais this code dates from the

Dear folks (or Dear Professor Bates), I'm quite confused as to the current status of some of the available functions applicable to lmer objects. Following the examples in Baayen, Davidson, Bates (2006), my plan is to run mcmcsamp on a random effect model created by lmer in package lme4, then use the (perhaps outdated) pvals to estimate p-value. But then I couldn't find pvals anywhere.

Hi, I am on a developpers list for mandrake: cooker@ and the rsync-servers broke which resulted in a lot of very unhappy people cause their rsync directories got empied. My idea is the following. rsync is a mirroring tool. The usual thing that happens in a mirroring tool is that packages do not disappear. You server that archive: like the linux-kernel server. And you got servers that update