Re: [R] syntactical meaning of fullstop in R functions

2013-07-12 Thread Kay Cichini
just found it myself: in '.GADM' the leading period designates an internal function - the source can be viewed with: getAnywhere('.GADM') 2013/7/12 Kay Cichini > hello, > > what is the meaning of a fullstop in the below R-Ffunction (like in > '.GADM&

[R] syntactical meaning of fullstop in R functions

2013-07-12 Thread Kay Cichini
} else if (name == "countries") { .countries(download = download, path = path, ...) } else { stop(name, " not recognized as a valid name.") } } thanks in advance, kay -- Kay Cichini, MSc Biol Grubenweg 22, 6071 Aldrans Tel.: 0650 9359101 E-Mai

Re: [R] Testing continuous zero-inflated response

2013-01-28 Thread Kay Cichini
Many thanks - this was very helpful! Regards, Kay Am 28.01.2013 13:19 schrieb "Achim Zeileis" : > On Sun, 27 Jan 2013, Kay Cichini wrote: > > That said, >> >> wilcox_test(x ~ factor(y), distribution = "exact") >>> >> >> or t

Re: [R] Testing continuous zero-inflated response

2013-01-27 Thread Kay Cichini
That said, > wilcox_test(x ~ factor(y), distribution = "exact") or the same with oneway_test, i.e would be ok? 2013/1/27 Achim Zeileis > On Sun, 27 Jan 2013, Kay Cichini wrote: > > Thanks for the reply! >> >> Still, aren't there issues with 2-sample

Re: [R] Testing continuous zero-inflated response

2013-01-27 Thread Kay Cichini
Thanks for the reply! Still, aren't there issues with 2-sample test vs y and excess zeroes (->many ties), like for Mann-Whitney-U tests? Kind regards, Kay 2013/1/26 Achim Zeileis > On Fri, 25 Jan 2013, Kay Cichini wrote: > > Hello, >> >> I'm searching for a

[R] Testing continuous zero-inflated response

2013-01-25 Thread Kay Cichini
anges from 75 to 97%.. I searched the web and found several modelling approaches but have the feeling that they are overly complex for my very simple dataset. Thanks in advance for any help! Kay -- Kay Cichini, MSc Biol Grubenweg 22, 6071 Aldrans Tel.: 0650 9359101 E-Mail: kay.cich...@gmai

Re: [R] bootstrapping quantile regression

2012-10-30 Thread Kay Cichini
sry, I forgot to replace rlm() - but actually I tried both and the question applies to both approaches.. Am 31.10.2012 00:19 schrieb "Kay Cichini" : > > HI everyone, > > I try to get some bootstrap CIs for coefficients obtained by quantile regression. I have influencial va

[R] bootstrapping quantile regression

2012-10-30 Thread Kay Cichini
"47", "48"), class = "data.frame") ## do quantile regression and bootstrap the coefficients, allowing for clustered data ## by putting "Gebiet" as strata argument (?), ## dv variation within clusters/Gebiet = 0! bs <- function(formula, data, indices) { d <- data

Re: [R] How can I export a paired t-test output table to an excel file?

2012-06-03 Thread Kay Cichini
If you wish to use the table in a write-up or something like this, there is a package for exporting to MS Word - see http://cran.r-project.org/web/packages/R2wd/index.html http://cran.r-project.org/web/packages/R2wd/index.html and these examples http://thebiobucket.blogspot.co.at/2012/04/reprod

[R] Test if a sample mean of integers with range -inf; inf is different from zero

2012-05-04 Thread Kay Cichini
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]] _

[R] replace double backslash with singel backslash

2011-11-04 Thread Kay Cichini
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 ----- ----

[R] Package/Function for Blending Images

2011-10-11 Thread Kay Cichini
Hi there, Does someone know a package/function for blending two pictures or to add transparency.. Thanks in advance, KC - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r

Re: [R] interaction contrasts

2010-10-15 Thread Kay Cichini
,0,0,0,0,0,1,0), "fac2-effect, A*D"=c(0,0,0,0,0,0,0,1), "fac2-effect, B*C"=c(0,0,0,0,0,-1,1,0), "fac2-effect, B*D"=c(0,0,0,0,0,-1,0,1), "fac2-effect, C*D"=c(0,0,0,0,0,0,-1,1)) summary(glht(mod,c2)) Kay Cichini wrote: > &g

