--------------------------------------------------------------------------------------------------------------
log: C:\Documents and Settings\HP_Owner\My Documents\My Classes\SOCY498C\Macros\SOCY498C-HW7-MACROS.tx
> t
log type: text
opened on: 30 Jan 2009, 12:34:29
.
. local problem 2+2
. local solution = 2+2
. di "`problem'"
2+2
. di "`solution'"
4
.
. clear
. capture confirm file "C:\Documents and Settings\HP_Owner\My Documents\My Research\GSS\GSS-Cumulative-72-06.d
> ta"
. if _rc {
. di "_rc"
. }
. else {
. local gss "C:\Documents and Settings\HP_Owner\My Documents\My Research\GSS\GSS-Cumulative-72-06.dta"
. }
.
. #delimit ;
delimiter now ;
. use year
> socrel
> socommun
> socfrend
> socbar
> age
> sex
> race
> educ
> sexfreq if race!=3 using "`gss'", clear;
(General Social Surveys, 1972-2006 [Cumulative File], Dataset 0001)
. #delimit cr
delimiter now cr
.
. #delimit ;
delimiter now ;
. recode sexfreq (0=0)
> (1=1.5)
> (2=12)
> (3=30)
> (4=52)
> (5=130)
> (6=237), gen(sexfreq1)
> ;
(15308 differences between sexfreq and sexfreq1)
. #delimit cr
delimiter now cr
. label variable sexfreq1 "Sexual Frequency (yearly)"
. gen agesqr=age^2
(177 missing values generated)
. label variable agesqr "Age-squared Interaction Term"
.
. sum age
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
age | 48620 45.75072 17.52199 18 89
. return list
scalars:
r(N) = 48620
r(sum_w) = 48620
r(mean) = 45.75071986836693
r(Var) = 307.0200276600631
r(sd) = 17.52198697808165
r(min) = 18
r(max) = 89
r(sum) = 2224400
.
. regress sexfreq1 age agesqr
Source | SS df MS Number of obs = 19679
-------------+------------------------------ F( 2, 19676) = 1699.93
Model | 12626793.4 2 6313396.68 Prob > F = 0.0000
Residual | 73075181.7 19676 3713.92466 R-squared = 0.1473
-------------+------------------------------ Adj R-squared = 0.1472
Total | 85701975.1 19678 4355.21776 Root MSE = 60.942
------------------------------------------------------------------------------
sexfreq1 | Coef. Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
age | -.9939651 .1422126 -6.99 0.000 -1.272714 -.7152163
agesqr | -.0049477 .0014046 -3.52 0.000 -.0077008 -.0021945
_cons | 112.4948 3.292289 34.17 0.000 106.0416 118.9479
------------------------------------------------------------------------------
. ereturn list
scalars:
e(N) = 19679
e(df_m) = 2
e(df_r) = 19676
e(F) = 1699.9258879198
e(r2) = .1473337499679749
e(rmse) = 60.94197785077314
e(mss) = 12626793.36547285
e(rss) = 73075181.69602849
e(r2_a) = .1472470792777907
e(ll) = -108800.9508199172
e(ll_0) = -110369.2399358887
macros:
e(cmdline) : "regress sexfreq1 age agesqr"
e(title) : "Linear regression"
e(vce) : "ols"
e(depvar) : "sexfreq1"
e(cmd) : "regress"
e(properties) : "b V"
e(predict) : "regres_p"
e(model) : "ols"
e(estat_cmd) : "regress_estat"
matrices:
e(b) : 1 x 3
e(V) : 3 x 3
functions:
e(sample)
.
. di e(mss)/(e(rss) + e(mss)) /* R-squared */
.14733375
. di (e(mss)/e(df_m))/(e(rss)/e(df_r)) /* F-statistic */
1699.9259
.
. matrix list e(b)
e(b)[1,3]
age agesqr _cons
y1 -.99396512 -.00494768 112.49477
. di _b[_cons] + (_b[age]*age) + (_b[agesqr]*agesqr)
50.427317
. generate yhat=_b[_cons] + (_b[age]*age) + (_b[agesqr]*agesqr)
(177 missing values generated)
.
.
. /* creturn */
. creturn list
System values
----------------------------------------------------------------------------------------------------------
c(current_date) = "30 Jan 2009"
c(current_time) = "12:34:31"
c(rmsg_time) = 0 (seconds, from set rmsg)
----------------------------------------------------------------------------------------------------------
c(stata_version) = 10.1
c(version) = 10.1 (version)
----------------------------------------------------------------------------------------------------------
c(born_date) = "11 Aug 2008"
c(flavor) = "IC"
c(SE) = 1
c(MP) = 0
c(processors) = 1 (Stata/MP, set processors)
c(processors_lic) = 1
c(processors_mach) = 1
c(processors_max) = 1
c(mode) = ""
c(console) = ""
----------------------------------------------------------------------------------------------------------
c(os) = "Windows"
c(osdtl) = ""
c(machine_type) = "PC"
c(byteorder) = "lohi"
c(username) = "HP_Owner"
----------------------------------------------------------------------------------------------------------
Directories and paths
----------------------------------------------------------------------------------------------------------
c(sysdir_stata) = "C:\Program Files\.." (sysdir)
c(sysdir_updates) = "C:\Program Files\.." (sysdir)
c(sysdir_base) = "C:\Program Files\.." (sysdir)
c(sysdir_site) = "C:\Program Files\.." (sysdir)
c(sysdir_plus) = "C:\Documents and .." (sysdir)
c(sysdir_personal) = "c:\ado\personal/" (sysdir)
c(sysdir_oldplace) = "c:\ado/" (sysdir)
c(tmpdir) = "C:\DOCUME~1\HP_Ow.."
----------------------------------------------------------------------------------------------------------
c(adopath) = "UPDATES;BASE;SITE.." (adopath)
c(pwd) = "C:\Documents and .." (cd)
c(dirsep) = "/"
----------------------------------------------------------------------------------------------------------
System limits
----------------------------------------------------------------------------------------------------------
c(max_N_theory) = 2147483647
c(max_k_theory) = 5000 (set maxvar)
c(max_width_theory) = 393204 (set maxvar)
----------------------------------------------------------------------------------------------------------
c(max_N_current) = 1638398 (set memory)
c(max_k_current) = 1070 (set memory)
c(max_width_current) = 1070 (set memory)
----------------------------------------------------------------------------------------------------------
c(max_matsize) = 11000
c(min_matsize) = 10
----------------------------------------------------------------------------------------------------------
c(max_macrolen) = 1081511
c(macrolen) = 165200 (set maxvar)
c(max_cmdlen) = 1081527
c(cmdlen) = 165216 (set maxvar)
c(namelen) = 32
----------------------------------------------------------------------------------------------------------
Numerical and string limits
----------------------------------------------------------------------------------------------------------
c(mindouble) = -8.9884656743e+307
c(maxdouble) = 8.9884656743e+307
c(epsdouble) = 2.22044604925e-16
c(smallestdouble) = 2.2250738585e-308
----------------------------------------------------------------------------------------------------------
c(minfloat) = -1.70141173319e+38
c(maxfloat) = 1.70141173319e+38
c(epsfloat) = 1.19209289551e-07
----------------------------------------------------------------------------------------------------------
c(minlong) = -2147483647
c(maxlong) = 2147483620
----------------------------------------------------------------------------------------------------------
c(minint) = -32767
c(maxint) = 32740
----------------------------------------------------------------------------------------------------------
c(minbyte) = -127
c(maxbyte) = 100
----------------------------------------------------------------------------------------------------------
c(maxstrvarlen) = 244
----------------------------------------------------------------------------------------------------------
Current dataset
----------------------------------------------------------------------------------------------------------
c(N) = 48797
c(k) = 14
c(width) = 24
c(changed) = 1
c(filename) = "C:\Documents and .."
c(filedate) = " 3 Nov 2007 13:59"
----------------------------------------------------------------------------------------------------------
Memory settings
----------------------------------------------------------------------------------------------------------
c(memory) = 52428800 (set memory)
c(maxvar) = 5000 (set maxvar)
c(matsize) = 400 (set matsize)
----------------------------------------------------------------------------------------------------------
Output settings
----------------------------------------------------------------------------------------------------------
c(more) = "off" (set more)
c(rmsg) = "off" (set rmsg)
c(dp) = "period" (set dp)
c(linesize) = 110 (set linesize)
c(pagesize) = 43 (set pagesize)
c(logtype) = "smcl" (set logtype)
----------------------------------------------------------------------------------------------------------
Interface settings
----------------------------------------------------------------------------------------------------------
c(dockable) = "on" (set dockable)
c(dockingguides) = "on" (set dockingguides)
c(locksplitters) = "off" (set locksplitters)
c(persistfv) = "off" (set persistfv)
c(persistvtopic) = "off" (set persistvtopic)
c(pinnable) = "on" (set pinnable)
c(reventries) = 5000 (set reventries)
c(xptheme) = "on" (set xptheme)
c(doublebuffer) = "on" (set doublebuffer)
c(linegap) = 1 (set linegap)
c(scrollbufsize) = 32000 (set scrollbufsize)
c(varlabelpos) = 32 (set varlabelpos)
c(maxdb) = 50 (set maxdb)
----------------------------------------------------------------------------------------------------------
Graphics settings
----------------------------------------------------------------------------------------------------------
c(graphics) = "on" (set graphics)
c(autotabgraphs) = "off" (set autotabgraphs)
c(scheme) = "s2color" (set scheme)
c(printcolor) = "automatic" (set printcolor)
c(copycolor) = "automatic" (set copycolor)
----------------------------------------------------------------------------------------------------------
Efficiency settings
----------------------------------------------------------------------------------------------------------
c(adosize) = 1000 (set adosize)
c(virtual) = "off" (set virtual)
----------------------------------------------------------------------------------------------------------
Network settings
----------------------------------------------------------------------------------------------------------
c(checksum) = "off" (set checksum)
c(timeout1) = 120 (set timeout1)
c(timeout2) = 300 (set timeout2)
----------------------------------------------------------------------------------------------------------
c(httpproxy) = "off" (set httpproxy)
c(httpproxyhost) = "" (set httpproxyhost)
c(httpproxyport) = 80 (set httpproxyport)
----------------------------------------------------------------------------------------------------------
c(httpproxyauth) = "off" (set httpproxyauth)
c(httpproxyuser) = "" (set httpproxyuser)
c(httpproxypw) = "" (set httpproxypw)
----------------------------------------------------------------------------------------------------------
Update settings
----------------------------------------------------------------------------------------------------------
c(update_query) = "on" (set update_query)
c(update_interval) = 7 (set update_interval)
c(update_prompt) = "on" (set update_prompt)
----------------------------------------------------------------------------------------------------------
Trace (program debugging) settings
----------------------------------------------------------------------------------------------------------
c(trace) = "off" (set trace)
c(tracedepth) = 32000 (set tracedepth)
c(tracesep) = "on" (set tracesep)
c(traceindent) = "on" (set traceindent)
c(traceexpand) = "on" (set traceexpand)
c(tracenumber) = "off" (set tracenumber)
c(tracehilite) = "" (set tracehilite)
----------------------------------------------------------------------------------------------------------
Mata settings
----------------------------------------------------------------------------------------------------------
c(matastrict) = "off" (set matastrict)
c(matalnum) = "off" (set matalnum)
c(mataoptimize) = "on" (set mataoptimize)
c(matafavor) = "space" (set matafavor)
c(matacache) = 400 (set matacache)
c(matalibs) = "lmatabase;lmataad.." (set matalibs)
c(matamofirst) = "off" (set matamofirst)
----------------------------------------------------------------------------------------------------------
Other settings
----------------------------------------------------------------------------------------------------------
c(type) = "float" (set type)
c(level) = 95 (set level)
c(maxiter) = 16000 (set maxiter)
c(searchdefault) = "local" (set searchdefault)
c(seed) = "X075bcd151f123bb5.." (set seed)
c(varabbrev) = "on" (set varabbrev)
----------------------------------------------------------------------------------------------------------
Other
----------------------------------------------------------------------------------------------------------
c(pi) = 3.141592653589793
c(alpha) = "a b c d e f g h i.."
c(ALPHA) = "A B C D E F G H I.."
c(Mons) = "Jan Feb Mar Apr M.."
c(Months) = "January February .."
c(Wdays) = "Sun Mon Tue Wed T.."
c(Weekdays) = "Sunday Monday Tue.."
c(rc) = 0 (capture)
----------------------------------------------------------------------------------------------------------
. di c(current_date)
30 Jan 2009
.
.
. quietly regress sexfreq1 age
. local rsqf e(r2)
. local rsqv = e(r2)
. di `rsqf' /* This has the current R-squared */
.14679606
. di `rsqv' /* as does this */
.14679606
.
. quietly regress sexfreq1 age sex race
. di `rsqf' /* the formula has the new R-squared */
.15021663
. di `rsqv' /* this guy has the old one */
.14679606
.
.
. foreach var of varlist age educ sexfreq1 {
2. summarize `var'
3. }
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
age | 48620 45.75072 17.52199 18 89
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
educ | 48650 12.67239 3.145478 0 20
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
sexfreq1 | 19711 55.69299 65.99505 0 237
.
. local varsum "age educ sexfreq1"
. foreach var of varlist `varsum' {
2. summarize `var'
3. }
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
age | 48620 45.75072 17.52199 18 89
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
educ | 48650 12.67239 3.145478 0 20
Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
sexfreq1 | 19711 55.69299 65.99505 0 237
.
. foreach var of varlist socrel socommun socfrend socbar {
2. recode `var' (1=1) (2/7=0) (else=.), gen(`var'1)
3. }
(25366 differences between socrel and socrel1)
(26104 differences between socommun and socommun1)
(26867 differences between socfrend and socfrend1)
(27295 differences between socbar and socbar1)
.
.
. /* Problems */
.
. log close
log: C:\Documents and Settings\HP_Owner\My Documents\My Classes\SOCY498C\Macros\SOCY498C-HW7-MACROS.tx
> t
log type: text
closed on: 30 Jan 2009, 12:34:33
--------------------------------------------------------------------------------------------------------------