similar to: plotting regression line

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

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, I'm looking for some ideas on how to reproduce the attached image in R. There are three samples, each of size n = 10. The first is drawn from a normal distribution with mean 60 and standard deviation 3. The second is drawn from a normal distribution with mean 65 and standard deviation 3. The third is drawn from a normal distribution with mean 70 and standard deviation 3.

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

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

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

Hello, I am trying to make a barplot with nested count data which is build like this: first there are several birds (n)laying 3 clutches composed of 2 eggs half of the second and third clutch received treatment and this treatment was tested to influence sex of offspring. I want a barplot showing counts for male and female for every egg of every clutch. can someone tell me what to do? I drew a

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

I made a plot, but after I made a second plot, the previous plot was gone. How can I save all the plots in a file (I do not manually copy and paste them one by one)? Thanks. Linda

Sorry, I have wrote a wrong subject in the first email! Regards, Linda ---------- Forwarded message ---------- From: linda Porz <linda.porz@gmail.com> Date: 2011/5/25 Subject: combined odds ratio To: r-help@r-project.org Cc: r-help-request@stat.math.ethz.ch Dear all, I am looking for an R function which does stepwise selection cox model in r (delta chisq likelihood ratio test) similar

I have two computers ROLAND (RedHat 5.2 and Samba) and LINDA (Windows95). I am already able to print on the printer connacted to LINDA. Yesterday, before I managed to get the printer working, I was alos able to see ROLAND's directory 'tmp' and 'winD' on LINDA. The problem that I saw this morning was that the printing still works but LINDA can no longer see the shares of

Hi there, In the description of command "pam", it mentions "For datasets larger than (say) 200 observations". Now my dataset is a "54732 by 5" dataframe named "test". When I try to run pam(test,4),it shows " cannot allocate vector of length 1497768547". Is it because the row too big that it can't handle? Thank you!

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")

Hi there, I have a question about one command sentence when I follow the example in the book of "Survival analysis in S": > aml1<-aml[aml$group==1] but I got the error warning: NULL data frame with 23 rows Thus, I couldn't keep going on the next command: esf.fit<-survfit(Surv(aml1,status)~1). and also when I try > aml1<-aml[aml$group==1,]

I have my XP-Pro machine setup as a member of an Samba-based Domain server. I usually run as "linda@<localmachine>", but wanted to try running a program using "RunAs" in the "Home", DOMAIN my machine belongs to. The runas command says: RUNAS [ [/noprofile | /profile] [/env] [/netonly] ] /user:<UserName> program with Notes: NOTE: Enter user's

Emmanuel Florac wrote: > id TESTAD\\testuser > returns "no such user" and > > getent passwd TESTAD\\testuser > > returns a "2" code. > ---- On linux, to get 'domain\user' to resolve, I had to have those entries in my /etc/passwd (and /etc/group for groups). I.e. *w/o krb*, (in samba 3.x), I had entries like: