similar to: Suggestions for R-intro manual (PR#10701)

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