.bib file: https://www.dropbox.com/s/r1yjqx0sni2pzjy/quasi.bib?dl=0
John Maindonald email: john.maindon...@anu.edu.au
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE d
No, it is not distribution free. Independent random sampling is assumed.
That is a non-trivial assumption, and one that is very often not true or not
strictly true.
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
On 21/03/2021, at 00:
. I do not know what Anamika’s
position is, but he may well be trying to build a business from a
much less privileged starting point than others of us may have
enjoyed. In any case, the comment strikes a tone that is out of
place in respectful discourse.
John Maindonald email
dom number generator),
small P-values shift the weight of evidence. An alternative that is apriori
highly unlikely takes a lot of shifting.
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
On 3/04/2016, at 22:00,
r-help-requ...@r-pr
What’s more, for pdf output one can use R Markdown and judiciously
sneak in html and/or LaTeX (consider however what the processing
steps might do to such markup).
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
On 19/11/2015, at 00:00,
estimated proportions
for binomial data. For the poisson or quasipoisson family, use of the sqrt link
avoids this problem.
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
On 28/07/2015, at 22:00,
r-help-requ...@r-project.org<mailto:r-help-
.
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
On 22/07/2015, at 22:00,
r-help-requ...@r-project.org<mailto:r-help-requ...@r-project.org> wrote:
Da: li...@dewey.myzen.co.uk<mailto:li...@dewey.myzen.co.uk>
Data: 2
9.899249 -1.65842 0.03797342
## Eureka!
## cf also
as.vector(pseudoinverse(t(Xcs0) %*% Xcs0) %*% t(Xcs0) %*% y)
## [1] -1.65842037 0.03797342
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
Centre for Mathematics & It
Thanks for that. Useful to have that question asked and to get that
information.
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
On 25/03/2015, at 0:00,
r-help-requ...@r-project.org<mailto:r-help-requ...@r-project.org> wrote:
monotonically with the number of splits).
The rel error column from printcp() can be hopelessly optimistic.
John Maindonald email:
john.maindon...@anu.edu.au<mailto:john.maindon...@anu.edu.au>
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics
nd y are local to the individual panel
NB also that layer() accepts a data argument. The following is an
alternative way to calculate one4all:
one4all <- layer(data=xdat, avy <- median(mortality),
panel.segments(0.1, avy, 0.3, avy, col="red",lwd=4),
panel.segments(
ph+addlayer
The code that produces the object gph would also be simpler
and easier to follow if some relevant part was separated out into
a separate layer.
See also my notices on layering of lattice objects at:
http://www.maths.anu.edu.au/%7Ejohnm/r-book/add-graphics.html
John Maindonald
ocial science
researchers and users of research.
"
The title for the planned session is:
"The R System as a Platform for Analysis and Development of Analysis
Methodology"
http://conference.acspri.org.au/index.php/rc33/2012/schedConf/trackPolicies
John Maindonald
It means that dat$v1, dat$v2, . . . are not columns in the data frame df
Specify the formula as v1 ~ v2+v2+v3+v4+v5+factor
Then (assuming that factor is a column object of the right length)
you should be fine.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473
You might look at the function confusion() in the DAAGxtras package.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Buildin
and
test.tex, irrespective of the 1 line or 2 issue. The first of the
\SweaveOpts lines had an effect, unless I am missing something.
John.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applicat
mysteriously appears in the second
line of code in output .tex file test2.tex
John.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building
Yes, that has fixed the problem. (2010-11-24 r53659)
Thanks.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Building 27)
Austra
-8 reshape_0.8.3
[7] plyr_1.2.1 proto_0.3-8
loaded via a namespace (and not attached):
[1] ggplot2_0.8.8 latticeExtra_0.6-14
Is there a workaround?
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
C
ssion(levlist), ...){
panel.text(x=0, y=0.5,
lab = as.expression(levlist[which.panel[which.given]]), adj=0)
}
e.g.
levlist <- list("A", "B")
xyplot(11:14 ~ 1:4 | rep(1:2,2), scales=list(x=list(alternating=c(1,1),
relation="sliced")), strip=stripfun, layou
e collection and
sharing of knowledge on animals and their environments for zoos, aquariums and
related organizationsvalues the use of objective data to benefit conservation,
science, animal welfare, education, and collection management."
John Maindonald email: john.maindon...@
at <- locatefun(tx)
backn <- regexpr("\\n",entrymat[,1],fixed=TRUE)
entrymat <- entrymat[backn < 0,]
entrymat[,2] <- entrymat[,2] - offset
entrymat[,1] <- gsub("_","\\_",entrymat[,1], fixed=TRUE)
nmatch <- matc
I had a large job time ago that ran fine under MacOS X.
I'd expect the same to be true under Linux. It would run
under Windows XP only if XP had been freshly rebooted.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centr
I am looking for information on experimentation with the use
of R in the teaching of statistics and science in schools. Any
leads would be very welcome. I am certain that there is such
experimentation.
I've made this inquiry on r-sig-teaching, with no response.
John.
John Maind
Please add
AUSTRALIA
Canberra: Canberra R Users Group (CANRUG) http://canrug.togaware.com/
Cheers
john.
John Maindonald email: john.maindon...@anu.edu.au
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John De
u, sep=""),
collapse="+")))
xyplot(form, data=df, outer=TRUE)
I'd hoped that the following would do the trick, but the first label
is repeated in each panel, and the variable names are added:
xyplot(form, data=df, outer=TRUE, strip=strip.c
e check of accuracy really has
then to be tuned to that application. Results in cases where the
"accuracy" beyond (optimistically) the sixth digit is of no
consequence are unlikely to give much clue on performance in such
specialist applications.
John Maindonald
es
at http://www.maths.anu.edu.au/%7Ejohnm/r-book/2edn/xtras/xtras.pdf
have brief comments, and code that can be used to illustrate the point.
I support Stephen Kolassa's suggestions re using simulation for
sensitivity analysis, though I think this can also be done analytically.
John Maindonald
e unfair to Ross and Robert, but no-one
complained. Perhaps no-one ever read that far through the preface!
There's an excellent brief summary of the history of R, and its
connections with S, in Section 1.4 of John Chambers' "Software for
Data Analysis".Appe
It would, I consider, be useful to note, on the help pages for Sweave
and RweaveLatex, that the path to Sweave.sty in the R installation can
be extracted as paste(R.home(),"share/texmf/", sep="/") or its
equivalent.
John Maindonald email: john.maindon...@an
; or suchlike package to which such datasets can be
submitted? All it would require is for someone to offer to act as
Keeper of the Miscellaneous Data".
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & I
land.ac.nz/uoa/about/news/articles/2008/11/rsnz_honours.cfm
12 November 2008.
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Building
2001" = 39, "2002" =
119,
"2003" = 276, "2004" = 523,"2005" = 945,"2006" = 1475, "2007" =
2015,
"2008"=1644)
These will not be all that accurate; there will be omissions
and duplications.
Growth is close to e
s=269+108, zoology=267+111,
plantSciences=250+108, biochemANDmolbio=229+200,
mathANDcompBIO=224+143,
biotechANDappliedmicrobiology=223+159,
evolutionaryBIO=210+117)
There's a great deal more summary information that might
be extracted. What is a good way, with readily av
later; why not do some
smoothing at the start?
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Building 27)
Australian National University
/machine random effects be independent,
but allows them to have different variances. I am wondering whether
it is possible to fit such a model using lmer().
[In this example the large estimated correlations suggest that it is not
a sensible model.]
John Maindonald email: [EMAIL PROTECTE
}
## Usage
glmm2 <- lmer(rcr ~ litter + (1 | Farm), family=binomial,
data=data.site)
ciplot(obj=glmm2)
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical
allelR,
visit
www.revolution-computing.com"
At present though, unless you have better success than me, all you
will get from the Revolution Computing site is a revolutionary song!
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for M
* .(x) * "," * .(y) * "]"),
zz[-length(zz)], zz[-1], SIMPLIFY=FALSE)
y <- rnorm(40)
x0 <- cut(y, breaks=breaks)
plot(x0, y, xaxt="n")
axis(1, at=1:3, labels=as.expression(lbl))
Thanks, Peter.
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (612
# Can alternatively use bquote()
axis(1, at=1:3, labels=as.expression(axlab), line=2.25, lty=0)
However I have been unable to modify axlab so that the infinity
symbol appears in place of 'Inf'. Is there is some relatively
straightforward way to do this? The issue is of course more
gener
d
to get an answer.
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Building 27)
Australian National University, Canberra ACT 0200.
On 27 Dec
n
extended paper.
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Building 27)
Australian National University, Canberra ACT 0200.
On
ing to
sea in a rowing boat to find out whether conditions are sufficiently
calm for an ocean liner to leave port".
John Maindonald email: [EMAIL PROTECTED]
phone : +61 2 (6125)3473fax : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Math
rt R in a working directory whose
.RData file, if any, has the proper R-2.6.0 pedigree, and source
the dumped file. I do not want to have to directly run the function
that creates the database each time I start a session.
Has anyone else hit this?
John Maindonald email: [EMAIL PRO
44 matches
Mail list logo
