[R] How to set up a Negative Exponential regression mixed effect model 'nlme' ?

2016-11-17 Thread Ronny Steen
Hi, I need help with fitting a non-linear mixed effects model (nested random effect). I look at the relationship between 'wingbeat frequency' (beats per sec) and 'wing length' in hummingbirds, and I will find a model with the best fit. I have posted the question at R-SIG-mixed-models list, but no

[R] Converting time zones in R using metadata file information of video files, help needed.

2015-11-22 Thread Ronny Steen
Hi, I have video files (FAT) that are taken in a different timezone than my current location. The modification date/time in the metafila data of the video file shows the time video was taken, although in the current timezone of my computer, if I understand right. I wish to convert the date/time t

[R] degrees of freedom in custom contrasts ananlysis

2013-08-12 Thread ronny
Hi, I would like to understand why the residual standard error, and the degrees of freedom are changing when I define custom contrasts, which are not orthogonal. For example: y <- rnorm(40) x <- factor(rep(1:10,4)) summary(lm(y~x)) #standard model: Residual standard error: 1.103 on 30 degrees of

[R] Lattice levelplot- remove unused levels per panel

2013-01-20 Thread ronny
Hi, I am using levelplot, and would like remove from each panel (condition) its unused x levels. e.g. Remove from panel vs=1 the cyl level=8. data(mtcars) levelplot(mpg~factor(cyl)*factor(gear)|factor(vs)) Thanks for your help, Ronny -- View this message in context: http://r.789695.n4

[R] problematic p values for Dixon Q test

2012-09-13 Thread ronny
I am using R dixon.test in order to perform Q test for detection of outliers. The p-values for high Q values are consistent with the table in Rorabacher, D.B. (1991) http://www.flworkshop.com/sscs/dixon_test.pdf, but it seems that small Q-values get problematic p-values. For example: dixon.test(c

[R] Bootstrap values for hierarchical tree based on distaance matrix

2011-11-15 Thread ronny
I would like to get an hierarchical clustering tree with bootstrap values indicated on the nodes, as in pvclust. The problem is that I have only distance matrix instead of the raw data, required for pvclust. Is there a way to get it? fit1 <- hclust(dist) # an object of class '"dist" plot(fit1) # d

Re: [R] Centering data frame by factor

2011-07-19 Thread ronny
Perfect! Made my day! -- View this message in context: http://r.789695.n4.nabble.com/Centering-data-frame-by-factor-tp3677609p3677665.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.c

[R] Centering data frame by factor

2011-07-19 Thread ronny
Hi, I would like to center P1 and P2 of the following data frame by the factor "Experiment", i.e. substruct from each value the average of its experiment, and keep the original data structure, i.e. the experiment and the group of each value. RAW= data.frame("Experiment"=c(2,2,2,1,1,1),"Group"=c(

Re: [R] How to convert S4 class slots into data.frame or how to assign variables of type 'Date'

2008-11-21 Thread Ronny Wölbing
ity. Do you have any suggestions? Is there maybe even a simple way to convert S4 slots into data.frames? Thanks in advance Ronny __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] How to convert S4 class slots into data.frame or how to assign variables of type 'Date'

2008-11-20 Thread Ronny Wölbing
the type of DateOfTransaction as 'Date'. But I don't see a way to use this assigning scheme in my method without using the actual slotnames and giving up a lot of flexibility. Do you have any suggestions? Is there maybe even a simple way to convert S4 slots into data.frames?

[R] Guidance for reporting results from lme test?

2008-01-29 Thread Ronny Steen
regards Ronny _ [[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