similar to: Newbie: how to calculate group averagege?

Hi Sherri, There are examples of topographic maps which you have been pointed to, however, I suspect that you want to know where you can obtain topographic data from rather than a canned example. There are quite a few intricacies to the process so I will go through them for you. (1) Topography files can be found in the geomapdata library. You will probably want to use the maps package too (if

I am simulating an ANOVA model and get a strange behavior from the summary function. To be more specific: please run the following code and see for yourself: the summary()[["r.squared"]] values of two identical models are quite different!! ## 3 x 3 ANOVA of two factors x and z on outcome y s.size <- 300 # the sample size p.z <- c(0.25, 0.5, 0.25) # the probabilities of factor z ##

I'm having trouble creating a looping variable and i can't see wher ethe problem arises from any hep gratfully appreciated First create a table x<-table(SURVEY$n_0,exposed) > x exposed False True Under 16 24 1 16-19 68 9 20-24 190 37 25-34 555 204 35-44 330 87 45-54 198 65 55-64 67 35 65+

Hi, I want to pick 3 "integer" random numbers from 1 to 10. How to do this? Thanks, Yu-Ling Wu __________________________________________________ Do You Yahoo!? Get email at your own domain with Yahoo! Mail. http://personal.mail.yahoo.com/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello all, to obtain a list of all objects in all search paths, I've found the following to work: > biglist <- sapply(1:length(search()), objects) This more obvious one, however, does not work: > biglist <- sapply(search(), objects) Error in pos.to.env(pos) : invalid argument Still, search() gives [1] ".GlobalEnv" "package:ctest" "Autoloads"

How can I access a function created in one data directory in another? I have been sourcing the text file that contains the function each time I go to another data dir. This is cumbersome, and I suspect there is a way to make functions global. Thanks in advance, S. David White Dept of Speech and Hearing Science Ohio State University Columbus, Ohio

Hello, first, I'd like to congratulate the core team to the new R version 1.1. I think it's a great update, with glimpses into an even greater future (tcltk!). While playing around with the new functions (on Win 95), I found the following: As the documentation states, when I generate a plot and save it with recordPlot, I can regenerate it by printing the variable: > plot(1:10) >

Currently (R-1.4.1 as well as R-devel, according to http://stat.ethz.ch/R-alpha/R-devel/library/base/html/par.html), the html version of help(par) shows [...] lab A numerical vector of the form c(x, y, len) which modifies the way that axes are annotated. The values of x and y give the (approximate) number of tickmarks on the x and y axes and len specifies the label size. The default

Hi, when I apply a glm() model in two ways, first with the response in a two column matrix specification with successes and failures y <- matrix(c( 5, 1, 3, 3, 2, 2, 0, 4), ncol=2, byrow=TRUE) X <- data.frame(x1 = factor(c(1,1,0,0)), x2 = factor(c(0,1,0,1))) glm(y ~ x1 + x2, data = X, family="binomial") second with a model matrix that

This is only partly an R question, but this seems a good place to ask. Like many, if not most, statistics programs, R has no post hoc tests for Kruskal Wallis or Friedman, and as far as I can make out they are not available in any of the R packages. I would be grateful if someone could point me in the right direction if I am wrong. As I said this is not uncommon. Parametric Anova procedures

Hi all, I have a problem which I am biting my teeth into unsuccessfully: (x <- data.frame(S1=c("a","b","d","F"),N1=c(2,4,6,9),N2=c(6,NA,0,6))) S1 N1 N2 1 a 2 6 2 b 4 NA 3 d 6 0 4 F 9 6 > is.na(x) S1 N1 N2 1 FALSE FALSE FALSE 2 FALSE FALSE TRUE 3 FALSE FALSE FALSE 4 FALSE FALSE FALSE No I want to be able to do: 1) Setting

I have a list with four fields and 3500 records, from a MySQL table. How can i have the total for each level using table, ftable, xtabs? ex: print(ftable(table[,"field1"]~table[,"field2"])) thanks Danilo Maurizio Veneto Lavoro - Area Osservatorio - -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello, does anyone know of a quick way to create a data frame "like" another, but with more rows? What I'd like to do is this: if mydata is a data.frame like a b c 1 TRUE yes 2 FALSE no 3 TRUE yes I'd like to get mydata2 with the same column names and column types, but without the values and with more rows. All I could think of was to manually do

I have a function in S, and I need to write the same function in R. Exist a program to translate from S to R ? Thank you Laura Pontiggia -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the

Hello, I'd like to propose an extension to the function summary.aov. In Splus (2000, I don't know about other versions), summary.aov allows a parameter ssType to be set to 1 or 3 (defaults to 1) to choose the type of Sums of Squares. I know I can get Type III SS in R with drop1(model), but including the functionality into summary.aov would, in my opinion, - yield a more usable table

On Thu, Jun 22, 2000 at 05:01:32PM +0200, p.dalgaard@biostat.ku.dk wrote: > I think this is pretty certainly a bug, so I'm cc'ing this to r-bugs. > > -p > > Kaspar Pflugshaupt <pflugshaupt@geobot.umnw.ethz.ch> writes: > > As the documentation states, when I generate a plot and save it with > > recordPlot, I can regenerate it by printing the variable:

On Thu, Jun 22, 2000 at 05:01:32PM +0200, p.dalgaard@biostat.ku.dk wrote: > I think this is pretty certainly a bug, so I'm cc'ing this to r-bugs. > > -p > > Kaspar Pflugshaupt <pflugshaupt@geobot.umnw.ethz.ch> writes: > > As the documentation states, when I generate a plot and save it with > > recordPlot, I can regenerate it by printing the variable:

Hi R experts, Can somebody suggest what can be wrong? PIII with 128 Mb RAM Using R 1.1 in Win98, the latest version downloaded yesterday. I run update.packages so 'survival', ' mass', and 'nlme' are the latest versions in CRAN > library(survival) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library

Hello all! I'd like to use a routine in R to find the zero of a function (like Newton-Raphson, Broyden, etc.) In the manual, I found a reference to Newton-Raphson method in the "qr" function, but I simply wonder how I could use it to implement such a method. Carlos -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello: I should first apologize for the very basic question, but I ought to start from somewhere [embarrassed grin]. I am using R-1.0.0 in a Windows 98 machine to plot some graphs. I tried to apply a "zoom in" in some portions of the graph using say 'par(yaxp=c(0.35,0.40,5))' but nothing happens, no error msg, no update in the active graphics device. What should be the