Re: [R] interaction contrasts

2010-10-15 Thread Kay Cichini
e point me there. yours, kay Kay Cichini wrote: > > 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=

[R] interaction contrasts

2010-10-13 Thread Kay Cichini
ac2, B*D" ## "fac2, C*D" ## to the simultanous tests to see whether the effects ## of fac2 within the levels of fac1 differ between ## each combination of the levels of fac1, or not ?? thanks for any advise! yours, kay - Kay Cichini Postgraduat

Re: [R] general question on binomial test / sign test

2010-09-02 Thread Kay Cichini
s,0.5,alternative="greater"). then this one of many tests for a sparse dataset came up where x=0, and n=1). there i thought the H0 is x is less than 0.5, and i then had my trouble interpreting the p-value of 1. best, kay Kay Cichini wrote: > > > thanks a lot for the elaborations

[R] general question on binomial test / sign test

2010-09-01 Thread Kay Cichini
best, kay - -------- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/general-question-on-binomial-test-sign-test-tp2419965p2419965.html Sent from the R help ma

Re: [R] lattice help required

2010-08-31 Thread Kay Cichini
re suited there should be only two - one for each response. can someone tell me how i have to change the code? best, kay > Kay, > > doe this do what you want? > > > dotplot(y1+y2 ~ facs$Treatment|facs$

Re: [R] lattice help required

2010-08-25 Thread Kay Cichini
___ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > - Kay Cic

Re: [R] lattice help required

2010-08-25 Thread Kay Cichini
> > > dotplot(y1+y2 ~ facs$Treatment|facs$Sites, > outer=TRUE, > scales = list(x = list(rot = 90, tck=c(1,0))), > ylab=c("y1", "y2"), > xlab=c("Site 1", "Site 2"), > strip=FALSE) > > > On Wed,

[R] lattice help required

2010-08-25 Thread Kay Cichini
nt(pl1, split=c(1,2,1,2), more=TRUE) print(pl2, split=c(1,1,1,2)) # but as said, ideally the plots should be stacked with only the lower plot giving the x-axis annotation and only the upper plot with text-panels. thanks

Re: [R] tick marks on both sides of axis

2010-08-24 Thread Kay Cichini
thanks a lot. sorry - i don't know what you mean by "restoring broken email chain", could not get an answer by re-reading the posting guide, either. yours, kay - ------------ Kay Cichini Postgraduate student Institute of Botany Uni

Re: [R] tick marks on both sides of axis

2010-08-24 Thread Kay Cichini
=F,tcl=-0.25) yours, 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-tp2336668p2336751.html Sent from the R help mailing l

[R] tick marks on both sides of axis

2010-08-24 Thread Kay Cichini
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

Re: [R] logistic regression tree

2010-08-22 Thread Kay Cichini
dear all, thank you everyone for the profound answers and the needful references! achim, thank you for the very kind offer!! sorrily i'm not around vienna in the near feature, otherwise i'd be glad to coming back to your invitation. yours, kay - -------- K

Re: [R] logistic regression tree

2010-08-20 Thread Kay Cichini
ifelse(dat$fac13[i]=="Ia"|dat$fac13[i]=="Ib"|dat$fac13[i]=="Ic", dat$Y[i]<-rbinom(1,1,0.75), dat$Y[i]<-rbinom(1,1,0)) } dat$Y<-as.factor(dat$Y) tr<-ctree(Y~fac1+fac2+fac3,data=dat) plot(tr) ########

Re: [R] logistic regression tree

2010-08-20 Thread Kay Cichini
tively. at the moment i don't know which this would be. regarding the meaning of absences at stands: as these species are frequent in the area and hence there is no limitation by propagules i guess absence is really due to unfavourable conditions. thanks a lot, kay - ---------

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Kay Cichini
..gabor gave the solution during i was typing - so please disregard this. yours, kay - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/paired-samples

Re: [R] paired samples, matching rows, merge()

2010-08-20 Thread Kay Cichini
14 [4,] 5 25 16 ...but wih the irregulary duplicated values in k1 i dead-ended and i guess it could get tricky to solve this. greetings, kay - ------------ Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in

