search for: buser

Hi I just started using RGui.exe under widnows. I have a text file containing date arranged in columns and rows, each column has the same format, each row with different formats. 3 of the columns are something like this 1/12/2006 3:59:45 PM I need to calculate the different in seconds between 2 selected periods using their row’s index My solution: Read the file in a data frame and

...list of factors, newly returns a single factor instead of an integer. I am happy with this change, but: Is it desirable to change apply so that it does not return a character matrix in the example above or include a warning for such a case? Thank you very much for an answer. Regards, Christoph Buser -------------------------------------------------------------- Christoph Buser <buser at stat.math.ethz.ch> Seminar fuer Statistik, LEO C13 ETH Zurich 8092 Zurich SWITZERLAND phone: x-41-44-632-4673 fax: 632-1228 http://stat.ethz.ch/~buser/

...he output is affected by the argument "adjustSigma", but not the estimate for residual standard error (see the artificial example below). May someone explain to me if there is an error in my understanding of the help page and the R code? Thank you very much. Best regards, Christoph Buser -------------------------------------------------------------- Christoph Buser <buser at stat.math.ethz.ch> Seminar fuer Statistik, LEO C13 ETH Zurich 8092 Zurich SWITZERLAND phone: x-41-44-632-4673 fax: 632-1228 http://stat.ethz.ch/~buser/ -------------------------------------------------...

Hello :) I have two data frames, one has properties taken on a piece by piece basis and the other has performance on a lot by lot basis. I wish to combine these two data frames but the problem is that each lot has multiple pieces and hence i need to take a mean of the properties of multiple pieces and match it to the row having data about the lot. I was wondering if there is a simple commmand,

Hello everbody, i'm currently using the anova()-test for a small data.frame of 40 rows and 2 columns. It works well, but is there any preconditions for a valid variance analysis, that i should consider? Thank you for your answer, Daniel

...ing with this argument main and creates this warning. You could replace ,... in the axis function by inserting all reasonable arguments of axis in the functions header (interaction.plot) and give over those arguments to axis. Then you shouldn't get this warning anymore. Best regards Christoph Buser -- Christoph Buser <buser@stat.math.ethz.ch> Seminar fuer Statistik, LEO D6 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-1-632-3505 fax: 632-1228 http://stat.ethz.ch/~buser/

...(similar to col = NA): > plot(1:9, col = factor(rep(1:3,3), labels = c("red", "blue", "black"))) My question: Shouldn't one get at least a warning (or an error) if one tries to use a factor as col argument? Thanks for an answer. Regards, Christoph Buser -------------------------------------------------------------- Christoph Buser <buser at stat.math.ethz.ch> Seminar fuer Statistik, LEO C13 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-44-632-4673 fax: 632-1228 http://stat.ethz.ch/~buser/

hi list, i have to ask you again, having tried and searched for several days... i want to do a TukeyHSD after an Anova, and want to get the adjusted p-values after the Tukey Correction. i found the p.adjust function, but it can only correct for "holm", "hochberg", bonferroni", but not "Tukey". Is it not possbile to get adjusted p-values after

...d like to use arms. But there seem to be a strong influence of the initial point y.start. To show the effect I constructed a demonstration example. It is reproducible without further information. Please note that my target density is not logconcave. Thanks for all comments or ideas. Christoph Buser ## R Code: library(HI) ## parameter for the distribution para <- 0.1 ## logdensity logDichteGam <- function(x, u = para, v = para) { -(u*x + v*1/x) - log(x) } ## density except for the constant propDichteGam <- function(x, u = para, v = para) { exp(-(u*x + v*1/x) - log(x)) } ## calc...

...t different DF for the denominator. Do I have to use lme in another way in the package lme4? I use R 2.0.1 under linux and Package: nlme Version: 3.1-53 Date: 2004-11-03 Package: lme4 Version: 0.6-11 Date: 2004-12-16 Thanks for help. Regards, Christoph Buser -- Christoph Buser <buser at stat.math.ethz.ch> Seminar fuer Statistik, LEO C11 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-1-632-5414 fax: 632-1228 http://stat.ethz.ch/~buser/

I have some doubts about the validity of my procedure to estimeate linear contrasts ina a factorial design. For sake of semplicity, let's imagine a one way ANOVA with three levels. I am interested to test the significance of the difference between the first and third level (called here contrast C1) and between the first and the seconda level (called here contrast C2). I used the following

....9 256.8 ##- Litter 3 63.6 3328.5 252.0 0.3508 0.78870 ##- Mother 3 775.1 4040.0 263.8 4.2732 0.00886 ** ##- --- ##- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 ##- Litter stays not signifikant. Best regards Christoph Buser -- Christoph Buser <buser at stat.math.ethz.ch> Seminar fuer Statistik, LEO D6 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-1-632-3505 fax: 632-1228 http://stat.ethz.ch/~buser/

Hello Could you please suggest a way to find out the r square values for each independent variable while using "lm" for developing a multiple linear regression model. Thank you avneet "I have no data yet. It is a capital mistake to theorize before one has data. Insensibly one begins to twist facts to suit theories instead of theories to suit facts." ~ Sir Arthur Conan Doyle

...status Under development (unstable) major 2 minor 1.0 year 2005 month 03 day 22 language R Regards, Christoph Buser -------------------------------------------------------------- Christoph Buser <buser at stat.math.ethz.ch> Seminar fuer Statistik, LEO C11 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-1-632-5414 fax: 632-1228 http://stat.ethz.ch/~buser/

Hi, I have a model with differents observations Xi. Each observation belongs to a group, either A or B. I would like to minimize a fonction like : sum( Xi - Z)^2 + sum (Xi - aZ -b)^2 A B The first sum contains all observations from group A and the second all observations from group B. I want to find the Z-value wich minimize this function. a and b are predefined parameters.

hi all, I wish to draw on the same graphic device 3 functions. But i don't want them to be on different graph, i want to compare them on the same I don't need mfrow or mfcol, I need something else... 1 graph on 1 device inside this graph 3 ploted function. I saw something unsing data.frame, but i think it's overkill, and something less complicated must exist, if not why? why not

Hi, I use anova() to find interaction effect and then I need to find data with main effect from those data with interaction effect. How to do that? I used : anova(lm(data~factor1*factor2)), then select data with interaction effect. Then I need to select those data also with main effect of factor1 or factor2, from previous selected data. How to do that? Many thanks for your time on my

Hi, I'm studying math, and i have to make an analysys using Randomised Block Design. I have two factors, i know how to do this in Statistica, but how to do this in R, i read some manuals but the only thing that i have found was 2 factor ANOVA. Please could someone help me, or give some usefull links ?? Krzytsztof Suwada

Hi All. I need a bit of help extracting the residual error variance from the VarCorr structure from lmer. #Here's a 2-way random effects model lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat) #Get the structure vc.fit <- VarCorr(lmer.1) #results in..... $person 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 0.7755392 $rater 1 x 1 Matrix

useRs, A no doubt simple question, but I am baffled. Indeed, I think I once knew the answer, but can't recover it. The default contrasts for aov (and lm, and...) are contr.treatment and contr.poly for unordered and ordered factors, respectively. But, how does one invoke the latter? That is, in a data.frame, how does one indicate that a factor is an *ordered* factor such that