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&
}
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
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
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
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
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
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
"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
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
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]]
_
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
-----
----
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
,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
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=
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
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
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 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$
___
> 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
>
>
> 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,
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
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
=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
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
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
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)
########
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
-
---------
..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
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
,
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
..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
>>
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
-
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
-
----
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
...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
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
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
"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.
, 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,
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
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
-
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
,
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
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
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
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
---
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
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
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
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
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
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
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
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
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
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.
_
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
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
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.
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
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
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
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
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
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
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
...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
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
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
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.
___
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:
###
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
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
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
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
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
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
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
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.
__
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
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
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
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"
...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.
___
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.
___
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
86 matches
Mail list logo