search for: drechsler

Hi, I was wondering if there is a ready made function or parameter for indicating the sample size in boxplots? Here's what I came up with so far: library(ISwR) data(energy) attach(energy) boxplot(expend~stature) sample.size <- tapply(expend, stature, length) sample.size <- paste("N=", sample.size, sep="") mtext(sample.size, at=1:length(unique(stature)), line=2,

...LC_IDENTIFICATION=C | | attached base packages: | [1] "stats" "graphics" "grDevices" "utils" "datasets" "methods" | [7] "base" | | other attached packages: | lattice | "0.14-17" `---- -- Patrick Drechsler Department of Zoology University of Cambridge Downing Street Cambridge CB2 3EJ, UK

...lot(foo) ## 3. plot "ancova" plot using HH: ## normal usage: ancova(y ~ x1.cont * x2.fac) ## 4. Trying to combine "ancova" with trellis: ## does not work: ancova(y ~ x1.cont * x2.fac * x3.fac) --8<---------------cut here---------------end--------------->8--- -- Patrick Drechsler Department of Zoology University of Cambridge Downing Street Cambridge CB2 3EJ, UK

Hi, is there a grayscale setting for lattice plots? I like the default color settings. I also like the settings that are available for setting black and white with something like this: --8<---------------cut here---------------start------------->8--- ltheme <- canonical.theme(color = FALSE) ## in-built B&W theme ltheme$strip.background$col <- "transparent" ##

Hi, I was wondering if it's possible to have fonts in plots to be autoscaled to the same font size used by LaTeX in a surrounding Sweave document. Here's a short example in which the fonts of the first plot are barely readable: --8<------------------------schnipp------------------------->8--- \documentclass{article} \usepackage{graphicx} \newcommand{\mytext}{Some normal sized

Hi, I am performing a two-way ANOVA (2 factors with 4 and 5 levels, respectively). If I'm interpreting the output of summary correctly, then the interaction between both factors is significant: ,---- | ## Two-way ANOVA with possible interaction: | > model1 <- aov(log(y) ~ xForce*xVel, data=mydataset) | | > summary(model1) | Df Sum Sq Mean Sq F value Pr(>F) |

Hi, sorry if this is OT here. I'm trying to compile StatDataML[1] to import Matlab data into R. The R part works fine but compiling the Matlab part exits with this error message (following the INSTALL instructions: autoconf -> ./configure -> make): ,----[ patrick at trurl:~/src/statdatml/StatDataML/MatOct> make ] | make[1]: Entering directory

Hi, I'm having trouble getting the function `npmc' working. I have a 2 column 30 row table called `mydata': > npmc(mydata, df=2, alpha=0.05) Error in sort(unique.default(x), na.last = TRUE) : `x' must be atomic > is.atomic(mydata) [1] TRUE Is there anything I have overlooked? Am I invoking the function correctly? I wasn't able to find any hints in the archives

Hi, can anybody tell me how to contact the maintainer of the "concord" package? The address given in the help file is not valid anymore. ,---- | > help(package=concord) | | Information f?r Paket 'concord' | | Description: | | Package: concord | Version: 1.4-2 | Date: 2005-05-10 | Title: Concordance and reliability | Author: Jim Lemon

Hi, what is the recommended way of adding a strip to a lattice plot? In the example below I would like to add the value of mean(y) to a new strip.: --8<---------------cut here---------------start------------->8--- library(lattice) ## Small sample data set: p0 <- xyplot(uptake ~ Type | Treatment, data = CO2) p1 <- update(p0, panel = function(x, y, ...) {

I've searched for 45 minutes, apparently in all the wrong places for a solution to a configuration issue I'm having. When I use Xemacs with ESS running in R-mode, and I type a "-" character, it autocompletes it to "<- ". How do I disable this annoying "feature"? Thanks much. [[alternative HTML version deleted]]

Hi, is there an R function to convert unixtime to one of the R time formats (using chron or POSIXct)? Example data: unixtime year month day hour 1183377301 2007 7 2 13 I would like to only use the first column. If such a function does not exist: What would be a good input format for R (I could parse the input using AWK or other GNU/Linux tools). TIA,

Hi, I just tried to compile the new release today and encountered some problems: --8<------------------------schnipp------------------------->8--- patrick@trurl:~/src/r-project/R-2.0.0> ./configure > pd.configure.output patrick@trurl:~/src/r-project/R-2.0.0> make > pd.make.output /usr/lib/gcc-lib/i486-suse-linux/3.3/../../../../i486-suse-linux/bin/ld: tcltk.so: undefined

Hi, I'm having trouble sweaving the following example: --8<------------------------schnipp------------------------->8--- \documentclass{article} \begin{document} \SweaveOpts{echo=FALSE} @ <<>>= rm(list=c(ls())) library(Hmisc) library(ISwR) data(energy) energy$log <- log(energy$expend) attach(energy) @ %def @ <<results=tex>>= mysum <- summary(stature ~

Hi, for testing if the slope of experimental data differs from a given slope I'm using the function "test_regression_against_slope" (see below). I am now confronted with the problem that I have data which requires a modelII regression (also called reduced major axes regression (RMA) or geometric mean regression). For this I use the function "modelII" (see below). What

Hi all, I've just started using the Friedrich Leisch's Sweave package to generate LaTeX reports with results of my analyses embedded as required. I've encountered a bit of a problem though in the processing of the resulting *.tex file that is as far as I can tell, down to the fact that my variable names have underscores ('_') in their names. The relevant section of code

--=-=-= Sorry, forgot to attach the file... --=-=-= Content-Type: application/zip Content-Disposition: attachment; filename=pdcompilelog.zip Content-Transfer-Encoding: base64 UEsDBBQAAAAIAImkRTGksC5Kaw8AACBOAAATABUAcGQuY29uZmlndXJlLm91dHB1dFVUCQADoeli QY3pYkFVeAQA9AFkAOVcbY/bNhL+nl8hFAc0vcR2Nt1kgwJ3QLLrtm73JdjNHRb3xaApymYskapE

...ror of mixed models (Bernardo Rangel Tura) 2. Where to ask general stats questions? (francogrex) 3. Re: structure vs. matrix (Mark Difford) 4. Problems with garch() function tseries package R 2.6.0 (Jos? Augusto Morais de Andrade J?nior - JAMAJ) 5. Re: Where to ask general stats questions? (Patrick Drechsler) 6. help on error message!!! (david csongor) 7. Re: structure vs. matrix (Duncan Murdoch) 8. Re: structure vs. matrix (Duncan Murdoch) 9. Re: help on error message!!! (Duncan Murdoch) 10. Re: help on error message!!! (Bernardo Rangel Tura) 11. Re: Problems with garch() function tseries package R 2....

Hi, I would like to change font sizes (more generally the cex) in plots produced by the function "interaction2wt" of the "HH" package. Here is an example (from the interaction2wt example) of what I have tried: --8<---------------cut here---------------start------------->8--- rm(list = ls(all = TRUE)) rm(list = c(ls())) library("HH") ## From the example: