Displaying 20 results from an estimated 2000 matches similar to: "Suggestions for R-intro manual (PR#10701)"
2010 Dec 06
1
waldtest and nested models - poolability (parameter stability)
Dear All,
I'm trying to use waldtest to test poolability (parameter stability) between
two logistic regressions. Because I need to use robust standard errors
(using sandwich), I cannot use anova. anova has no problems running the
test, but waldtest does, indipendently of specifying vcov or not. waldtest
does not appear to see that my models are nested. H0 in my case is the the
vector of
2003 Mar 04
4
writing several command line in R console
Hi R lovers
I would like to know how to step to the next line in the R console editor
without breaking the continuity of my code
more clearly : if for example I write a function, so far i have to write the
all code inside on the same line wich may become obscure as the function is
more and more complex.
I would like to do like in the example of the manuels:
>twosam <- function(y1, y2) {
2002 Nov 24
1
Understanding function residuals()
Hello:
I am trying to understand why glm() does not replicate the results in
Dobson, "Introduction to Generalized Linear Models," pp. 17-20.
I set up the following model. The variable CONDT is assumed as Poisson and
the objective is to estimate the expected value.
The data (chronic medical conditions among women in Australia) is as
follows:
CONDT <- c(0, 1, 1, 0, 2, 3, 0, 1,
2001 Oct 23
1
FTP-Access from R
I want to access a file on a ftp-server with R. But it doesn't work. I
suppose the reason is the username and password. Here the non working
file path
ftp://woudc:woudc*@ftp.tor.ec.gc.ca/Archive-NewFormat/totalozone_1.0_1/stn035/dobson/1929/19290301.dobson.beck.002.smi.csv
I tested other files on this server, same result: a crash of R. I tested
it with other 'normal' sites (on other
2007 Nov 07
1
mixed model testing
Is there a formal way to prove the need of a mixed model, apart from e.g. comparing the intervals estimated by lmList fit?
For example, should I compare (with AIC ML?) a model with seperately (unpooled) estimated fixed slopes (i.e.using an index for each group) with a model that treats this parameter as a random effect (both models treat the remaining parameters as random)?
Thank you!
2005 Jul 15
4
Can't get sample function from "An Introduction to R" to work
I'm trying to figure out R, a piece at a time, hours at a time... I
was trying to copy the sample function in, "An Introduction to R" (for
version 2.1.0) by W. N. Venables, D. M. Smith, page 42. Section 10.1
"Simple examples" provides a sample function which I tried to duplicate
(I'm using Mac OS X 10.3.9, and "R for Mac OS X Aqua GUI v1.11). The
following
2004 Aug 23
1
Two factor ANOVA with lm()
The following is a data frame
> "jjd" <- structure(list(Observations = c(6.8, 6.6, 5.3, 6.1,
7.5, 7.4, 7.2, 6.5, 7.8, 9.1, 8.8, 9.1), LevelA = structure(c(1,
1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3), .Label = c("A1", "A2",
"A3"), class = "factor"), LevelB = structure(c(1, 1, 2, 2,
1, 1, 2, 2, 1, 1, 2, 2), .Label =
2011 Jan 09
1
question about the chow test of poolability
Good day R-listers,
My question is more a statistical question than an R related question,
so please bear with me
i'm currently applying the chow test of poolability
in fact i'm working with panel N=17 T=5 , and my model looks like
this : Yit= a0+B1X1+B2X2+B3X3+B4X4+eit
My question is the following when i'm Testing for the equality of
the coefficients of the unpooled data (the
2010 May 18
0
CFP: 2nd International Conference on Cloud Computing (CloudComp2010)
[Apologies if you receive multiple copies of this announcement]
Reminder:
-------
* Paper submission deadline: May 31, 2010
========================================================================
Call for Papers
The 2nd International Conference on Cloud Computing
(CloudComp2010)
Barcelona, Spain, October
2010 May 18
0
CFP: 2nd International Conference on Cloud Computing (CloudComp2010)
[Apologies if you receive multiple copies of this announcement]
Reminder:
-------
* Paper submission deadline: May 31, 2010
========================================================================
Call for Papers
The 2nd International Conference on Cloud Computing
(CloudComp2010)
Barcelona, Spain, October
2010 Jun 09
0
CloudComp2010 Deadline Extended to June 14 (2nd International Conference on Cloud Computing)
Update:
-------
* Paper submission deadline extended to: June 14, 2010
========================================================================
Call for Papers
The 2nd International Conference on Cloud Computing
(CloudComp2010)
Barcelona, Spain, October 25 - 28, 2010
2010 Jun 09
0
CloudComp2010 Deadline Extended to June 14 (2nd International Conference on Cloud Computing)
Update:
-------
* Paper submission deadline extended to: June 14, 2010
========================================================================
Call for Papers
The 2nd International Conference on Cloud Computing
(CloudComp2010)
Barcelona, Spain, October 25 - 28, 2010
2005 Sep 01
3
xen2.0 stable periodic machine freeze
Hello,
I have a Dell PowerEdge 2850 that periodically (6 times in the last 24
hours) freezes, requiring a power cycle in order to come back.
The machine is running a Xen 2.0 stable source install. Here is the
appropriate Grub entry:
title Xen 2.0 / XenLinux 2.6
kernel /boot/xen-2.0.gz dom0_mem=131072
module /boot/vmlinuz-2.6-xen0 root=/dev/sda1 ro console=tty0
Dom0 boots happily, and new
2001 Feb 26
1
glm: family=exponential() or family=Gamma(alpha=1)?
In glm, is there a way to specify the family
to be exponential, or Gamma but with alpha=1?
This is required by Dobson Exercise 4.3d, page 48.
But without specifying alpha=1, the asnwer is the same
as in the book.
If reply to the list, please cc me a copy.
Thanks.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read
2007 Nov 26
1
pass lm( ) a char vector as the variables to be included
Dear Everyone in list:
I am writing some codes to automate the process of
fitting linear models where the names of variables of
models are produced and stored in character vectors.
But I have problems to pass the vectors to the lm( )
because I don't know how to strip the quotation marks
automatically.
Here are the codes of the example of lm( ):
## Annette Dobson (1990) "An Introduction
2004 Mar 23
1
influence.measures, cooks.distance, and glm
Dear list,
I've noticed that influence.measures and cooks.distance gives different
results for non-gaussian GLMs. For example, using R-1.9.0 alpha
(2003-03-17) under Windows:
> ## Dobson (1990) Page 93: Randomized Controlled Trial :
> counts <- c(18,17,15,20,10,20,25,13,12)
> outcome <- gl(3,1,9)
> treatment <- gl(3,3)
> glm.D93 <- glm(counts ~ outcome +
2011 Feb 08
1
Error in example Glm rms package
Hi all!
I've got this error while running
example(Glm)
library("rms")
> example(Glm)
Glm> ## Dobson (1990) Page 93: Randomized Controlled Trial :
Glm> counts <- c(18,17,15,20,10,20,25,13,12)
Glm> outcome <- gl(3,1,9)
Glm> treatment <- gl(3,3)
Glm> f <- glm(counts ~ outcome + treatment, family=poisson())
Glm> f
Call: glm(formula = counts ~
1998 Aug 28
1
R-beta: repeated measures GLM binomial data?
I don't know much about GLM in general or glm in R. Can anyone tell me
how to do the following (in R or some other stat system) or refer me to a
textbook discussion?
Two factor ANOVA repeated measures design.
Each subject gives a percent correct.
I am assuming the correct way to proceed is to fit a generalized linear
model with binomial responses and logistic link. But I have not found a
2008 Jun 08
1
exponential distribution
Dear all,
I've tried to solve the Es. 12, cap 4 of "Introduction to GLM" by Annette Dobson.
It's about the relationship between survival time of leukemia patients and blood cell count.
I tried to fit a model with exponential distribution, first by glm (family gamma and then dispersion parameter fixed to 1) and
then with survreg.
They gave me the same point estimates but the
2013 May 29
1
quick question about glm() example
I don't have a copy of Dobson (1990) from which the glm.D93 example is
taken in example("glm"), but I'm strongly suspecting that these are
made-up data rather than real data; the means of the responses within
each treatment are _identical_ (equal to 16 2/3), so two of the
parameters are estimated as being zero (within machine tolerance). (At
this moment I don't understand