search for: arcriswell

Displaying 20 results from an estimated 38 matches for "arcriswell".

2004 Nov 25
1
Error in anova(): objects must inherit from classes
...x-gnu system i586, linux-gnu status major 2 minor 0.0 year 2004 month 10 day 04 language R -- Andrew R. Criswell, Ph.D. Graduate School, Bangkok University mailto:andrew.c at bu.ac.th <http://email.bu.ac.th/src/compose.php?send_to=andrew.c%40bu.ac.th> mailto:andrew at arcriswell.com <http://email.bu.ac.th/src/compose.php?send_to=andrew%40arcriswell.com>
2002 Dec 20
5
Getting graphs into LaTeX
Hello ALL: I ran with success the following commands in R getting a file saved ------------------------------------------------------------------------------------ postscript() postscript('~/data/st202/2003/lecture00/lecture00-graph-01.eps', horizontal = FALSE, height = 6, pointsize = 10) hist(trial.outcome.5, breaks = 5, main = '1000 Replications of 5 Trials of a
2002 Mar 11
2
Help with Python, R and RPY
Hi All: Walter Moreira wrote a small extension module for using the R programming language from within Python. Tim Church's example at http://www.cmat.edu.uy/~walterm/rpy was so compelling, I could not resist installing it on my Linux Mandrake 8.1 box. But I ran into problems. I have installed on Mandrake 8.1 both python and R: R-base-1.4.1-1mdk.i586.rpm
2002 May 21
6
Point labels
Hello, Is there an option I can easily set to put labels next to the points in a simple plot? The closest I got to it was to use the option pch, but this replaces the point symbol with the character specified, and it is not what I want. Thanks -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read
2002 May 07
1
Problem with ties in rank()
Hello All: I have a vector of data, z > z [1] 0.1 0.1 0.1 0.1 0.2 0.2 0.3 0.3 0.3 0.4 0.5 0.5 0.5 0.7 0.7 0.7 0.9 0.9 1.1 [20] 1.1 1.2 1.3 1.4 The first 4 elements have values of 0.1 followed 2 elements with values 0.2. When I invoke rank(z), I expected to get (1+2+3+4)/4 = 2.5 for the first 4 elements in the ranking and (5+6)/2 = 5.5 for elements 5 and 6. But what I do
2002 Jul 08
1
R Libraries for ORDINAL categorical data
Hello All: I know the function loglin() and loglm() from librarary(MASS) performs analysis on nominal categorical data. Are there any libraries, functions or examples available for analysis of ordinal categorical data? I have in mind procedures that can replicate results in Alan Agresti (1984) "Analysis of Ordinal Categorical Data." Thanks, ANDREW
2004 Apr 30
2
Arrow movements in R
Hello: I apologize for asking this question again. I recently upgraded to R-1.9.1 but now the arrow keys do not scroll up and down the command lines in R. Instead I get ^[[A My SuSe 9 platform is missing some program. What is it? Thanks, ANDREW
2004 Nov 24
0
problem with anova and glmmPQL
..._ platform i586-mandrake-linux-gnu arch i586 os linux-gnu system i586, linux-gnu status major 2 minor 0.0 year 2004 month 10 day 04 language R -- Andrew R. Criswell, Ph.D. Graduate School, Bangkok University mailto:andrew.c at bu.ac.th mailto:andrew at arcriswell.com
2004 Nov 25
0
MASS problem -- glmmPQL and anova
...F-value p-value (Intercept) 1 2031 5.70343 0.0170 day 1 2031 213.21673 <.0001 envir 1 2031 12.50388 0.0004 stereotypy 1 2031 0.27256 0.6017 -- Andrew R. Criswell, Ph.D. Graduate School, Bangkok University mailto:andrew.c at bu.ac.th mailto:andrew at arcriswell.com
2001 Oct 29
1
Help with 'get.hist.quote' on tseries
Hi ALL: I am trying to use get.help.quote from library(tseries). I tried to run the example from help(get.hist.quote) but R complained. Here is the command I used and the response: ibm <- get.hist.quote(instrument = "ibm", start = "1998-01-01") trying URL
2002 May 07
1
More on ties with rank()
I am grateful to Prof. Ripley for his explanation. Indeed, rounding explains it all. I take the difference between two vectors and call it "z" method.a <- c(6.3, 6.3, 3.5, 5.1, 5.5, 7.7, 6.3, 2.8, 3.4, 5.7, 5.6, 6.2, 6.6, 7.7, 7.4, 5.6, 6.3, 8.4, 5.6, 4.8, 4.3, 4.2, 3.3, 3.8, 5.7, 4.1) method.b <- c(5.2, 6.6, 2.3, 4.4, 4.1, 6.4, 5.4, 2.3, 3.2, 5.2, 4.9,
2002 Dec 11
2
ordering x's and y's
Hello ALL: How do I get R to list all possible orderings of 2 x's and 3 y's? It should look like this (which rows appear first is unimportant): x x y y y x y x y y x y y x y x y y y x y x x y y y x y x y y x y y x y y x x y y y x y x y y y x x Thanks, ANDREW
2004 Nov 26
1
help with glmmPQL
Hello: Will someone PLEASE help me with this problem. This is the third time I've posted it. When I appply anova() to two equations estimated using glmmPQL, I get a complaint, > anova(fm1, fm2) Error in anova.lme(fm1, fm2) : Objects must inherit from classes "gls", "gnls" "lm","lmList", "lme","nlme","nlsList", or
2002 Sep 22
2
Daglaard's book
Hello all: Have any of you succeeded in getting hold of Peter Daglaard's book, "Introductory Statistics with R"? I have placed an order some weeks ago with AMAZON.CO.UK and it appears I must wait another 6 weeks before delivery. ANDREW -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read
2004 Jul 12
5
help with paste
Hello All: Suppose the following little data frame: > x <- data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1)) > > x$cat [1] 8 2 3 6 1 > How can I get the paste() function to do the same thing. The command below is obviously wrong > paste(x, cat, sep = "$") > Thanks, ANDREW
2003 Oct 02
6
how calculate mean for each group
Hello, R experts: I got data like this: group duplicate treatment A Y 5 A Y 3 A N 6 B Y 2 B N 4 B Y 1 How to sort the data and calculate the average treatment value for each group in two level of duplicate. Results like this: group duplicate treatment A Y 4 A N
2005 Sep 04
1
specification for glmmPQL
Hello All, I have a question regarding how glmmPQL should be specified. Which of these two is correct? summary(fm.3 <- glmmPQL(cbind(response, 100 - response) ~ expt, data = data.1, random = ~ 1 | subject, family = binomial)) summary(fm.4 <- glmmPQL(response ~ expt, data = data.2, random = ~ 1 | subject, family =
2001 Oct 11
2
Where's MVA?
Hi All: Package TSERIES is stated to depend on MVA. However, there is no MVA package to be found under the list of package sources. Best wishes, ANDREW tseries: Package for time series analysis Package for time series analysis with emphasis on non-linear and non-stationary modelling Version: 0.7-6 Depends: ts, mva, quadprog Date: 2001-08-27 Author: Compiled by Adrian
2001 Nov 20
3
Sampling from a population
Hi ALL: Suppose you have a population of N <- 5 observations, x <- c(43, 28, 7, 61, 39). From that you can draw a maximum of 10 samples without replacement of size n <- 3. (Command choose(N,n) yields 10). For instance the samples I seek are 43, 61, 7 39, 7, 28 ...etc How can I get R to do that for me, to get an exhaustive list of samples of size n drawn without replacement
2003 Dec 16
6
Resampling Stats software
Hi, I am new to R (I have most of my experience in SAS and SPSS). I was wondering if anyone has used both Resampling Stats and R, and could comment on strengths/relationships. Also, I have no clue on how to do the various examples from the book "Resampling: The New Statistics" in R. Can anyone give me some possible starting points? Or websites/books? Thanks, Brandon