similar to: Testing continuous zero-inflated response

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,

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

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

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

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 *

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

Hi R users, Could anyone let me know how to export a paired t-test output table (see below) to an excel file? Jason, with(score2,pairwise.t.test(values,ind, + p.adjust.method="holm", paired=T)) Pairwise comparisons using paired t tests data: values and ind test1 test2 test3 test4 test5 test6 test1 1.0000 - - -

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 =

Hola Javier, Si, cuando hablo de valor observado me refiero al valor real en campo y el predicho al que estiman los modelos. Disculpa, que no lo detallase así desde el principio. En mi caso trabajo con dos diferentes: Zero inflated y Binomial Negativo y me gustaría comprobar que diferencia (distancia) existe entre cada uno de ellos y la realidad. Estoy trabajando con los siguientes paquetes:

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

Buenas tardes: Soy Pedro José Martínez, del Departamento de Economía Financiera y Contabilidad de la Universidad de Murcia y tengo una duda sobre R. Mi investigación se basa en calcular la eficiencia de los servicios municipales con DEA (packages deaR), y posteriormente identificar los determinantes de dicho nivel con una regresión truncada (packages truncreg). Para calcular la eficiencia

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

Post-hoc test for a zero inflated continuous data set with a tweedie distribution? I have a zero inflated continuous data set of aphid feeding duration on 10+ species of plant. I have fitted a glm model with a tweedie distribution and used anova() function to show that there is significance between the plant species. However, I would now like to perform of post-hoc test, ideally a Tukey-Kramer