similar to: lattice help required

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

Dear all, let consider the following function: Fun1 <- function() { library(lattice) plot1 <- 1:10~1:10 pl1 <- xyplot(plot1) return(pl1$call$x) } In R 2.5.0 (or older version) we have > Fun1() plot1 but starting from R 2.5.1 until the latest R 2.8.0 we obtain instead > Fun1() NULL because pl1$call seems to be equal to xyplot() without arguments. Something like

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

Posted on behalf of chenwj <stop4optimal at hotmail.com> who is having trouble posting. hi all, I'm completely new in R, and here are my questions (under windows): 1. i want to run my four R files (pl1.r -- pl4.r) in batch mode. i'd like to write a batch file using Rterm.exe, since i dont get perl installed on my computer. My experimental batch file (*.bat) didnt work, which

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

Greetings R Help Group, How does one effect a multiphasic logistic growth model with 4 phases (e.g. Koops 1986; Weigel, Craig, Bidwell and Bates 1992; Grossman and Koops 2003) with R. Before writing to the group, the R help archives were searched, the web was searched with Google, Venables and Ripley 2002 was consulted, Pinheiro and Bates 2000 was consulted, Bates and Watts 2007 was bought and

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, I'm trying to generate a log-log plot with ggplot. Within this plot, I would like to draw a straight line with the help of ggabline which does not work. The code is pl1 <- qplot(correlFunc, correlFunc.ref, data=all, log=c("xy")) pl1 <- ggabline(pl1, slope=1, intercept=0) print(pl1) ## no line in plot :( I can produce plots as I want, if I transform all values before

Hi R Users, When I use package lme4 for mixed model analysis, I can't distinguish the significant and insignificant variables from all random independent variables. Here is my data and result: Data: Rice<-data.frame(Yield=c(8,7,4,9,7,6,9,8,8,8,7,5,9,9,5,7,7,8,8,8,4,8,6,4,8,8,9), Variety=rep(rep(c("A1","A2","A3"),each=3),3),

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