similar to: tick marks on both sides of axis

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 ------------------------

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

hello, i did several binomial tests and noticed for one sparse dataset that binom.test(1,1,0.5) gives a p-value of 1 for the null, what i can't quite grasp. that would say that the a prob of 1/2 has p-value of 0 ?? - i must be wrong but can't figure out the right interpretation.. best, kay ----- ------------------------ Kay Cichini Postgraduate student Institute of Botany Univ. of

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

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 R-help, I'm trying to make a fairly simple plot axis that goes something like this: plot(-10:10,-10:10, yaxt='n') axis(side=2, las=1, hadj=1, tck=-.01, cex.axis=.6) ...but as you can see, the labels are not close enough to the y-axis (where I want them... to save space for publication). Can anybody help me figure out how to move these labels over the the right a bit?

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 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

Hello, A made a xyplot using the lattice library in R (latest version). The publisher of our paper has requested: 1. all tick marks should point inwards instead of outwards. 2. All lines should be thicker (lines, axes, boxes, etc. Everything). Lines is easy...I used: lwd=1.5 but what about the lines of the axes, and the lines that build up the plot itself?....? Any

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 all, I have created a barplot that shows change in hardwood/softwood density from 1965 to 2005 in 5 year periods (1965,1970, etc). I would like to have an X-axis where the labels for the years line up after every two bars in the plot (there is one bar for hardwood, and another for softwood). Below is my script: density<-read.table("F:\\Megan\\Vtest.csv", header=TRUE,

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, I try to increase the distance between tick labels and ticks in a lattice wireframe plot. Here's a minimal example: ## Minimal example x <- y <- z <- c(1,2,3) df <- data.frame(x, y, z) wireframe(z ~ x*y, df, scales = list(arrows = FALSE, col = "black", font = 1, tck=0.6)) I tried the axis.components option

Dear R group I am making multiple xyplot, and would like to have tick marks on bottom and left in EACH panel, but only tick labels at the bottom and left of the whole graph. I have browsed the internet, as well as the help page without success. If anyone could help me find the path to the solution I would appreciate. Here is an example, this is the best I could do: lset(col.whitebg()) x.data