[R] logistic regression tree

2010-08-19 Thread Kay Cichini
, kay - ---- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/logistic-regression-tree-tp2331847p2331847.html Sent from the R help mailing list archive at Nabbl

Re: [R] how to paste a vector to expression?

2010-08-18 Thread Kay Cichini
..of course that works, too. thanks a lot, kay Zitat von "Gavin Simpson [via R]" : > > > > On Tue, 2010-08-17 at 13:49 -0700, Kay Cichini wrote: >> ..thank you! >> >> this generally works well - but after saving the result as jpeg i noticed >>

Re: [R] how to paste a vector to expression?

2010-08-17 Thread Kay Cichini
ion are much smaller). maybe there is still another solution? - really no way using the paste command? curiously paste(expression(bar(x)),"=",dat) does not work.. thanks, kay ----- ---- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -

[R] how to paste a vector to expression?

2010-08-17 Thread Kay Cichini
hello, i need to annotate a plot with an expression and a variable value - like: plot(1:4) dat<-1:2 text(2:3,2:3,expression(paste(bar(x)==dat))) ... but with the dat values plotted, 1 at x=2,y=2, and 2 at x=3,y=3. maybe someone can help? yours, kay - ----

[R] two sample permutation test by sample(..., replacement = T)

2010-07-23 Thread Kay Cichini
thout replacement? thanks a lot for any hints, kay - ---- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/two-sample-permutation-test-by-sample-replacem

Re: [R] stripplot, lattice

2010-05-27 Thread Kay Cichini
...layout=..., of course. - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/stripplot-lattice-tp2232908p2232959.html Sent from the R help mailing list

[R] stripplot, lattice

2010-05-27 Thread Kay Cichini
D7NW1 4 D D7 D7SO1 5 D D7 D7SW 1 5 - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/stripplot-latt

Re: [R] contrasts for lmer model

2010-05-19 Thread Kay Cichini
4.78e-05 *** 3 == 0 -0.528470.09133 -5.786 2.88e-08 *** 4 == 0 -0.443960.15549 -2.855 0.0171 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Adjusted p values reported -- single-step method) ----- ---- Kay Cichini Postgraduate student Institute

Re: [R] contrasts for lmer model

2010-05-19 Thread Kay Cichini
"gap" within each stage simultanously so that the p-values will be adjusted accordingly, isn't it? is this possible? - if so, how would this be done? anyone who can help with this.. best, kay - Kay Cichini Postgraduate student Institute of Botany Univ.

[R] contrasts for lmer model

2010-05-19 Thread Kay Cichini
, 2L, 1L, 4L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 3L, 1L, 7L, 4L, 2L, 4L, 2L, 2L, 3L, 1L, 3L, 1L, 1L, 1L, 2L, 2L, 3L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 0L, 3L, 2L, 2L, 3L, 1L, 2L, 2L, 4L, 5L, 4L, 5L, 7L, 6L, 4L, 9L, 5L, 2L, 2L, 0L, 3L, 6L, 1L, 4L, 3L, 3L, 4L, 2L, 5L, 7L, 4L,

Re: [R] post-hoc t-tests, for lme / lmer models

2010-05-10 Thread Kay Cichini
no one? any pointers would really be greatly appreciated! thanks, kay - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/post-hoc-t-tests-for-lme-lmer

[R] post-hoc t-tests, for lme / lmer models

2010-05-07 Thread Kay Cichini
not seem to be valid - at least i did not find any threads where it was applied for such models. can anyone give some advise on this? thanks, kay - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this

[R] confidence interval for small non iid sample

2010-05-06 Thread Kay Cichini
- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/confidence-interval-for-small-non-iid-sample-tp2132804p2132804.html Sent from the R help mailing list archive at Nabble.com

[R] how to extract from obj. created by plot function stressplot without printing plot?

2010-05-06 Thread Kay Cichini
, kay - ---- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/how-to-extract-from-obj-created-by-plot-function-stressplot-without-printing-plot-tp2132366p2132366.html Sent from the R help ma

Re: [R] superscript

2010-05-04 Thread Kay Cichini
thanks a lot! - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/superscript-tp2125341p2125386.html Sent from the R help mailing list archive at

