similar to: Side by side strip charts

Hi, We were wondering how we could make a stacked frequency diagram such as this one: http://msemac.redwoods.edu/~darnold/math15/liz.pdf We don't necessarily need the shaded "balls", other characters would be fine, such as stacks of x's. David

My friend sent an Excel file: http://msemac.redwoods.edu/~darnold/temp/cyu01_iqscores.xls http://msemac.redwoods.edu/~darnold/temp/cyu01_iqscores.xls I opened it in Excel, saved is as cyu01_iqscores.csv, then imported it into R with: iqscores=read.csv('cyu01_iqscores.csv',header=TRUE) The result was: > head(iqscores) IQ.Scores X 1 145 NA 2 101 NA 3 123 NA 4

Hi, Typical things you read when new to stats are cautions about using a t-statistic when comparing independent samples. You are steered toward a pooled test or welch's approximation of the degrees of freedom in order to make the distribution a t-distribution. However, most texts give no information why you have to do this. So I thought I try a little experiment which is outlined here.

Hi, I am encountering a difficulty I don't understand. Be patient, I'm very new to analysis of variance. If I load this data: example12_7=read.table("http://msemac.redwoods.edu/~darnold/math15/data/chapter12/example12_7.dat",header=TRUE) The run the oneway.test: oneway.test(time~drug,data=example12_7,var.equal=TRUE) I get these results: data: time and drug F = 4.1881, num

Hello, I'm trying to plot an exponential regression line through my data (scatterplot with trend line), but can't find a way to do that. it is to be the best fit possible. -- View this message in context: http://www.nabble.com/plotting-regression-line-tp17837683p17837683.html Sent from the R help mailing list archive at Nabble.com.

All, Consider: BagA <- c(-1000,10,10,10,10,10,10, 10,20,20,20,20,20,20,30, 30,40,40,50,60) BagB <- c(10,20,30,30,40,40,50,50, 50,50,50,50,60,60,60,60, 60,60,60,1000) layout(c(2,1)) barplot(table(BagB)) barplot(table(BagA)) At this point, I'd like to arrange the plots so that the 10-bars are aligned, the 20-bars are aligned, etc. So, I started

All, I am looking at an example in Aliaga's Interactive Statistics. Bag A has the following vouchers. BagA <- c(-1000,10,10,10,10,10,10, 10,20,20,20,20,20,20,30, 30,40,40,50,60) Bag B has the following vouchers. BagB <- c(10,20,30,30,40,40,50,50, 50,50,50,50,60,60,60,60, 60,60,60,1000) Two values are selected (from BagA or BagB) without

Hi, I hope that folks can give me some simple approaches to taking the data set below, which is accumulated in two columns called "long" and "group", then arrange the data is the "long" column into a data frame containing five variables: "Group 1", "Group 2", "Group 3", "Group 4", and "Group 5". I am hoping for a few

Hi, Reading about a "Heads and Tails" game in http://www.dartmouth.edu/~chance/teaching_aids/books_articles/probability_book/amsbook.mac.pdf Introduction to Probability (Example 1.4, pp. 5-8). You toss a coin 40 times. If heads, Peter wins $1, tails, he loses $1. I think I can do that ok with: winnings <- sum(sample(c(-1,1), 40, replace=TRUE)) But I have to do it 10,000 times

Hi, Can someone explain this: > options(digits=20) > 1/3 [1] 0.33333333333333331483 Why the 1483 at the end? Thanks, David. david-arnolds-macbook-pro-2:~ darnold$ R --version R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) -- View this

All, Is this typical of how people will print a summary of results? CoinTosses <- function(n) { x <- sample(c(0,1), n, replace=TRUE) y <- x y[y==0] <- "T" y[y==1] <- "H" numHeads <- sum(x) numTails <- n-sum(x) p <- numHeads/n cat(cat(y,sep=""),"\n") cat("Number of heads: ", numHeads, "\n")

All, I believe I am running the latest version of rshape2 (1.2.1). But this code: library(reshape2) tmp <- melt(smiths, id.vars=1:2, measure.vars=c("age","weight","height"), variable.name="myvars", value.name="myvals" ) names(tmp) Produces this output: > names(tmp) [1] "subject" "time"

Hi, Entering data from a given file is the hardest part of learning R for me. For example, I have this datafile from Moore's text: Live Age Count Parents Age19 324 Another Age19 37 OwnPlace Age19 116 Group Age19 58 Other Age19 5 Parents Age20 378 Another Age20 47 OwnPlace Age20 279 Group Age20 60 Other Age20 2 Parents Age21 337 Another Age21 40 OwnPlace Age21 372 Group Age21 49 Other Age21 3

Hi, I managed to use the attached data set and figure out the following: flies <- read.table("example12_1.dat",header=TRUE,sep="\t") boxplot(long ~ group, data = flies, horizontal = TRUE, col = "red") I'm very new to R and would like some help with the following: 1. Change the order on the y-axis from 1, 2, 3, 4, 5 to 5, 4, 3, 2, 1.

All, How come i=1 in the first case, but i=2 in the second case. The second case seems to work, but the first case does not. David. > findruns <- function(x,k) { + n <- length(x) + runs <- NULL + for (i in 1:(n-k+1)) { + if (all(x[i:(i+k-1)]==1)) runs <- c(runs,i) + browser(i>1) + } + return(runs) + } > x=c(1,1,0,0,1,1,1) > findruns(x,2) Called from:

Hi, I see a lot of folks verify the regression identity SST = SSE + SSR numerically, but I cannot seem to find a proof. I wonder if any folks on this list could guide me to a mathematical proof of this fact. Thanks. David. -- View this message in context: http://r.789695.n4.nabble.com/Regression-Identity-tp4636829.html Sent from the R help mailing list archive at Nabble.com.

All, 1. Is there a command that tells you what libraries are loaded? 2. How do you remove a library that has been loaded? 3. If you have loaded a library that masks previously loaded objects, how can you tell which object is being run? In Matlab, I would run "which f" to determine the location of the file f that will be run. Is there a similar procedure in R? David. -- View this

All, Can someone explain why this does not print the contents of x when I source this file? CoinTosses <- function(n,print=TRUE) { x <- sample(c(0,1), n, replace=TRUE) y <- x y[y==0] <- "T" y[y==1] <- "H" p <- sum(x)/n p } x <- CoinTosses(40) x On the other hand, if I source this file: CoinTosses <- function(n,print=TRUE) { x <-

Hi, Currently I have a density plot generated with this snippet. Is there a way I can add a line curve on top of it? I mean in one figure __BEGIN__ myhist <- hist(x col="blue", main = "Density Plot", xlab = "Exp Level", ) __END__ - Gundala Viswanath Jakarta - Indonesia

All, I have: sales <- c(2300,900,155,102,42,10) names(sales) <- c("Christmas","Valentine's Day", "Mother's Day","Father's Day", "Thanksgiving","New Year's Day") barplot(sales,ylim=c(0,2500)) But it doesn't place all of the name labels on the plot. So I tried: sales <-