On 01/02/15 02:17, John Sorkin wrote:
I am trying to understand the Error function and its use in ANOVA. In
particular I want to understand the difference between two models
that differ only with respect to the Error statement:
aovsubj - aov(value~group+time+Error(subject),data=dataRMANOVA)
Hello Mr. Holtman,
Thank you very much for your reply and suggestion. This is what each Year's
data looks like;
tmp1 - structure(list(FIPS = c(1001L, 1003L, 1005L), X2026.01.01.1 =
c(285.5533142,
285.5533142, 286.2481079), X2026.01.01.2 = c(283.4977112, 283.4977112,
285.0860291),
I am trying to use Poisson regression to model count data with four
explanatory variables: ratio, ordinal, nominal and dichotomous – x1, x2, x3
and x4. After playing around with the input for a bit, I have formed – what
I believe is – a series of badly fitting models probably due to
overdispersion
I run the code below successfully on Mac and Ubuntu successfully.
When I run on Windows, I get the results shown. How do I get the code
to work on Windows? I've googled extensively with no success. Thanks
in advance.
require(twitteR)
Loading required package: twitteR
Loading required package:
Honestly? Did you try rcurl https windows (without the quotes)?
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...
Try
indx - which(!!mat, arr.ind=TRUE)
v1 -unname(sapply(split(indx[,2], indx[,1]),toString))
cat(paste(v1, collapse=\n), sep=\n)
1, 2, 3, 6, 7, 8, 9
1, 2, 3, 6, 8, 9
1, 3, 4, 6, 7, 8, 9
1, 8
1, 3, 6, 7, 8, 9
1, 3, 4, 6, 8, 9
1, 3, 5, 9
A.K.
Hi,
Is there a way to map data from Binary
Dear all:
I have a list like that,which is a standard str_locate_all() function
(stringr package) output:
$K
start end
$GSEGTCSCSSK
start end
[1,] 6 6
[2,] 8 8
$GFSTTCPAHVDDLTPEQVLDGDVNELMDVVLHHVPEAK
start end
[1,] 6 6
$LVECIGQELIFLLPNK
start end
[1,] 4
There are two straightforward ways of modelling overdispersion:
1) Use glm as in your example but specify family=quasipoisson.
2) Use glm.nb in the MASS package, which fits a negative binomial model.
On 1 February 2015 at 16:26, JvanDyne e283...@trbvm.com wrote:
I am trying to use Poisson
Dear group,
I have many plots and numeric results in my R program, kindly how can I save
them all sequently on one file.
thanks in advance
RAI
[[alternative HTML version deleted]]
__
A third, and often preferable, way is to add an observation-level random effect:
library(lme4)
data1$obs - factor(seq_len(nrow(data1)))
model - glmer(y ~ x1 + x2 + (1 | obs), family=poisson(link=log), data=data1)
See http://glmm.wikidot.com/faq and search for individual-level
random effects.
But in addition to what Jeff noted, see ?save and ?save.image
(noting that that the resulting .Rdata file can only be read by R).
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
Data is not information. Information is not knowledge. And knowledge
is certainly not
If you're happy with outputting to a multi-page PDF, then you can just
set the default graphics device to pdf(), i.e.
options(device=pdf)
and the start plotting:
plot(1:10, col=0)
plot(10:1, col=1)
plot((1:10)^2, col=2)
plot((10:1)^2, col=3)
and at the end make sure to close the device:
In general this depends what you plan to do with those results. I suspect you
are looking for something like knitr with rmarkdown (.Rmd files to create HTML
or Word) or LaTeX (.Rnw files to create PDF).
---
Jeff Newmiller
On Feb 1, 2015, at 8:26 AM, JvanDyne wrote:
I am trying to use Poisson regression to model count data with four
explanatory variables: ratio, ordinal, nominal and dichotomous – x1, x2, x3
and x4. After playing around with the input for a bit, I have formed – what
I believe is – a series of
I've got to remember to use more spaces. Here's the basic problem:
These are the same:
v 1
v1
But these are extremely different:
v -1
v-1
This mistake can get you even inside of a function call like this:
v - -2:2
which( v1 )
[1] 1 2 3
which( v-1 ) # oops, I meant v -1 not v-1 a HUGE
Using = has it's problems too.
For example,
print(fit - lm(...))
Assigns the result of the lm call to fit and prints the results. This is
quite a useful trick actually.
print(fit = lm(...))
Throws an error.
Moral of story, computers do what you tell them, not what you meant.
Kevin
On
Mike Miller mbmiller+l at gmail.com writes:
I've got to remember to use more spaces. Here's the basic problem:
These are the same:
v 1
v1
But these are extremely different:
v -1
v-1
This is indeed documented, in passing, in one of the pages you listed:
All the more reason to use = instead of -
-Original Message-
From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Ben Bolker
Sent: Monday, 2 February 2015 2:07p
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] the less-than-minus gotcha
Mike Miller mbmiller+l at gmail.com
Just what is meant by dummy points as referred to by the help for the
deldir() function? I understood they indicated the boundary beyond
which triangulation would cease.
I thought I would need the x/y elements (as described in the help file
at the end of the description of the use of the dpl
Hi
Just a warning. Do you intend it for naming some objects. If yes do not do it,
use list instead. If not, just discard my comment.
Cheers
Petr
-Original Message-
From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Knut
Hansen
Sent: Friday, January 30, 2015 2:34 PM
20 matches
Mail list logo