Re: [R] superscript

2010-05-04 Thread Kay Cichini
thanks a lot! - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/superscript-tp2125341p2125384.html Sent from the R help mailing list archive at

[R] superscript

2010-05-04 Thread Kay Cichini
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 ---

[R] mblm, confint problem

2010-05-03 Thread Kay Cichini
d) # gives: > mod$coeff (Intercept) x 0.3187684 0.6324555 > > confint(mod) 0.0250.975 (Intercept) -0.02332502 0.5550533 x0.77481683 1.0266857 (?) ----- ---- Kay Cichini Postgraduate student

[R] model selection, lme

2010-04-28 Thread Kay Cichini
h appreciate help with this. it may be a trivial problem - but though i could not figure out. thanks in advance, kay - -------- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r

Re: [R] glmer with non integer weights

2010-04-23 Thread Kay Cichini
hello, krebs (1995) states MH as prob., but yes it's rather a ratio of probs. at each site i had 4 blocks with 2 treatments (treat vs. control) - after treating i looked for similarity between each of those pairs. it is of interest if changes in similarity due to treatment differ between stages

Re: [R] multiple plots problem

2010-04-20 Thread Kay Cichini
hi thomas, thanks a lot- of course that's it. i knew why the middle plot is stretcht, but didn't figure out that settting par(oma) is the key.. greetings, kay -- View this message in context: http://n4.nabble.com/multiple-plots-problem-tp2017326p2017551.html Sent from the R help mailing list

[R] multiple plots problem

2010-04-20 Thread Kay Cichini
hello, i try to plot 3 graphs which have the same x.axis underneath each other. i'd like the plots to be aligned without margings between the boxes and draw a single x axis beneath the lowest plot. i managed to get the alignment by setting par(mar), but the middle box is stretched and i cant't

Re: [R] glmer with non integer weights

2010-04-20 Thread Kay Cichini
hello, it's the Morisita Horn Index, which is an ecological index for similarity between two multivariate objects (vegetation samples with species and its abundance) where a value of one indicates completely same relative importance of species in both samples and 0 denotes total absence of any sa

Re: [R] glmer with non integer weights

2010-04-19 Thread Kay Cichini
hi emmanuel, thanks a lot for your extensive answer. do you think using the asin(sqrt()) transf. can be justified for publishing prurpose or do i have to expect criticism. naivly i excluded that possibility, because of violated anova-assumptions, but if i did get you right the finite range rathe

Re: [R] glmer with non integer weights

2010-04-16 Thread Kay Cichini
thank you thomas for the helpful hint! yours, kay -- View this message in context: http://n4.nabble.com/glmer-with-non-integer-weights-tp1837179p1965827.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing lis

Re: [R] glmer with non integer weights

2010-04-16 Thread Kay Cichini
thanks thierry, i considered this transformations already, but variance is not stabilized and/or normality is neither achieved. i guess i'll have to look out for non-parametrics? best regards, kay -- View this message in context: http://n4.nabble.com/glmer-with-non-integer-weights-tp1837179p19

r-help@r-project.org

2010-04-13 Thread Kay Cichini
hello heinz, my mistake, i thought this threads dealt only with tests for 2-level grouping factor. many thanks for the hint! kay -- View this message in context: http://n4.nabble.com/kruskal-wallis-stratified-tp1838210p1838666.html Sent from the R help mailing list archive at Nabble.com. _

[R] standard errors, graphical representation of glmer

2010-04-13 Thread Kay Cichini
hi, i did a binomial glmer for a repeated measure design. now i wanted to use some indicator of variance for the fixed eff. of the model for the graphical representations. sd's of the plain incident rates are huge and misleading in that context. thus i thought of extracting se's of the coefficien

Re: [R] glmer with non integer weights

2010-04-13 Thread Kay Cichini
thanks thierry, my problem is that the index is a propability which is not derived from incidents per nr. of observations, thus i don't have those numbers but only the plain index, which i want to test. greatings, kay -- View this message in context: http://n4.nabble.com/glmer-with-non-integer

r-help@r-project.org

