# With index p between 1 and 2, this produces continuous
# data with exact zeros
x <- rnorm( length(y), 0, 1) # Unrelated predictor
# With exact zeros, index p must be between 1 and 2
# Fit the tweedie distribution; expect p about 1.5
out <- tweedie.profile( y~1, p.vec=se
year 2010
month 10
day15
svn rev53317
language R
version.string R version 2.12.0 (2010-10-15)
Peter Dunn: Biostatistician (Room T4.15
utput is the same as for termplot(m1) produced, and I
don't really learn anything new.
So the original questions remain: If this is intended, what is the logic? Or am
I missing something?
Perhaps I could add: How do I explain and justify this behaviour to someone?
Thanks.
P.
Peter Dunn
win9.8.0
locale:
[1] en_AU.UTF-8/en_AU.UTF-8/C/C/en_AU.UTF-8/en_AU.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
Peter Dunn: Biostatistician (Room T4.12)
School of Health and Sport Science
Faculty of Science, Health and Education ML-34
C;LC_MEASUREMENT=en_AU.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] GLMsData_0.7
loaded via a namespace (and not attached):
[1] tools_2.9.0
>
Peter Dunn: Biostatist
tex
This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6)
%&-line parsing enabled.
**
And it waits for input.
So that doesn't seem to be the problem.
Further suggestions?
P.
Peter Dunn: Biostatistician
School of Health and Sport Science
Faculty of Science, Health and Educ
/enc/dvips/base/8r.enc}
Output written on GLMsData-manual.pdf (79 pages, 405091 bytes).
Transcript written on GLMsData-manual.log.
pdu...@pdunnubuntu:~/DSdata/GLMsData.Rcheck$
--
Peter Dunn: Biostatistician
School of Health and Sport Science
Faculty of Science, Health and Education ML-
ove choices are easily adopted.
If no easy solutions exist or emerge, I'm happy to run with \url.
Thanks again.
P.
Peter Dunn
Biostatistician
School of Health and Sport Science
Faculty of Science, Health and Education
University of the Sunshine Coast
Tel: +61 7 5456 5085
Fax: +61 7
it will never changed.
(By my definition, the square of the deviance residuals is not quite the unit
deviance as the weights need to be incorporated differently, but others may
disagree.)
Of course, I could be way off track. I would appreciate enlightenment.
P.
Peter Dunn
Biostatistician
about it? That makes me quite
angry...
P.
Peter Dunn
Biostatistician
School of Health and Sport Science
Faculty of Science, Health and Education
University of the Sunshine Coast
Tel: +61 7 5456 5085
Fax: +61 7 5430 2896
Email: [EMAIL PROTECTED]
www.usc.edu.au
Peter Dunn
Biostatistician
Scho
phics grDevices utils datasets methods base
loaded via a namespace (and not attached):
[1] tools_2.7.0
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Sciences, USQ; http://www.sci.usq.edu.au/staff/dunn
Aust. Centre for Sustainable Catchments: www.usq.edu.au/acsc
This email (including an
nderstand!)
Also, since the above obviously doesn't do what I hoped
(consistently printing two sig figs), could someone also explain
how I can consistently get two significant figures in situation
like above?
Thanks as always.
P.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Scie
nderstand!)
Also, since the above obviously doesn't do what I hoped
(consistently printing two sig figs), could someone also explain
how I can consistently get two significant figures in situation
like above?
Thanks as always.
P.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Scie
t.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.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Sciences, USQ; http:/
> What is quantile residuals?
See
Dunn, Peter K. and Smyth, Gordon K. (1996). Randomized Quantile
Residuals. Journal of Computational and Graphical Statistics,
volume 5, issue 3, 236–244.
Or see ?qresid in R package statmod.
P.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Sciences,
;=
plot( c(-1,1), c(-1,1), type="n",
main="Graph of Eq (\label{EQ:equation})")
abline(-1, 3)
@
Thanks as always.
P.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Sciences, USQ; http://www.sci.usq.edu.au/staff/dunn
Aust. Centre
e Tweedie distributions have mu
greater than zero. In using the identity link, you may be asking R
to fit the impossible.
Perhaps the question you need to ask is this: Why do you need an
identity link?
P.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Sciences, USQ; http://www.sci.usq.edu
tions
exist such as hurdle models, Tweedie glms (package tweedie) and so
forth.
P.
--
Dr Peter Dunn | dunn usq.edu.au
Faculty of Sciences, USQ; http://www.sci.usq.edu.au/staff/dunn
Aust. Centre for Sustainable Catchments: www.usq.edu.au/acsc
This email (including any attached fi
=English_Australia.1252;LC_NUMERIC=C;LC_TIME=English_Australia.1252
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "datasets"
"methods" "base"
other attached packages:
car xtable statmod
&
R version 2.5.0 (2007-04-23)
i486-pc-linux-gnu
locale:
LC_CTYPE=en_AU.UTF-8;LC_NUMERIC=C;LC_TIME=en_AU.UTF-8;LC_COLLATE=en_AU.UTF-8;LC_MONETARY=en_AU.UTF-8;LC_MESSAGES=en_AU.UTF-8;LC_PAPER=en_AU.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_AU.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1
20 matches
Mail list logo