similar to: general question on binomial test / sign test

hello, i want to stack two lattice plots beneath each other using one x-axis and sharing the same text-panels, like: ##################################################################### library(lattice) y1 <- rnorm(100,100,10) y2 <- rnorm(100,10,1) facs<-expand.grid(Sites=rep(c("Site I","Site II"),25),Treatment=c("A","B")) pl1<-dotplot(y1 ~

hello, i need to add legend text: "4th-root transformation", with the "th" superscripted - tried much - but nothing worked.. thanks for any hints, kay ----- ------------------------ Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck ------------------------ -- View this message in context:

I want to replace \\ with \ in: str <- "C:\\DOKUME~1\\u0327336\\LOKALE~1\\Temp\\RtmpQ5NJ8X\\TIRIS_PICS\\1_Img.jpg" and tried: gsub("\\\\", "\\", str) but this removes the \\ without replacing them by \ Any help much appreciated, Kay ----- ------------------------ Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck ------------------------

Hallo List, I'm trying to implemement a restricted permutation scheme in permutest(). More precisely I have dependence in my data that should be allowed for in the permutation - I simulated the problem in the example of the vegan documentation p.24: library(vegan) data(varespec) ## Bray-Curtis distances between samples dis <- vegdist(varespec) ## First 16 sites grazed, remaining 8 sites

hello, all my attempts to get a plot with tick marks on both sides of axis failed - can someone please help me with this? thanks, kay ----- ------------------------ Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck ------------------------ -- View this message in context: http://r.789695.n4.nabble.com/tick-marks-on-both-sides-of-axis-tp2336668p2336668.html Sent from the

hello, i'd like to evaluate the 2 factors (f1 with 4, f2 with 2 levels) and its interactions in my model for writing up a summary table with likelihood-ratio tests for publication purpose. now i'm unsure about which terms of my model may be dropped from the full model and which not. mfull<-lme(resp~f1*f2,random=...) m1<-lme(resp~f1+f2,..) m2<-lme(resp~f1,..)

hello, i'd like to extract values from stressplot, but as this should happen within a long run of code, which also does a multiple plot, i need to prevent stressplot to print to the device. i tried ..,plot=F), but this is useless. any advise on this? thanks, kay ----- ------------------------ Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck

hi, i'd like to install the package "RGoogleDocs ". i downloaded to path "E:/R/R-2.13.0/library/RCurl_0.91-0.tar.gz" i run R from an usb-stick and can't get the install.packages() prompt to run correctly - can anyone help with this? thanks, kay > sessionInfo() R version 2.13.0 (2011-04-13) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1]

hello list, i'd very much appreciate help with setting up the contrast for a 2-factorial crossed design. here is a toy example: library(multcomp) dat<-data.frame(fac1=gl(4,8,labels=LETTERS[1:4]), fac2=rep(c("I","II"),16),y=rnorm(32,1,1)) mod<-lm(y~fac1*fac2,data=dat) ## the contrasts i'm interressted in: c1<-rbind("fac2-effect in

hello everyone, i sampled 100 stands at 20 restoration sites and presence of 3 different invasive plant species. i came across logistic regression trees and wonder if this is suited for my purpose - predicting presence of these problematic invasive plant species (one by one) by a set of recorded ecological / geographical parameters. i'd be glad if someone would comment on applying this

Hi everyone! I'm matching two samples to create one sample that have pairs of observations equal for the k1 variable. Merge() doesn't work because I dont't want to recycle the values. x <- data.frame(k1=c(1,1,2,3,3,5), k2=c(20,21,22,23,24,25)) x y <- data.frame(k1=c(1,1,2,2,3,4,5,5), k2=c(10,11,12,13,14,15,16,17)) y merge(x,y,by="k1") k1 k2.x k2.y 1 1 20

HI everyone, I try to get some bootstrap CIs for coefficients obtained by quantile regression. I have influencial values and thus switched to quantreg.. The data is clustered and within clusters the variance of my DV = 0.. Is this sensible for the below data? And what about the warnings? Thanks in advance for any guidance, Kay > dput(d) structure(list(Porenfläche = c(4990L, 7002L, 7558L,

hello, i found it most convenient to use package contrast for planned comparisons on mixed models. for instance i have a model with 2 fixed factors, one with 4 levels (stage) and one with 2 levels (gap) and a nested random factor (site) and i tested gap within level A of factor stage: library(contrast) library(nlme) m1<-lme(rich ~ stage*gap, random=~1|site,data=richness) contrast(m1, a =

Hello, I'm searching for a test that applies to a dataset (N=36) with a continuous zero-inflated dependent variable and only one nominal grouping variable with 2 levels (balanced). In fact there are 4 response variables of this kind which I plan to test seperately - the amount of zeroes ranges from 75 to 97%.. I searched the web and found several modelling approaches but have the feeling

Hi all, how would you test if a sample mean of integers with range -inf;inf is different from zero: # my sample of integers: c <- c(-3, -1, 0, 1, 0, 3, 4, 10, 12) # is mean of c <> 0?: mean(c) Thanks, Kay [[alternative HTML version deleted]]

Hello list, I' d very much appreciate some help with a two sample repeated measures ANOVA. I did the analysis yielding sign. main effects (between subj.=site, within subj.=cover) and a sign. interaction: Univariate Type II Repeated-Measures ANOVA Assuming Sphericity SS num Df Error SS den Df F Pr(>F) site 18.7620 1 18.831 10 9.9631 0.010220 *

hello, i'd be happy if someone could provide help with the following problem: i have a dist.matrix that comes from vegdist() function of the vegan package. the used method = "horn" is not accepted as argument in hvcluster(...,dist.method="..."). is there a way to incorporate the method "horn" in hvcluster()? thanks in advance! yours, kay -- View this

hello, i can't figure out how to set position of panels of my stripplot - i`d like the panels of one level of the factor stage (nr. of panels within each stage, A: 12, B: 12, C: 12, D: 4, each panel representing a site) to be in one column, with A to D from left to right and with descending site.nr at each row. like: A1 B1 C1 D1 A2 B2 .. .. A3 .. .. .. how is this achieved? any help

I'm trying to compute sample size requirements for a binomial exact test. we want to show that the proportion is at least 90% assuming that it is 95%, with 80% power so any asymptotic approximations are out of the questions. I was planning on using binom.test to perform the simple test against a prespecified value, but cannot find any functions for computing sample size. do any exist?

hello, i can't figure out how to change the vertical position of my x axis labels.. boxplot(c(1:12)~c(rep("1",6),rep("2",6)),at=c(1,2), col=c(0,"grey"),las=1,xaxt="n") ### i put paragraphs in the x-labels because of limited horizontal space axis(1,at=c(1,2),adj=1,labels=c("Salix Scrub","Tall Forb")) ...the labels are one line