2010-04-13 Thread Kay Cichini
hello heinz, i read the thread already. i think it applies only to 2-sample problems. greatings, kay -- View this message in context: http://n4.nabble.com/kruskal-wallis-stratified-tp1838210p1838261.html Sent from the R help mailing list archive at Nabble.com.

[R] kruskal-wallis, stratified

2010-04-13 Thread Kay Cichini
hello everyone, can anybody tell me if there is a kruskal-wallis, or another non-parametric test, that can deal with multiple samples that are stratified? thanks, kay -- View this message in context: http://n4.nabble.com/kruskal-wallis-stratified-tp1838210p1838210.html Sent from the R help mai

[R] glmer with non integer weights

2010-04-12 Thread Kay Cichini
hello, i'd appreciate help with my glmer. i have a dependent which is an index (MH.index) ranging from 0-1. this index can also be considered as a propability. as i have a fixed factor (stage) and a nested random factor (site) i tried to model with glmer. i read that it's possible to use a quasib

[R] graphical data representation after glmm, glmer

2010-04-12 Thread Kay Cichini
hello, i ran a glmer for binomial data, with one (ordered) between and one within subjects factor. now i was thinking of how to plot my dependent. i plotted the probablities (p=X(incidents)/n(observations)) for each group of the between factor in a line chart, with its corresponding standard dev

Re: [R] bootstrap confidence intervals, non iid

2010-04-09 Thread Kay Cichini
hi glen, i need conf.intervals for blocked data, as described in the first place. i've learned in the meantime, that the boot() function can handle this. i had to formulate the function for the boot command, put "sites" to the strata argument and resample from each subsetted level of the fact

Re: [R] bootstrap confidence intervals, non iid

2010-04-06 Thread Kay Cichini
not yet, but it would either not solve the non iid problem, or would it? -- View this message in context: http://n4.nabble.com/bootstrap-confidence-intervals-non-iid-tp1751619p1752583.html Sent from the R help mailing list archive at Nabble.com. __ R

Re: [R] NMDS Ordination Graphics Problem

2010-04-06 Thread Kay Cichini
hi, maybe you should see: http://cc.oulu.fi/~jarioksa/softhelp/FAQ-vegan.html#How-to-use-different-plotting-symbols-in-ordination-graphics_003f greetings, kay -- View this message in context: http://n4.nabble.com/NMDS-Ordination-Graphics-Problem-tp1751845p1752557.html Sent from the R help mail

[R] bootstrap confidence intervals, non iid

2010-04-05 Thread Kay Cichini
hello, i need to calculate ci's for each of 4 groups within a dataset, to be able to infere about differences in the variable "similarity". the problem is that data within groups is dependent, as assigned by the blocking-factor "site". my guess was to use a block bootstrap but samples within in t

Re: [R] getting CI's for certain y of nls fitted curve

2010-03-30 Thread Kay Cichini
...it's of course simply using the desired x in the predict function. in this case: predict(mod1,data.frame(press = x_tenth[1]). it must have been a trivial syntax error, why this didn't work in the first place. kay -- View this message in context: http://n4.nabble.com/getting-CI-s-for-certai

Re: [R] inverse function of SSlogis

2010-03-29 Thread Kay Cichini
hello gabor, thanks a lot for that. i'll try. best wishes, kay -- View this message in context: http://n4.nabble.com/inverse-function-of-SSlogis-tp1694995p1695061.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

[R] getting CI's for certain y of nls fitted curve

2010-03-29 Thread Kay Cichini
hello, i managed to get CI's for my curve - but now I need the intervall for a certain y point (y_tenth) of the curve.. can anyone help me with this? #data: por<-data.frame(list(structure(list(run = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), .Label = c("1", "3", "4"), class

[R] inverse function of SSlogis

2010-03-29 Thread Kay Cichini
hello, is there a more convenient way to get the inverse of SSlogis than transforming the function by hand? thanks, kay -- View this message in context: http://n4.nabble.com/inverse-function-of-SSlogis-tp1694995p1694995.html Sent from the R help mailing list archive at Nabble.com. ___

[R] predict

2010-03-26 Thread Kay Cichini
hello, sorry to have hassled you with this - in the end it is all quite simple and after digging deeper into the help files, also of predict.lm, i finally managed to do it. here is the worked example with solution for anyone beeing clumsy as i was: ###

