[R] syntactical meaning of fullstop in R functions

{ stop(name, not recognized as a valid name.) } } bytecode: 0x10f6e5ac environment: namespace:raster thanks in advance, kay -- Kay Cichini, MSc Biol Grubenweg 22, 6071 Aldrans Tel.: 0650 9359101 E-Mail: kay.cich...@gmail.com Web: www.theBioBucket.blogspot.co.athttp

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

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 kay.cich...@gmail.com hello, what is the meaning of a fullstop in the below R-Ffunction (like in '.GADM') library(dismo) Lade

Re: [R] Testing continuous zero-inflated response

Many thanks - this was very helpful! Regards, Kay Am 28.01.2013 13:19 schrieb Achim Zeileis achim.zeil...@uibk.ac.at: On Sun, 27 Jan 2013, Kay Cichini wrote: That said, wilcox_test(x ~ factor(y), distribution = exact) or the same with oneway_test, i.e would be ok? Yep, exactly

Re: [R] Testing continuous zero-inflated response

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 achim.zeil...@uibk.ac.at On Fri, 25 Jan 2013, Kay Cichini wrote: Hello, I'm searching for a test that applies

Re: [R] Testing continuous zero-inflated response

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 achim.zeil...@uibk.ac.at On Sun, 27 Jan 2013, Kay Cichini wrote: Thanks for the reply! Still, aren't there issues with 2-sample test vs y and excess zeroes

[R] Testing continuous zero-inflated response

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...@gmail.com

Re: [R] bootstrapping quantile regression

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 kay.cich...@gmail.com: HI everyone, I try to get some bootstrap CIs for coefficients obtained by quantile regression. I have influencial values

[R] bootstrapping quantile regression

+ P.Form) # get 99% confidence intervals boot.ci(results, type=bca, index=1, conf = .99) # intercept boot.ci(results, type=bca, index=2, conf = .99) # P.Perimeter boot.ci(results, type=bca, index=3, conf = .99) # P.Form -- Kay Cichini, MSc Biol Grubenweg 22, 6071 Aldrans E-Mail: kay.cich...@gmail.com

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

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

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

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

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

[R] Package/Function for Blending Images

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

=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: hello, i was shortly asking the list for help with some interaction contrasts (see below) for which i had to change

Re: [R] interaction contrasts

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=LETTERS[1:4]), fac2=rep(c(I,II),16),y=rnorm(32,1,1

[R] interaction contrasts

of the levels of fac1, or not ?? thanks for any advise! yours, kay - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/interaction-contrasts

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

,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. your explanations clearly

[R] general question on binomial test / sign test

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

Re: [R] lattice help required

, doe this do what you want? 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, Aug 25, 2010 at 11:04 AM, Kay Cichini kay.cich

[R] lattice help required

)) # 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 a lot, kay - Kay Cichini Postgraduate student

Re: [R] lattice help required

= list(rot = 90, tck=c(1,0))), ylab=c(y1, y2), xlab=c(Site 1, Site 2), strip=FALSE) On Wed, Aug 25, 2010 at 11:04 AM, Kay Cichini kay.cich...@uibk.ac.atwrote: hello, i want to stack two lattice plots beneath each other using one x-axis and sharing the same text

Re: [R] lattice help required

]] __ 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 Cichini Postgraduate student

[R] tick marks on both sides of axis

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

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

- 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 list archive at Nabble.com

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

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 Univ. of Innsbruck

Re: [R] logistic regression tree

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

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

...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 context: http://r

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

..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] logistic regression tree

. 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 - Kay Cichini

Re: [R] logistic regression tree

, 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) ## - Kay Cichini Postgraduate student Institute

[R] logistic regression tree

- 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 Nabble.com

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

..of course that works, too. thanks a lot, kay Zitat von Gavin Simpson [via R] ml-node+2329504-916404301-40...@n4.nabble.com: 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

[R] how to paste a vector to expression?

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

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

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 -- View this message

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

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

[R] stripplot, lattice

D7SW1 5 - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/stripplot-lattice-tp2232908p2232908.html Sent from the R help mailing list

Re: [R] stripplot, lattice

...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] contrasts for lmer model

- Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck -- View this message in context: http://r.789695.n4.nabble.com/contrasts-for-lmer-model-tp682p682.html Sent from the R help mailing list archive at Nabble.com

Re: [R] contrasts for lmer model

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. of Innsbruck -- View this message in context

Re: [R] contrasts for lmer model

-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 of Botany Univ. of Innsbruck

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

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

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

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

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

[R] confidence interval for small non iid sample

- 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

Re: [R] superscript

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

Re: [R] superscript

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

[R] mblm, confint problem

(Intercept) x 0.3187684 0.6324555 confint(mod) 0.0250.975 (Intercept) -0.02332502 0.5550533 x0.77481683 1.0266857 (?) - Kay Cichini Postgraduate student Institute of Botany Univ. of Innsbruck

[R] model selection, lme

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.789695.n4.nabble.com/model-selection

Re: [R] glmer with non integer weights

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

Re: [R] glmer with non integer weights

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

[R] multiple plots problem

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] multiple plots problem

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

Re: [R] glmer with non integer weights

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

Re: [R] glmer with non integer weights

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:

Re: [R] glmer with non integer weights

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

[R] kruskal-wallis, stratified

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

Re: [R] kruskal-wallis, stratifiedhttp://n4.nabble.com/forum/NewNode.jtp?tpl=replynode=1838232

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.

Re: [R] glmer with non integer weights

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:

[R] standard errors, graphical representation of glmer

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

Re: [R] kruskal-wallis, stratifiedhttp://n4.nabble.com/forum/NewNode.jtp?tpl=replynode=1838232

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] graphical data representation after glmm, glmer

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

[R] glmer with non integer weights

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

Re: [R] bootstrap confidence intervals, non iid

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

Re: [R] NMDS Ordination Graphics Problem

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

Re: [R] bootstrap confidence intervals, non iid

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] bootstrap confidence intervals, non iid

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

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

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

[R] inverse function of SSlogis

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] getting CI's for certain y of nls fitted curve

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 =

Re: [R] inverse function of SSlogis

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

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

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,

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

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

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

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

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

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

Re: [R] adonis(), design

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

Re: [R] repeated measures anova, car package

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

Re: [R] repeated measures anova, car package

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

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

[R] adonis(), design

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

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

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:

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

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

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

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

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

...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()?

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

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

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.