v <- 23:37 v v <- seq( 24.2, 37.8, by=.2 ) v v <- seq( .001,.999, by=.01 ) v help(dnorm) help(dbeta) dens <- dbeta( v, 2, 7) dens help( plot ) plot( v, dens ) plot( v, dens, type = "l" ) plot( v, dens, type = "p" ) plot( v, dens, type = "b" ) repvect <- rep( 1:3, 4 ) repvect repvect <- rep( 1:3, each=4 ) repvect logvect <- c(F,F,F,T,T) logvect v1 <- c( 1,2,3,NA,5) v1 repvect mean( repvect ) mean( v1 ) help(mean) mean( v1, na.rm=T ) is.na(v1) ls() save.image() save( repvect, file="repvect.Rdata") q() rm(dens) rm(logvect) rm(repvect) rm(v) rm(v1_ ) rm(v1) ls() q() ls() load("repvect.Rdata") ls() repvect cuberepvect <- repvect^3 cuberepvect bothvects <- data.frame( rawdata = repvect, cube = cuberepvect ) bothvects <- data.frame( rawdata = repvest, cube = cuberepvect ) bothvects <- data.frame( rawdata = repvect, cube = cuberepvect ) bothvects for( i in 1:12 ) print(repvect[i]^3) cubes <- rep(0,5000) cubes length(cubes) for( i in 1:length(cubes) ) cubes[i] <- i^3 cubes system.time( for( i in 1:length(cubes) ) cubes[i] <- i^3 ) cubes <- rep(0,50000) system.time( for( i in 1:length(cubes) ) cubes[i] <- i^3 ) system.time( cubes <- (1:50000)^3 ) repvect repchar <- rep( c("g","w","p"), 6) repchar is.vector(repchar) repchar3 <- rep( "gwp", 6) repchar3 repfact <- as.factor( repchar ) ls() is.vector( repfact ) is.factor( repfact ) repfact repvect repvect10 <- 10 * repvect repvect10 repfact10 <- as.factor( repvect10 ) repfact10 mean( repfact10 ) help(read.table) newlist <- list( myscalar = 7, myvector = c("a","b","c","a"), myfactor=repfact10) newlist newlist$myvector q()