• If you are citizen of an European Union member nation, you may not use this service unless you are at least 16 years old.

  • Stop wasting time looking for files and revisions. Connect your Gmail, DriveDropbox, and Slack accounts and in less than 2 minutes, Dokkio will automatically organize all your file attachments. Learn more and claim your free account.


Catherine, Dima, Toni Discussion Page

Page history last edited by Tony Leonardi 9 years, 6 months ago

This our re-do of the dependent sample assignment from 2/10 Driving Reaction Time Under Distractive Conditions.docx  


Monday 2/7/2011





number of elements/components 


dimensions of an object 


structure of an object 


class or type of an object 


how an object is stored 


names of components in an object 

c(object, object,...) 

combines objects into a vector 


cbind(object, object, ...) 

combines objects as columns 

rbind(object, object, ...) 

combines objects as rows 


prints the object 


list the first part of the object 


list the last part of the object 


list current objects 


delete an object 

newobject <- edit(object

edit object and save as newobject 


edit in place 



To import data from an EXCEL file  - RODBC package to access the files 

install.packages ("RODBC")


Webscraping - the user extracts info imbedded in a web page using:

readlines ( )         grep ( ) gsub ( ) functions


Sunday 2/6/2011

Example for Linear Least Squares Regression:


When we type varialbename<-lm(what ever values we pass), we should issue summary(variablename) to get Residuals, Coefficients, Residual standard error, ...etc. issuing variablename will only show the coefficients! plot(myvariable) will give multiple plots (should discuss how to read them!!)


runif(n, min=0, max=1): It let us generate random deviates; however, unlike mvrnorm we can't specify the mean and the covariance matrix. We forgot to ask how we can specify min and max limit so we don't give unwanted values!


useful commands:

help.search("a topic")

example("a topic")

data() shows all datasets in the loaded packages

sink(file = "filename", append = FALSE, type = c("output", "message"), split = FALSE) to write the output to the file filename :) "split = TRUE write to both the file and the screen. Not quite sure what type is?

sink(): return the output back to the screen.

source("myfile") will submit the code from myfile to the current R session.

png("myfile.png") to redirect graphic output to myfile. We have to issue the drawing command (plot() for example) and then dev.off(). we can also use jpeg(), bmp() or pdf().

help(package="packagename"): to list the functions and datasets included in packagename






Shoe Wear Graphics



granova.ds(data.frame(shoes),revc=TRUE, main="Shoe Wear",xlab="Sole material B",ylab="Sole material A")


 "More about comparing Two independent samples." Assignment:

To create the dataframe: Check http://msenux.redwoods.edu/math/R/dataframe.php for more information.


musicTypeB=c(18, 23, 21, 20, 20, 29, 20, 16, 20, 26, 21, 25, 17, 18, 19)



The updated function "es.twogps" that includes calculating the CI: 


es.twogps <- function(x, y)


{#generates standardized ES for input  


 # x,y (indep.) sample data


mx <- mean(x)


my <- mean(y)


nx <- length(x)


#nb: nx will generally NOT equal ny


ny <- length(y)


md <- mx - my


vpoold <-((nx-1)*var(x) +(ny-1)* 


var(y))/ (nx + ny - 2)


sdpoold <- sqrt(vpoold)


ES <- md/sdpoold















  1. What is: (in simple linear regression)
    • influence statistics
    • Residuals standard error
    • F-statistic
    • Multiple R-squared
    • Adjusted R-squared

Comments (5)

bob pruzek said

at 6:56 pm on Feb 6, 2011

I will provide an alternative way, maybe two, for calculating a ( 95% ) confidence interval for the population mean difference. The direct method you've used is not recommended (tho' it might, as we shall see, be close
numerically to what I shall provide) More on this Monday. bp

Dima Kassab said

at 7:33 pm on Feb 6, 2011

Thanks Bob!

bob pruzek said

at 1:05 am on Feb 7, 2011

Dima (or ??), You asked: "Questions:
What is: (in simple linear regression)
influence statistics
Residuals standard error
Multiple R-squared
Adjusted R-squared"
We have not covered regression yet, so these questions will have to wait a bit for detailed answers. But I will bring a book or two that should help in the meantime. And I can provide many URLs that will also help
if you want a head start.As for some of the other things on this page, I'm not sure if there are q's you want me to answer, but if any of you (in the whole class) do, just make the q's explicit. BP

bob pruzek said

at 1:18 am on Feb 7, 2011

BTW, to simulate data it is probably best to use rnorm or rt functions, and not runif(n, min=0, max=1):
Real data rarely conform to what runif generates! bp

Dima Kassab said

at 1:36 am on Feb 7, 2011

Indeed, these questions are topics any one of our small group (Catherine, Dima and Tony) encounters and doesn't have a clear idea about! We post them here to keep track about areas we need to learn more about. Eventually if no one of our group has an answer and the questions aren't covered in class we'll make the q's explicit.
Thanks for your concern Bob.

You don't have permission to comment on this page.