search for: eawag

...temp0$z, xlab="x", ylab="y", main="z") temp <- kde2d.weighted(x=x, y=y, w=z, n = 100, lims=c(range(x),range(y))) contour(x=temp$x, y=temp$y, z=temp$z, xlab="x", ylab="y", main="z", col="red", add=TRUE) ??? Christoph Ort Eawag - aquatic research Environmental Engineering ?berlandstrasse 133 CH-8600 D?bendorf phone: +41-44-823-5041 fax: +41-44-823-5389 cell: +41-79-218-9242 mailto:christoph.ort at eawag.ch http://www.eawag.ch/~ortchris/ http://www.eawag.ch ??? Christoph Ort Eawag - aquatic research...

...? (3) By the way, how to deal with low frequency cells? r <- c(10, 100, 500, 1000, 2000, 5000) v <- c(35, 40, 45, 45, 40, 35) sapply(list(r), function (x) { chisq.test(v, p=c(rep.int(40, 6)), rescale.p=T, simulate.p.value=T, B=x)$p.value }) Thank you, S?ren -- S?ren Vogel, PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...nted: outputs the string properly. Wanted: assigns the list properly. How can I get rid of the *unwanted* part? That is, how do I define what the functions prints and -- on the other hand -- what it returns without printing? Thanks S?ren -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...applied over each var according to the grouping? That is: Desired output table: A B C r count count count s ... ... ... or likewise transposed. I'd like to use the table for textual display and/or barplot creation. Thx, S?ren -- S?ren Vogel, PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...t; ### example:end What I want in the end is the following: m$locA: 0, 0, 0, 0, 0, 1, 1, 1, 0, 1 m$locB: 0, 0, 0, 0, 1, 0, 0, 0, 0, 0 m$locC: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 m$locD: 1, 1, 1, 1, 0, 0, 0, 0, 1, 0 Instead of 0, NA's would also be fine. Thanks, S?ren -- S?ren Vogel, PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...ows an s-shaped distribution when using QQ-plots (some zero values and some very high values). hence, which transformation should i apply that Y follows a normal distribution? any r-function/package available to do this? thanks for any hint, regards, lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220 38 25 Fax : +41 (0) 44 823 53 15 Email: lukas.indermaur at eawag.ch www.lukasindermaur.ch <http://www.lukasindermaur.ch/&...

...a batch job: # e.g. model 1 preference<- formula(Y~Lwd + N + Sex + YY) # e.g. model 2 preference_heterogeneity<- formula(Y~Ri + Lwd + N + Sex + YY) etc. etc. I appreciate any hint Cheers Lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220 38 25 Fax : +41 (0) 44 823 53 15 Email: lukas.indermaur at eawag.ch www.lukasindermaur.ch

...practical for me since the number of bytes before the specific line varies between the files, therefore I use reanLines() and match the right line instead.) My Systems: OS: Windows Server 2003 & Linux Red Hat (interchangeably) R version: 2.7.2 Thanks for any suggestions! ooo Jafet Andersson Eawag - The Swiss Federal Institute of Aquatic Science and Technology Ueberlandstrasse 133 P.O. Box 611 CH-8600 Duebendorf Switzerland Phone: +41 (0)44 823 5358 Fax: +41 (0)44 823 5028 http://www.eawag.ch/index_EN

Hello Instead of testing against 0 i would like to test regression slopes against -1. Any idea if there's an R script (package?) available. Thanks for any hint. Cheers Lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220 38 25 Fax : +41 (0) 44 823 53 15 Email: lukas.indermaur at eawag.ch www.lukasindermaur.ch

...matrix cell: m <- matrix(1:9, nrow=3, dimnames=list(LETTERS[1:3], letters[1:3])) l <- list(m1=m, m2=m*2, m3=m*3) l[[3]] # works l[[3]][1:2, ] # works l[[1:3]][1, 1] # does not work How can I slice all C-c combinations in the list? S?ren -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM +41 76 233 3637, http://www.eawag.ch, http://sozmod.eawag.ch

...el) } #Writes results to data file x = data.frame( modelnames, n, AICA) write.table(x, paste(t.url, file="Results.txt"), sep="\t", quote=F) ----------------------------------------END R-CODE----------------------------------- ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220 38 25 Fax : +41 (0) 44 823 53 15 Email: lukas.indermaur at eawag.ch www.lukasindermaur.ch

...e variables "mama" and "papa" correctly. But how do I remove all "Mama" list entries in "a" in the same run, that is, shrink the list by what was already matched? Thank you for your help! S?ren Vogel -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...re stored in a dataframe called "mydataframe". How can I implement my script simply into a bootstrap and sum up the frequencies of variables of the 10'000 bootstraped results? I am not so good in R and would appreciate any hint. Regards Lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220 38 25 Fax : +41 (0) 44 823 53 15 Email: lukas.indermaur at eawag.ch www.lukasindermaur.ch

### example:start v <- sample(rnorm(200), 100, replace=T) k <- rep.int(c("locA", "locB", "locC", "locD"), 25) tapply(v, k, summary) ### example:end ... (hopefully) produces 4 summaries of v according to k group membership. How can I transform the output into a nice table with the croups as columns and the interesting statistics as lines? Thx,

...sample(1:3, 14, repl=T)), 'x3'=round(rnorm(20, 4, 2), 0) ) df[17:18, c(2, 3)] <- df[17:18, c(3, 2)] glm(y ~ ., data=df[, -3], family=binomial("logit")) glm(y ~ ., data=df, family=binomial("logit")) Thanks, S?ren -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

Hello The "Craddock-Flood Test" is recommended for large tables with small degrees of freedom and low-frequency cells. Is there an R procedure and/or package which does the test? Thank you for your help! S?ren Vogel -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...pression (http://www.mail-archive.com/r-help at r-project.org/msg55480.html), which in turn is quite outdated. Is the whole idea dead or did I miss something important when googling? *S* (R version 2.11.1 (2010-05-31), x86_64-apple-darwin9.8.0) -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...;=Rcs, 'Rnagel'=Rnagel) class(out) <- c("list", "table") return(out) } y <- sample(c(T, F), 50, repl=T) x <- sample(1:7, 50, repl=T) mod <- glm(y ~ x, family=binomial("logit")) Rcsnagel(mod) -- S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM http://www.eawag.ch, http://sozmod.eawag.ch

...andom = ~ 1 | ANIMAL, correlation=corCAR1(0.2, form = ~1 | ANIMAL ), method='ML', na.action=na.omit)). Is there a possability to easly compute an R-square from the output of the model summary? I would appreciate any hint. Best regards Lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220 38 25 Fax : +41 (0) 44 823 53 15 Email: lukas.indermaur at eawag.ch www.lukasindermaur.ch

Dear All Where can I lookup good methods to compute p from q=bin(m,n)p^n*(1-p)^(m-n) such that q<=alfa, alfa small. Are there such libs, code and source in R? Best Regards -- Jan Burse SIAM, EAWAG Scheuchzerstr. 67 ?berlandstr. 133 8006 Z?rich 8600 D?bendorf tel: +41-1-364 17 66 tel: +41-1-823 55 34 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R...