On Thu, Nov 13, 2008 at 5:31 PM, Albyn Jones <[EMAIL PROTECTED]> wrote: > I have no clue what "autorecode" means. Selecting a subset of a file is > easy. Suppose your data frame is named X, assuming REGION is numeric: > > attach(X) > X17 = X[REGION == 17,]
I think it is better style to use X17 <- subset(X, REGION == 17) Using attach/detach is discouraged these days. The "subset" function, the"with" function and the general use of the formula/data arguments in model fitting functions and the lattice graphics functions provide cleaner approaches than using attach. The problem with attach is that it has the side-effect of changing the search path, which is, of course, what it is supposed to do. However, any functions with side-effects should be used only when you can't find a better alternative and in this case there are better alternatives. > > albyn > > On Thu, Nov 13, 2008 at 05:18:24PM -0600, Richardson, Paul wrote: >> >> Greetings Fellow Teachers of R, >> >> >> * >> I am working on a tutorial for a class to be taught at my college. >> * >> First, I am going to provide a PowerPoint seminar for the faculty, and then >> I will create a series of Camtasia screen recordings with narration. >> * >> The course is aimed at teaching those familiar with introductory level stats >> how to perform similar actions in R, that would normally be performed in >> SPSS. >> * >> Is there an autorecode capability in R which is analogous to autorecode in >> SPSS? I'm having trouble with item with item f below (second half) >> "autorecode". >> * >> I am also having trouble with item "i" below, because I don't know how to >> split a data file, like you normally do in SPSS. >> * >> Is there some equivalent way to do this in R? >> >> thanks much, >> Paul Richardson >> cell 806-543-0087 >> >> >> a. Go to AEIS data site: >> http://www.tea.state.tx.us/perfreport/aeis/2007/DownloadData.html >> b. Select Campus, Reference File, and Layout, and download; then do the >> same for the Reference data file. Download Campus Financial Statistics, >> both layout and download. Download Student Statistics, both layout and >> download. >> c. Import the data files into R. >> d. Using the layout file, create variable labels; create value labels >> (if possible). Assign 0 as a missing value for variable BPETALLC. >> e. Merge the Reference, Financial, and Student files >> f. Recode variable C_RATING into a new, numeric variable named >> campus.rating in R. Recode GRDTYPE in a new, numeric variable called >> grade.type in R. Use autorecode (if there is something similar in R) to >> recode CNTYNAME into county.name. Create a new variable called >> log.cperra1r, which is the logarithm of CPERRA1R. >> g. Do summary statistics (minimum, maximum, range, mean, median, >> variance, standard deviation, skewness, kurtosis). >> h. Run histograms, boxplots, scatterplots, scatterplot matrices >> i. Select only cases from Region 17 (using REGION variable) and do >> the analyses for j-r below; using the equivalent of split file (if possible) >> to do same analysis for each region separately (using REGION variable). >> j. Conduct 1-way ANOVA using campus.rating as IV and CPETECOP as DV >> k. Conduct repeated measures ANOVA using BPFEAOPRT and CPFEAOPRT >> (although these are measuring different things, I am having you use them for >> purely instructional purposes) >> l. Conduct factorial ANOVA using campus.rating and grade.type as IV's >> and CPETECOP as DV >> m. Conduct correlation using CPFPGSPET and CPFPGGIFT. >> n. Conduct bivariate regression using CPFPGGIFK and CPFPGSPEK >> o. Conduct multiple regression using CPFPGREGT as criterion/outcome >> variable and CPFPAVOCT, CPFPASPET, and CPFPABILT as predictor variables >> p. Conduct the following regression diagnostics >> i. Numerical diagnostics >> 1. Variance inflation factors >> 2. Breusch-Pagan test for heteroscedasticity >> 3. Durbin-Watson test for autocorrelation >> 4. RESET test for nonlinearity >> 5. Bonferonni outlier test >> ii. Graphs >> 1. Basic diagnostic plots >> 2. Residual quantile comparison plot >> 3. Component+residuals plot >> 4. Added variables plot >> 5. Influence plot >> 6. Effects plot >> q. Conduct chi-square using campus.rating and grade.type. >> r. Aggregate all the variables in the file by region (using REGION) >> and save to a new file called cref_cfin_cstud_agg.rdat >> >> >> >> >> Below is the actual syntax for the above tasks: >> ======================// beginning of R syntax //======================= >> # Import CREF >> CREF <- sqlQuery(channel = 1, select * from [cref_tabdelimited$]) >> attach(CREF) >> # SORT >> CREF<- order(+CREF$campus) >> CREF >> # variable labels >> library(Hmisc) >> label(CREF$DISTNAME)<- "District Name" >> label(CREF$CAMPNAME)<- "Campus Name" >> label(CREF$CACKDTL)<- "Campus 2007 Ratings: GPA Acknowledgment Detail" >> label(CREF$CAMPUS)<- "CURRENT 9-DIGIT CDC NUMBER" >> label(CREF$DISTRICT)<- "6 Digit County District Number" >> label(CREF$COUNTY)<- "3 Digit County Number" >> label(CREF$REGION)<- "2 Digit ESC Region" >> label(CREF$CNTYNAME)<- "County Name" >> label(CREF$GRDTYPE)<- "Campus Type (B,E,M,S) Total Enrollment" >> label(CREF$GRDSPAN)<- "Grade Span (Total Enrollment)" >> label(CREF$CFLCHART)<- "CAMPUS 2007 Flag-Chater School (Y/N)" >> label(CREF$PAIRCAMP)<- "CAMPUS 2007 Flag-Campus/District No. of TAKS Giver" >> label(CREF$PAIRNAME)<- "CAMPUS 2007 Flag-Campus/District Name of TAKS Giver" >> label(CREF$C_RATING)<- "Campus 2007 Ratings: Final Accountability Rating" >> #Import CSTUD >> CSTUD <- read.table("C:/Users/Matt/Desktop/TTU/Courses/EPSY 5385_Stats >> Packages/DATA FILES/cstud.dat", header=TRUE, sep=",", na.strings="NA", >> dec=".", strip.white=TRUE) >> attach(CSTUD) >> # SORT >> CSTUD<- order(+CSTUD$campus) >> CSTUD >> # variable labels >> library(Hmisc) >> label(CSTUD$CA0GR06N)<- "Campus 2006 Graduates: All Students Count" >> label(CSTUD$CA0GH06N)<- "Campus 2006 Graduates: Recom HS Pgm All Students >> Count" >> label(CSTUD$CA0GM06N)<- "Campus 2006 Graduates: Min HS Pgm All Students >> Count" >> label(CSTUD$CB0GR06N)<- "Campus 2006 Graduates: African American Count" >> label(CSTUD$CH0GR06N)<- "Campus 2006 Graduates: Hispanic Count" >> label(CSTUD$CW0GR06N)<- "Campus 2006 Graduates: White Count" >> label(CSTUD$CI0GR06N)<- "Campus 2006 Graduates: Native American Count" >> label(CSTUD$CP0GR06N)<- "Campus 2006 Graduates: Asian/Pacific Islander Count" >> label(CSTUD$CS0GR06N)<- "Campus 2006 Graduates: Special Ed Count" >> label(CSTUD$BPETALLC)<- "Group 2007 Student: All Students Count" >> label(CSTUD$BA0GR06N)<- "Group 2006 Graduates: All Students Count" >> label(CSTUD$BB0GR06N)<- "Group 2006 Graduates: African American Count" >> label(CSTUD$BH0GR06N)<- "Group 2006 Graduates: Hispanic Count" >> label(CSTUD$BW0GR06N)<- "Group 2006 Graduates: White Count" >> label(CSTUD$BI0GR06N)<- "Group 2006 Graduates: Native American Count" >> label(CSTUD$BP0GR06N)<- "Group 2006 Graduates: Asian/Pacific Islander Count" >> label(CSTUD$BS0GR06N)<- "Group 2006 Graduates: Special Ed Count" >> label(CSTUD$BA0GH06N)<- "Group 2006 Graduates: Recom HS Pgm All Students >> Count" >> label(CSTUD$BA0GM06N)<- "Group 2006 Graduates: Min HS Pgm All Students Count" >> label(CSTUD$BPETDISC)<- "Group 2006 DAEP: Count" >> label(CSTUD$CPERRAKR)<- "Campus 2006 Retention: Regular Ed KG Rate" >> label(CSTUD$CPERRA1R)<- "Campus 2006 Retention: Regular Ed 01 Rate" >> label(CSTUD$CPERRA2R)<- "Campus 2006 Retention: Regular Ed 02 Rate" >> label(CSTUD$CPERRA3R)<- "Campus 2006 Retention: Regular Ed 03 Rate" >> label(CSTUD$CPERRA4R)<- "Campus 2006 Retention: Regular Ed 04 Rate" >> label(CSTUD$CPERRA5R)<- "Campus 2006 Retention: Regular Ed 05 Rate" >> label(CSTUD$CPERRA6R)<- "Campus 2006 Retention: Regular Ed 06 Rate" >> label(CSTUD$CPERRA7R)<- "Campus 2006 Retention: Regular Ed 07 Rate" >> label(CSTUD$CPERRA8R)<- "Campus 2006 Retention: Regular Ed 08 Rate" >> label(CSTUD$CPERSAKR)<- "Campus 2006 Retention: Special Ed KG Rate" >> label(CSTUD$CPERSA1R)<- "Campus 2006 Retention: Special Ed 01 Rate" >> label(CSTUD$CPERSA2R)<- "Campus 2006 Retention: Special Ed 02 Rate" >> label(CSTUD$CPERSA3R)<- "Campus 2006 Retention: Special Ed 03 Rate" >> label(CSTUD$CPERSA4R)<- "Campus 2006 Retention: Special Ed 04 Rate" >> label(CSTUD$CPERSA5R)<- "Campus 2006 Retention: Special Ed 05 Rate" >> label(CSTUD$CPERSA6R)<- "Campus 2006 Retention: Special Ed 06 Rate" >> label(CSTUD$CPERSA7R)<- "Campus 2006 Retention: Special Ed 07 Rate" >> label(CSTUD$CPERSA8R)<- "Campus 2006 Retention: Special Ed 08 Rate" >> label(CSTUD$CPETDISC)<- "Campus 2006 DAEP: Count" >> label(CSTUD$CPETDISP)<- "Campus 2006 DAEP: Percent" >> label(CSTUD$CPETLEPP)<- "Campus 2007 Student: LEP Percent" >> label(CSTUD$CPETECOP)<- "Campus 2007 Student: ECO Percent" >> label(CSTUD$CPETECOC)<- "Campus 2007 Student: ECO Count" >> label(CSTUD$CPETLEPC)<- "Campus 2007 Student: LEP Count" >> label(CSTUD$CPETSPEC)<- "Campus 2007 Student: SPE Count" >> label(CSTUD$CPETGIFC)<- "Campus 2007 Student: GIF/TAL Count" >> label(CSTUD$CPETBILC)<- "Campus 2007 Student: BIL Count" >> label(CSTUD$CPETALLC)<- "Campus 2007 Student: All Students Count" >> label(CSTUD$CPETVOCC)<- "Campus 2007 Student: VOC Count" >> label(CSTUD$CPETINDC)<- "Campus 2007 Student: Native American Count" >> label(CSTUD$CPETPACC)<- "Campus 2007 Student: Asian/Pacific Islander Count" >> label(CSTUD$CPETBLAC)<- "Campus 2007 Student: African American Count" >> label(CSTUD$CPETHISC)<- "Campus 2007 Student: Hispanic Count" >> label(CSTUD$CPETWHIC)<- "Campus 2007 Student: White Count" >> label(CSTUD$CPETGEEC)<- "Campus 2007 Student: EE Count" >> label(CSTUD$CPETGPKC)<- "Campus 2007 Student: PK Count" >> label(CSTUD$CPETGKNC)<- "Campus 2007 Student: KG Count" >> label(CSTUD$CPETG01C)<- "Campus 2007 Student: 01 Count" >> label(CSTUD$CPETG02C)<- "Campus 2007 Student: 02 Count" >> label(CSTUD$CPETG03C)<- "Campus 2007 Student: 03 Count" >> label(CSTUD$CPETG04C)<- "Campus 2007 Student: 04 Count" >> label(CSTUD$CPETG05C)<- "Campus 2007 Student: 05 Count" >> label(CSTUD$CPETG06C)<- "Campus 2007 Student: 06 Count" >> label(CSTUD$CPETG07C)<- "Campus 2007 Student: 07 Count" >> label(CSTUD$CPETG08C)<- "Campus 2007 Student: 08 Count" >> label(CSTUD$CPETG09C)<- "Campus 2007 Student: 09 Count" >> label(CSTUD$CPETG10C)<- "Campus 2007 Student: 10 Count" >> label(CSTUD$CPETG11C)<- "Campus 2007 Student: 11 Count" >> label(CSTUD$CPETG12C)<- "Campus 2007 Student: 12 Count" >> label(CSTUD$CPEMALLC)<- "Campus 2007 Student: 06 Mobility Numerator (Campus >> Only)" >> label(CSTUD$CPETWHIP)<- "Campus 2007 Student: White Percent" >> label(CSTUD$CPETBLAP)<- "Campus 2007 Student: African American Percent" >> label(CSTUD$CPETHISP)<- "Campus 2007 Student: Hispanic Percent" >> label(CSTUD$CPETINDP)<- "Campus 2007 Student: Native American Percent" >> label(CSTUD$CPETPACP)<- "Campus 2007 Student: Asian/Pacific Islander Percent" >> label(CSTUD$CPETSPEP)<- "Campus 2007 Student: SPE Percent" >> label(CSTUD$CPETVOCP)<- "Campus 2007 Student: VOC Percent" >> label(CSTUD$CPETBILP)<- "Campus 2007 Student: BIL Percent" >> label(CSTUD$CPETGIFP)<- "Campus 2007 Student: GIF/TAL Percent" >> label(CSTUD$CPETGEEP)<- "Campus 2007 Student: EE Percent" >> label(CSTUD$CPETGPKP)<- "Campus 2007 Student: PK Percent" >> label(CSTUD$CPETGKNP)<- "Campus 2007 Student: KG Percent" >> label(CSTUD$CPETG01P)<- "Campus 2007 Student: 01 Percent" >> label(CSTUD$CPETG02P)<- "Campus 2007 Student: 02 Percent" >> label(CSTUD$CPETG03P)<- "Campus 2007 Student: 03 Percent" >> label(CSTUD$CPETG04P)<- "Campus 2007 Student: 04 Percent" >> label(CSTUD$CPETG05P)<- "Campus 2007 Student: 05 Percent" >> label(CSTUD$CPETG06P)<- "Campus 2007 Student: 06 Percent" >> label(CSTUD$CPETG07P)<- "Campus 2007 Student: 07 Percent" >> label(CSTUD$CPETG08P)<- "Campus 2007 Student: 08 Percent" >> label(CSTUD$CPETG09P)<- "Campus 2007 Student: 09 Percent" >> label(CSTUD$CPETG10P)<- "Campus 2007 Student: 10 Percent" >> label(CSTUD$CPETG11P)<- "Campus 2007 Student: 11 Percent" >> label(CSTUD$CPETG12P)<- "Campus 2007 Student: 12 Percent" >> label(CSTUD$CPEMALLP)<- "Campus 2007 Student: 06 Mobility Percent (Campus >> Only)" >> label(CSTUD$BPETWHIP)<- "Group 2007 Student: White Percent" >> label(CSTUD$BPETBLAP)<- "Group 2007 Student: African American Percent" >> label(CSTUD$BPETHISP)<- "Group 2007 Student: Hispanic Percent" >> label(CSTUD$BPETINDP)<- "Group 2007 Student: Native American Percent" >> label(CSTUD$BPETPACP)<- "Group 2007 Student: Asian/Pacific Islander Percent" >> label(CSTUD$BPETECOP)<- "Group 2007 Student: ECO Percent" >> label(CSTUD$BPETLEPP)<- "Group 2007 Student: LEP Percent" >> label(CSTUD$BPETSPEP)<- "Group 2007 Student: SPE Percent" >> label(CSTUD$BPETVOCP)<- "Group 2007 Student: VOC Percent" >> label(CSTUD$BPETBILP)<- "Group 2007 Student: BIL Percent" >> label(CSTUD$BPETGIFP)<- "Group 2007 Student: GIF/TAL Percent" >> label(CSTUD$BPETGEEP)<- "Group 2007 Student: EE Percent" >> label(CSTUD$BPETGPKP)<- "Group 2007 Student: PK Percent" >> label(CSTUD$BPETGKNP)<- "Group 2007 Student: KG Percent" >> label(CSTUD$BPETG01P)<- "Group 2007 Student: 01 Percent" >> label(CSTUD$BPETG02P)<- "Group 2007 Student: 02 Percent" >> label(CSTUD$BPETG03P)<- "Group 2007 Student: 03 Percent" >> label(CSTUD$BPETG04P)<- "Group 2007 Student: 04 Percent" >> label(CSTUD$BPETG05P)<- "Group 2007 Student: 05 Percent" >> label(CSTUD$BPETG06P)<- "Group 2007 Student: 06 Percent" >> label(CSTUD$BPETG07P)<- "Group 2007 Student: 07 Percent" >> label(CSTUD$BPETG08P)<- "Group 2007 Student: 08 Percent" >> label(CSTUD$BPETG09P)<- "Group 2007 Student: 09 Percent" >> label(CSTUD$BPETG10P)<- "Group 2007 Student: 10 Percent" >> label(CSTUD$BPETG11P)<- "Group 2007 Student: 11 Percent" >> label(CSTUD$BPETG12P)<- "Group 2007 Student: 12 Percent" >> label(CSTUD$BPEMALLP)<- "Group 2007 Student: 06 Mobility Percent (Campus >> ONLY)" >> label(CSTUD$BPERRAKR)<- "Group 2006 Retention: Regular Ed KG Rate" >> label(CSTUD$BPERRA1R)<- "Group 2006 Retention: Regular Ed 01 Rate" >> label(CSTUD$BPERRA2R)<- "Group 2006 Retention: Regular Ed 02 Rate" >> label(CSTUD$BPERRA3R)<- "Group 2006 Retention: Regular Ed 03 Rate" >> label(CSTUD$BPERRA4R)<- "Group 2006 Retention: Regular Ed 04 Rate" >> label(CSTUD$BPERRA5R)<- "Group 2006 Retention: Regular Ed 05 Rate" >> label(CSTUD$BPERRA6R)<- "Group 2006 Retention: Regular Ed 06 Rate" >> label(CSTUD$BPERRA7R)<- "Group 2006 Retention: Regular Ed 07 Rate" >> label(CSTUD$BPERRA8R)<- "Group 2006 Retention: Regular Ed 08 Rate" >> label(CSTUD$BPERSAKR)<- "Group 2006 Retention: Special Ed KG Rate" >> label(CSTUD$BPERSA1R)<- "Group 2006 Retention: Special Ed 01 Rate" >> label(CSTUD$BPERSA2R)<- "Group 2006 Retention: Special Ed 02 Rate" >> label(CSTUD$BPERSA3R)<- "Group 2006 Retention: Special Ed 03 Rate" >> label(CSTUD$BPERSA4R)<- "Group 2006 Retention: Special Ed 04 Rate" >> label(CSTUD$BPERSA5R)<- "Group 2006 Retention: Special Ed 05 Rate" >> label(CSTUD$BPERSA6R)<- "Group 2006 Retention: Special Ed 06 Rate" >> label(CSTUD$BPERSA7R)<- "Group 2006 Retention: Special Ed 07 Rate" >> label(CSTUD$BPERSA8R)<- "Group 2006 Retention: Special Ed 08 Rate" >> label(CSTUD$BPETDISP)<- "Group 2006 DAEP: Percent" >> CSTUD$BPETALLC [CSTUD$BPETALLC == 0] <- NA >> CSTUD$BPETALLC [is.na(CSTUD$BPETALLC)] <- 0 >> # Import CFIN >> CFIN <- read.table("C:/Users/Matt/Desktop/TTU/Courses/EPSY 5385_Stats >> Packages/DATA FILES/cfin.dat", header=TRUE, sep=",", na.strings="NA", >> dec=".", strip.white=TRUE) >> attach(CFIN) >> # SORT >> CFIN<- order(+CFIN$campus) >> CFIN >> # variable labels >> library(Hmisc) >> label(CFIN$CPFEGADST)<- "Campus 2006 Finance: Expenditure by Function-School >> Leadership, General Fund" >> label(CFIN$CPFEGADIT)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional Leadership, General Fund" >> label(CFIN$CPFEAADST)<- "Campus 2006 Finance: Expenditure by Function-School >> Leadership, All Funds" >> label(CFIN$CPFEAADIT)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional Leadership, All Funds" >> label(CFIN$CPFPGREGT)<- "Campus 2006 Finance: Expenditure by >> Program-Regular, GeneralFund" >> label(CFIN$CPFPGVOCT)<- "Campus 2006 Finance: Expenditure by Program-Career >> & Technology, General Fund" >> label(CFIN$CPFPGSPET)<- "Campus 2006 Finance: Expenditure by Program-Special >> Education, General Fund" >> label(CFIN$CPFPGGIFT)<- "Campus 2006 Finance: Expenditure by Program-Gifted >> & Talented, General Fund" >> label(CFIN$CPFPGBILT)<- "Campus 2006 Finance: Expenditure by >> Program-Bilingual, General Fund" >> label(CFIN$CPFPAREGT)<- "Campus 2006 Finance: Expenditure by >> Program-Regular, All Funds" >> label(CFIN$CPFPAVOCT)<- "Campus 2006 Finance: Expenditure by Program-Career >> & Technology, All Funds" >> label(CFIN$CPFPASPET)<- "Campus 2006 Finance: Expenditure by Program-Special >> Education, All Funds" >> label(CFIN$CPFPAGIFT)<- "Campus 2006 Finance: Expenditure by Program-Gifted >> & Talented, All Funds" >> label(CFIN$CPFPABILT)<- "Campus 2006 Finance: Expenditure by >> Program-Bilingual, All Funds" >> label(CFIN$CPFEGINST)<- "Campus 2006 Finance: Expenditure by >> Function-Instruction, General Fund" >> label(CFIN$CPFEAINST)<- "Campus 2006 Finance: Expenditure by >> Function-Instruction, All Funds" >> label(CFIN$CPFEGOPRT)<- "Campus 2006 Finance: Expenditure by Function-Total >> Operating,General Fund" >> label(CFIN$CPFEAOPRT)<- "Campus 2006 Finance: Expenditure by Function-Total >> Operating, All Funds" >> label(CFIN$CPFPGALLT)<- "Campus 2006 Finance: Expenditure by Program-Total, >> General Fund" >> label(CFIN$CPFPAALLT)<- "Campus 2006 Finance: Expenditure by Program-Total, >> All Funds" >> label(CFIN$CPFEGINSK)<- "Campus 2006 Finance: Expenditure by >> Function-Instruction PerPupil, General Fund" >> label(CFIN$CPFEGADIK)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional Leadership Per Pupil, General Fund" >> label(CFIN$CPFEGADSK)<- "Campus 2006 Finance: Expenditure by Function-School >> Leadership Per Pupil, General Fund" >> label(CFIN$CPFEGOPRK)<- "Campus 2006 Finance: Expenditure by Function-Total >> OperatingPer Pupil, General Fund" >> label(CFIN$CPFPGREGK)<- "Campus 2006 Finance: Expenditure by Program-Regular >> Per Pupil, General Fund" >> label(CFIN$CPFPGVOCK)<- "Campus 2006 Finance: Expenditure by Program-Career >> & Technology Per Pupil, General Fund" >> label(CFIN$CPFPGBILK)<- "Campus 2006 Finance: Expenditure by >> Program-Bilingual Per Pupil, General Fund" >> label(CFIN$CPFPGGIFK)<- "Campus 2006 Finance: Expenditure by Program-Gifted >> & Talented Per Pupil, General Fund" >> label(CFIN$CPFPGSPEK)<- "Campus 2006 Finance: Expenditure by Program-Special >> Education Per Pupil, General Fund" >> label(CFIN$CPFPGALLK)<- "Campus 2006 Finance: Expenditure by Program-Total >> Per Pupil,General Fund" >> label(CFIN$CPFEAINSK)<- "Campus 2006 Finance: Expenditure by >> Function-Instruction PerPupil, All Funds" >> label(CFIN$CPFEAADIK)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional Leadership Per Pupil, All Funds" >> label(CFIN$CPFEAADSK)<- "Campus 2006 Finance: Expenditure by Function-School >> Leadership Per Pupil, All Funds" >> label(CFIN$CPFEAOPRK)<- "Campus 2006 Finance: Expenditure by Function-Total >> OperatingPer Pupil, ALL Funds" >> label(CFIN$CPFPAREGK)<- "Campus 2006 Finance: Expenditure by Program-Regular >> Per Pupil, All Funds" >> label(CFIN$CPFPAVOCK)<- "Campus 2006 Finance: Expenditure by Program-Career >> & Technology Per Pupil, All Funds" >> label(CFIN$CPFPABILK)<- "Campus 2006 Finance: Expenditure by >> Program-Bilingual Per Pupil, All Funds" >> label(CFIN$CPFPAGIFK)<- "Campus 2006 Finance: Expenditure by Program-Gifted >> & Talented Per Pupil, All Funds" >> label(CFIN$CPFPASPEK)<- "Campus 2006 Finance: Expenditure by Program-Special >> Education Per Pupil, All Funds" >> label(CFIN$CPFPAALLK)<- "Campus 2006 Finance: Expenditure by Program-Total >> Per Pupil,All Funds" >> label(CFIN$CPFEGRELT)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional-Related Services, General Fund" >> label(CFIN$CPFEARELT)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional-Related Services, All Funds" >> label(CFIN$CPFEGSUPT)<- "Campus 2006 Finance: Expenditure by >> Function-Support Services-Student, General Fund" >> label(CFIN$CPFEASUPT)<- "Campus 2006 Finance: Expenditure by >> Function-Support Services-Student, All Funds" >> label(CFIN$CPFPGCOMT)<- "Campus 2006 Finance: Expenditure by >> Program-Compensatory, General Fund" >> label(CFIN$CPFPACOMT)<- "Campus 2006 Finance: Expenditure by >> Program-Compensatory, All Funds" >> label(CFIN$CPFPGOTHT)<- "Campus 2006 Finance: Expenditure by Program-Other, >> General Fund" >> label(CFIN$CPFPAOTHT)<- "Campus 2006 Finance: Expenditure by Program-Other, >> All Funds" >> label(CFIN$CPFEGOPRP)<- "Campus 2006 Finance: Expenditure by Function-Total >> OperatingPercent, General Fund" >> label(CFIN$CPFEGINSP)<- "Campus 2006 Finance: Expenditure by >> Function-Instruction Percent, General Fund" >> label(CFIN$CPFEGRELP)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional-Related Services Percent, General Fund" >> label(CFIN$CPFEGADIP)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional Leadership Percent, General Fund" >> label(CFIN$CPFEGADSP)<- "Campus 2006 Finance: Expenditure by Function-School >> Leadership Percent, General Fund" >> label(CFIN$CPFEGSUPP)<- "Campus 2006 Finance: Expenditure by >> Function-Support Services-Student Percent, General Fund" >> label(CFIN$CPFEAOPRP)<- "Campus 2006 Finance: Expenditure by Function-Total >> OperatingPercent, All Funds" >> label(CFIN$CPFEAINSP)<- "Campus 2006 Finance: Expenditure by >> Function-Instruction Percent, All Funds" >> label(CFIN$CPFEARELP)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional-Related Services Percent, All Funds" >> label(CFIN$CPFEAADIP)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional Leadership Percent, All Funds" >> label(CFIN$CPFEAADSP)<- "Campus 2006 Finance: Expenditure by Function-School >> Leadership Percent, All Funds" >> label(CFIN$CPFEASUPP)<- "Campus 2006 Finance: Expenditure by >> Function-Support Services-Student Percent, All Funds" >> label(CFIN$CPFPGALLP)<- "Campus 2006 Finance: Expenditure by Program-Total >> Percent, General Fund" >> label(CFIN$CPFPGBILP)<- "Campus 2006 Finance: Expenditure by >> Program-Bilingual Percent, General Fund" >> label(CFIN$CPFPGVOCP)<- "Campus 2006 Finance: Expenditure by Program-Career >> & Technology Percent, General Fund" >> label(CFIN$CPFPGCOMP)<- "Campus 2006 Finance: Expenditure by >> Program-Compensatory Percent, General Fund" >> label(CFIN$CPFPGGIFP)<- "Campus 2006 Finance: Expenditure by Program-Gifted >> & Talented Percent, General Fund" >> label(CFIN$CPFPGREGP)<- "Campus 2006 Finance: Expenditure by Program-Regular >> Percent,General Fund" >> label(CFIN$CPFPGSPEP)<- "Campus 2006 Finance: Expenditure by Program-Special >> Education Percent, General Fund" >> label(CFIN$CPFPGOTHP)<- "Campus 2006 Finance: Expenditure by Program-Other >> Percent, General Fund" >> label(CFIN$CPFPAALLP)<- "Campus 2006 Finance: Expenditure by Program-Total >> Percent, All Funds" >> label(CFIN$CPFPABILP)<- "Campus 2006 Finance: Expenditure by >> Program-Bilingual Percent, All Funds" >> label(CFIN$CPFPAVOCP)<- "Campus 2006 Finance: Expenditure by Program-Career >> & Technology Percent, All Funds" >> label(CFIN$CPFPACOMP)<- "Campus 2006 Finance: Expenditure by >> Program-Compensatory Percent, All Funds" >> label(CFIN$CPFPAGIFP)<- "Campus 2006 Finance: Expenditure by Program-Gifted >> & Talented Percent, All Funds" >> label(CFIN$CPFPAREGP)<- "Campus 2006 Finance: Expenditure by Program-Regular >> Percent,All Funds" >> label(CFIN$CPFPASPEP)<- "Campus 2006 Finance: Expenditure by Program-Special >> Education Percent, All Funds" >> label(CFIN$CPFPAOTHP)<- "Campus 2006 Finance: Expenditure by Program-Other >> Percent, All Funds" >> label(CFIN$CPFEGRELK)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional-Related Services Per Pupil, General Fund" >> label(CFIN$CPFEARELK)<- "Campus 2006 Finance: Expenditure by >> Function-Instructional-Related Services Per Pupil, All Funds" >> label(CFIN$CPFEGSUPK)<- "Campus 2006 Finance: Expenditure by >> Function-Support Services-Student Per Pupil, General Fund" >> label(CFIN$CPFEASUPK)<- "Campus 2006 Finance: Expenditure by >> Function-Support Services-Student Per Pupil, All Funds" >> label(CFIN$CPFPGCOMK)<- "Campus 2006 Finance: Expenditure by >> Program-Compensatory PerPupil, General Fund" >> label(CFIN$CPFPACOMK)<- "Campus 2006 Finance: Expenditure by >> Program-Compensatory PerPupil, All Funds" >> label(CFIN$CPFPGOTHK)<- "Campus 2006 Finance: Expenditure by Program-Other >> Per Pupil,General Fund" >> label(CFIN$CPFPAOTHK)<- "Campus 2006 Finance: Expenditure by Program-Other >> Per Pupil,All Funds" >> label(CFIN$BPFEAOPRT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Total Operating, All Funds" >> label(CFIN$BPFEAINST)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instruction, All Funds" >> label(CFIN$BPFEARELT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instructional-Related Services, All Funds" >> label(CFIN$BPFEAADIT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instructional Leadership, All Funds" >> label(CFIN$BPFEAADST)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-School Leadership, All Funds" >> label(CFIN$BPFEASUPT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Support Services--Student, All Funds" >> label(CFIN$BPFPAALLT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Total, All Funds" >> label(CFIN$BPFPABILT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Bilingual, All Funds" >> label(CFIN$BPFPAVOCT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Career & Technology, All Funds" >> label(CFIN$BPFPACOMT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Compensatory, All Funds" >> label(CFIN$BPFPAGIFT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Gifted & Talented, All Funds" >> label(CFIN$BPFPAREGT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Regular, All Funds" >> label(CFIN$BPFPASPET)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Special Education, All Funds" >> label(CFIN$BPFPAOTHT)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Other, All Funds" >> label(CFIN$BPFPAALLP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Total Percent, All Funds" >> label(CFIN$BPFPAREGP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Regular Percent, All Funds" >> label(CFIN$BPFPABILP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Bilingual Percent, All Funds" >> label(CFIN$BPFPACOMP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Compensatory Percent, All Funds" >> label(CFIN$BPFPAGIFP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Gifted & Talented Percent, All Funds" >> label(CFIN$BPFPAVOCP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Career & Technology Percent, All Funds" >> label(CFIN$BPFPASPEP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Special Education Percent, All Funds" >> label(CFIN$BPFPAOTHP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Other Percent, All Funds" >> label(CFIN$BPFEAOPRP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Total Operating Percent, All Funds" >> label(CFIN$BPFEAINSP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instruction Percent, All Funds" >> label(CFIN$BPFEARELP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instructional-Related Services Percent, All Funds" >> label(CFIN$BPFEAADIP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instructional Leadership Percent, All Funds" >> label(CFIN$BPFEAADSP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-School Leadership Percent, All Funds" >> label(CFIN$BPFEASUPP)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Support Services--Student Percent, All Funds" >> label(CFIN$BPFEAOPRK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Total Operating Per Pupil, All Funds" >> label(CFIN$BPFEAINSK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instruction Per Pupil, All Funds" >> label(CFIN$BPFEARELK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instructional-Related Services Per Pupil, All Funds" >> label(CFIN$BPFEAADIK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Instructional Leadership Per Pupil, All Funds" >> label(CFIN$BPFEAADSK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-School Leadership Per Pupil, All Funds" >> label(CFIN$BPFEASUPK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Function-Support Services--Student Per Pupil, All Funds" >> label(CFIN$BPFPAALLK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Total Per Pupil, All Funds" >> label(CFIN$BPFPABILK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Bilingual Per Pupil, All Funds" >> label(CFIN$BPFPAVOCK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Career & Technology Per Pupil, All Funds" >> label(CFIN$BPFPACOMK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Compensatory Per Pupil, All Funds" >> label(CFIN$BPFPAGIFK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Gifted & Talented Per Pupil, All Funds" >> label(CFIN$BPFPAREGK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Regular Per Pupil, All Funds" >> label(CFIN$BPFPASPEK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Special Education Per Pupil, All Funds" >> label(CFIN$BPFPAOTHK)<- "Campus 2006 Finance: Group Total-Expenditure by >> Program-Other Per Pupil, All Funds" >> # merge two datasets by CAMPUS >> total <- merge(CREF,CSTUD, by="CAMPUS", all=T) >> attach(total) >> total2 <- merge(total,CFIN, by="CAMPUS", all=T) >> attach(total2) >> # Sort >> total2 <- order(+total2$campus) >> total2 >> #Recodes >> total2$campus.rating <- recode(total2$C_RATING, '"A" = 1; "E" = 2; "R" = 3; >> "X" = 4; "L" = 5; ', as.factor.result=TRUE) >> total2$campus.rating <- >> factor(total2$campus.rating,levels=c("1","2","3","4","5"),labels=c("Acceptable","Exemplary","Recognized","Unacceptable","Unknown")) >> total2$grade.type <- recode(total2$GRDTYPE, '"B" = 1; "E" = 2; "M" = 3; "S" >> = 4; ', as.factor.result=TRUE) >> total2$grade.type <- >> factor(total2$grade.type,levels=c("1","2","3","4"),labels=c("Both","Elementary","Middle","Senior")) >> total2$CPERRA1R <- as.numeric(total2$CPERRA1R) >> total2$log.cperra1r <- with(total2, log10(CPERRA1R)) >> # auto-recode -- county.name >> #Summary statistics >> summary(total2) >> # Histogram >> Hist(total2$BA0GH06N, scale="frequency", breaks="Sturges", col="darkgray") >> #Scatterplot >> scatterplot(BA0GR06N~BA0GM06N, reg.line=FALSE, smooth=FALSE, labels=FALSE, >> boxplots='xy', span=0.5, data=total2) >> #Scatterplot matrices >> scatterplot.matrix(~BA0GH06N+BA0GR06N+BH0GR06N+BP0GR06N+BPETDISC, >> reg.line=FALSE, smooth=FALSE, span=0.5, diagonal = 'density', data=total2) >> # Split File -- CAN'T IN R >> # 1-way ANOVA >> total2$CPETECOP <- as.numeric(total2$CPETECOP) >> .Anova<- lm(CPETECOP ~ campus.rating, data=total2) >> anova(.Anova) >> tapply(total2$CPETECOP, total2$campus.rating, mean, na.rm=TRUE) # means >> tapply(total2$CPETECOP, total2$campus.rating, sd, na.rm=TRUE) # std. >> deviations >> tapply(total2$CPETECOP, total2$campus.rating, function(x) sum(!is.na(x))) # >> counts >> remove(.Anova) >> # Repeated measures ANOVA? >> GLM.1 <- glm(BPFEAOPRT ~ CPFEAOPRT, family=binomial(logit), data=total2) >> summary(GLM.1) >> # Factorial ANOVA >> # Correlation >> library(RODBC) >> total2$CPFPGSPET <- as.numeric(total2$CPFPGSPET) >> total2$CPFPGGIFT <- as.numeric(total2$CPFPGGIFT) >> cor.test(total2$CPFPGSPET , total2$CPFPGGIFT, use="pairwise.complete.obs", >> alternative="two.sided", method="pearson") >> # Bivariate Regression >> total2$CPFPGGIFK <- as.numeric(total2$CPFPGGIFK) >> total2$CPFPGSPEK <- as.numeric(total2$CPFPGSPEK) >> RegModel.1 <- lm(CPFPGGIFK~CPFPGSPEK, data=total2) >> summary(RegModel.1) >> anova(RegModel.1) >> plot(RegModel.1) >> termplot(RegModel.1) >> bptest(CPFPGGIFK~CPFPGSPEK, varformula = ~ fitted.values(RegModel.1), >> studentize=FALSE, data=total2) >> dwtest(CPFPGGIFK~CPFPGSPEK, alternative="greater", data=total2) >> resettest(CPFPGGIFK~CPFPGSPEK, power=2:3, type="regressor", data=total2) >> outlier.test(RegModel.1) >> #Graphs >> oldpar <- par(oma=c(0,0,3,0), mfrow=c(2,2)) >> plot(RegModel.1) >> par(oldpar) >> qq.plot(RegModel.1, simulate=TRUE, labels=FALSE) >> cr.plots(RegModel.1, ask=FALSE) >> av.plots(RegModel.1, ask=FALSE, identify.points=TRUE) >> influencePlot(RegModel.1) >> trellis.device(theme="col.whitebg") >> plot(all.effects(RegModel.1), ask=FALSE) >> # Multiple Regression >> total2$CPFPGREGT <- as.numeric(total2$CPFPGREGT) >> total2$CPFPAVOCT <- as.numeric(total2$CPFPAVOCT) >> total2$CPFPASPET <- as.numeric(total2$CPFPASPET) >> total2$CPFPABILT <- as.numeric(total2$CPFPABILT) >> RegModel.2 <- lm(CPFPGREGT~CPFPAVOCT+CPFPASPET+CPFPABILT, data=total2) >> summary(RegModel.2) >> anova(RegModel.2) >> plot(RegModel.2) >> termplot(RegModel.2) >> vif(RegModel.1) >> bptest(CPFPGREGT~CPFPAVOCT+CPFPASPET+CPFPABILT, varformula = ~ >> fitted.values(RegModel.1), studentize=FALSE, data=total2) >> dwtest(CPFPGREGT~CPFPAVOCT+CPFPASPET+CPFPABILT, alternative="greater", >> data=total2) >> resettest(CPFPGREGT~CPFPAVOCT+CPFPASPET+CPFPABILT, power=2:3, >> type="regressor", data=total2) >> outlier.test(RegModel.1) >> #Graphs >> oldpar <- par(oma=c(0,0,3,0), mfrow=c(2,2)) >> plot(RegModel.2) >> par(oldpar) >> qq.plot(RegModel.2, simulate=TRUE, labels=FALSE) >> cr.plots(RegModel.2, ask=FALSE) >> av.plots(RegModel.2, ask=FALSE, identify.points=TRUE) >> influencePlot(RegModel.2) >> trellis.device(theme="col.whitebg") >> plot(all.effects(RegModel.2), ask=FALSE) >> # Chi-square >> chisq.test(total2$campus.rating, total2$grade.type)) >> # Aggregate >> attach(total2) >> cref_cfin_cstud_agg <-aggregate(total2, by=list(REGION), >> FUN=mean, na.rm=TRUE) >> >> _______________________________________________ >> [email protected] mailing list >> https://stat.ethz.ch/mailman/listinfo/r-sig-teaching >> > > _______________________________________________ > [email protected] mailing list > https://stat.ethz.ch/mailman/listinfo/r-sig-teaching > _______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-teaching