[R] nls, predict() problem

2010-03-25 Thread Kay Cichini
hello, can anyone help with this: ### ###data: measurments (response = trans) run several times at the same predictor value level (press) por<-data.frame(list(structure(list(run = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L

Re: [R] Odp: boxplot, vertical position of x-axis labels

2010-03-17 Thread Kay Cichini
thanks a lot! both ways work. greetings - kay -- View this message in context: http://n4.nabble.com/boxplot-vertical-position-of-x-axis-labels-tp1594857p1596456.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org m

[R] boxplot, vertical position of x-axis labels

2010-03-16 Thread Kay Cichini
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 F

[R] lme for repeated measures, one within, one between factor

2010-03-07 Thread Kay Cichini
hello list, the topic is covered extensively but from none of the postings i could conclude the correct statement for my design: a 2-level within and a 2-level between subjects factor, both fixed, subjects as random factor. i want to test wheter there is a within effect and if it is different fo

Re: [R] adonis(), design

2010-03-04 Thread Kay Cichini
inclusion of repeated measures of course should gain power, but here I guess one would have to restrict permutations and that's what may reduce power drastically if sample size is small, at least that's how I understood it. I have also dug out a thread where someone asked for random factors in

Re: [R] repeated measures anova, car package

2010-03-03 Thread Kay Cichini
sorry, it should say "I copied the model from the car manual.." more precise, I did the same as there in doing the lm(cbind()) call first an then the av.ok, etc. -- View this message in context: http://n4.nabble.com/repeated-measures-anova-car-package-tp1573721p1577432.html Sent from the R h

Re: [R] repeated measures anova, car package

2010-03-02 Thread Kay Cichini
Hello John, As you said, I could also take a means model and test linear hypothesis for the desired effects - would this also be the case for the repeated measure i did in the first place. I copied the model from the car model where you first call: > modx<-lm(cbind(div_h, div_l) ~ site, divrep

[R] repeated measures anova, car package

2010-03-01 Thread Kay Cichini
level, but I guess that's not the best solution. Thanks for any advise, Kay Cichini -- View this message in context: http://n4.nabble.com/repeated-measures-anova-car-package-tp1573721p1573721.html Sent from the R help mailing list archive at Nabble.com. __

[R] adonis(), design

2010-02-25 Thread Kay Cichini
hello list, i'd appreciate help regarding formula expression in adonis(), package vegan. i'm interested in the effect of an impact on species composition and if this effect is the same for different sites. i have an impact survey with plot pairs, one plot with impact the other with similar ecolg

Re: [R] pvclust() with distance method from vegdist(), package = vegan

2010-02-03 Thread Kay Cichini
thank you for the information. would be great to have the possibility to check significance of clusters coming from other distance-measures, too - but i'll be satisified with what pvclust() can do. yours, kay -- View this message in context: http://n4.nabble.com/hvcluster-with-distance-method-f

Re: [R] hvcluster() with distance method from vegdist(), package = vegan

2010-02-02 Thread Kay Cichini
sorry, it should say pvclust(), of course... -- View this message in context: http://n4.nabble.com/hvcluster-with-distance-method-from-vegdist-package-vegan-tp1459859p1459865.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pr

[R] hvcluster() with distance method from vegdist(), package = vegan

2010-02-02 Thread Kay Cichini
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"

Re: [R] restricted permutations in permutest()?

2010-01-20 Thread Kay Cichini
...original data and design were posted at r-forge vegan forum today. greetings, kay -- View this message in context: http://n4.nabble.com/restricted-permutations-in-permtest-tp1017422p1018531.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] restricted permutations in permutest()?

2010-01-19 Thread Kay Cichini
i'll post in r-forge vegan help forum and appreciate your help very much. greetings, kay -- View this message in context: http://n4.nabble.com/restricted-permutations-in-permtest-tp1017422p1017865.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] restricted permutations in permutest()?

2010-01-19 Thread Kay Cichini
hello gavin, you are right, i didn't get into the documentation to deep and i'm also a beginner, that's why i'm just about to get into the logical part of the syntax. now, the output from perm.disp() says: #No. of permutations: 999 #Permutation type: free #Permutations are unstratified #Mirr