[R-sig-teaching] R Equivalent of SPSS "Split File" and "AutoRecode"

Fri Nov 14 00:31:55 CET 2008

 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,]

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)
> 
> _______________________________________________
> R-sig-teaching at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-teaching
>