Skip to content

Commit

Permalink
s
Browse files Browse the repository at this point in the history
  • Loading branch information
@friosavila
friosavila committed Mar 12, 2024
1 parent 332f1c5 commit 279e978
Show file tree
Hide file tree
Showing 2 changed files with 84 additions and 63 deletions.
32 changes: 16 additions & 16 deletions csdid3/csdid.mata
View file
Original file line number Diff line number Diff line change
Expand Up @@ -271,7 +271,7 @@ void csdid::csdid_setup(){
foid = uniqrows(oid)
}
if (length(wvar)==0) wvar=J(rows(yvar),1,1)
//if (length(wvar)==0) wvar=J(rows(yvar),1,1)
}
// Justs Puts all into Running Data
void csdid::makeid(){
Expand Down Expand Up @@ -465,7 +465,7 @@ void csdid::csdid(){
// To make it faster
real matrix index1
frwt=frif=J( ((type_data==1) ? max(oid) : rows(oid)) ,
frif=J( ((type_data==1) ? max(oid) : rows(oid)) ,
rows(fgtvar),.)
index1 = range(1,rows(oid),1)
Expand Down Expand Up @@ -504,17 +504,17 @@ void csdid::csdid(){
convar[i,]=1
if (shortx==0) frif[drdid.oid,i]=drdid.rif:*sign(eventvar[i]+.01)
else frif[drdid.oid,i]=drdid.rif
frwt[drdid.oid,i]=drdid.wtrt
//frwt[drdid.oid,i]=drdid.wtrt
}
dots = dots-convar[i,]
stata(sprintf("_dots %f %f",i,dots))
}
/// fixes missing in rifs
fixrif()
_editmissing(frwt,0)
//_editmissing(frwt,0)
/// Extra clean up
frwt = select(frwt , convar')
//frwt = select(frwt , convar')
frif = select(frif , convar')
eventvar = select(eventvar, convar)
fgtvar = select(fgtvar , convar)
Expand All @@ -529,8 +529,8 @@ void csdid::csdid(){
if (type_data==1) {
real matrix aux
if (length(cvar)>0) aux = oid, gvar, wvar, cvar
else aux = oid, gvar, wvar
if (length(cvar)>0) aux = oid, gvar, wvar, cvar
else aux = oid, gvar, wvar
aux=uniqrows(aux)
oid = aux[,1]
gvar= aux[,2]
Expand All @@ -539,15 +539,15 @@ void csdid::csdid(){
}
/// Very last step. Sort important variables by Cvar?
if (length(cvar)>0) {
ord = order( (cvar,oid), (1,2) )
oid = oid[ord,]
cvar = cvar[ord,]
gvar = gvar[ord,]
wvar = wvar[ord,]
frif = frif[ord,]
frwt = frwt[ord,]
}
//if (length(cvar)>0) {
// ord = order( (cvar,oid), (1,2) )
// oid = oid[ord,]
// cvar = cvar[ord,]
// gvar = gvar[ord,]
// wvar = wvar[ord,]
// frif = frif[ord,]
// frwt = frwt[ord,]
//}
aux = J(0,0,.)
}
Expand Down
115 changes: 68 additions & 47 deletions csdid3/csdid_stats.mata
View file
Original file line number Diff line number Diff line change
Expand Up @@ -88,13 +88,15 @@ void csdid_estat::bvcv_clus(real matrix rif,
nobs= rows(rif)
// sort //ord = order(cvar,1) //rif = rif[ord,] //cvar= cvar[ord,]
// Standard Errors
cvar = csdid.cvar[csdid.sortcvar,]
rif = rif[csdid.sortcvar,]
// cvar = csdid.cvar
//[csdid.sortcvar,]
//rif = rif
//[csdid.sortcvar,]
// Consider Cleaning data based on Rows. How? Would it be more Time consuming?
// need to do this at an earlier stage.
// Smaller sample may have faster Run.
info = panelsetup(cvar,1)
info = panelsetup(csdid.cvar,1)
nclust= rows(info)
real matrix sumrif
sumrif= panelsum(rif:-bb,info)
Expand All @@ -105,12 +107,14 @@ void csdid_estat::bvcv_clus(real matrix rif,
real matrix csdid_estat::aggte(real matrix rif,| real matrix wgt ) {
real matrix mn_all, mn_rif, mn_wgt
real matrix mn_all, mn_rif, mn_wgt, bwgt
if (args()==2) {
//wgt = J(1,cols(rif),1)
// Avg Effect
bwgt = wgt:*(rif:!=.)
fixrif(rif)
mn_rif = mean(rif)
mn_wgt = mean(wgt)
mn_wgt = mean(bwgt)
mn_all = sum(mn_rif:*mn_wgt):/sum(mn_wgt)
// gets agg rif
real matrix wgtw, attw
Expand All @@ -119,12 +123,13 @@ real matrix csdid_estat::aggte(real matrix rif,| real matrix wgt ) {
// r1 r2 r3
real matrix r1 , r2 , r3
r1 = (wgtw:*(rif :-mn_rif))
r2 = (attw:*(wgt :-mn_wgt ))
r3 = (wgt :- mn_wgt) :* (mn_all :/ sum(mn_wgt) )
r2 = (attw:*(bwgt :-mn_wgt ))
r3 = (bwgt :- mn_wgt) :* (mn_all :/ sum(mn_wgt) )
// Aggregates into 1
return(rowsum(r1):+rowsum(r2):-rowsum(r3):+mn_all)
}
else {
fixrif(rif)
mn_rif = mean(rif)
mn_all = mean(mn_rif')
// gets agg rif
Expand All @@ -133,17 +138,16 @@ real matrix csdid_estat::aggte(real matrix rif,| real matrix wgt ) {
// r1 r2 r3
real matrix r1 , r2 , r3
r1 = ((rif :-mn_rif):/cols(mn_rif))
r2 = (mn_rif:*(wgt :-mn_wgt ):/cols(mn_rif))
//r3 = (wgt :- mn_wgt) :* (mn_all :/ sum(mn_wgt) )
// Aggregates into 1
return(rowsum(r1):+rowsum(r2):+mn_all)
return(rowsum(r1):+mn_all)
}
}
//
real matrix csdid_estat::spaggte(class csdid matrix csdid ,
real matrix toselect ) {
real matrix toselect ) {
// csdid contains the Index for spsdid
// spsdid contains the RIFs
Expand All @@ -152,35 +156,52 @@ real matrix csdid_estat::spaggte(class csdid matrix csdid ,
// How many Cols will be needed. Ideally less than FULL sample
ntoselect = cols(toselect)
mn_wgt = mn_rif = J(1,ntoselect,0)
// over all toselect.
for(i=1;i<=ntoselect;i++) {
mn_wgt[,i] = csdid.spcsdid[i].mn_wattgt
mn_rif[,i] = csdid.spcsdid[i].mn_attgt
}
mn_all=mean(mn_rif',mn_wgt')
// Stays as is
real matrix wgtw, attw
wgtw = (mn_wgt ) :/sum(mn_wgt)
attw = (mn_rif ) :/sum(mn_wgt)
real matrix rr1
rr1=J(csdid.nobs,1,0)
//rr2=rr3 = J(10,1,0)
for(i=1;i<=ntoselect;i++) {
rr1[csdid.spcsdid[i].index] = rr1[csdid.spcsdid[i].index]:+ wgtw[i]:*(csdid.spcsdid[i].attgt )
arr2 = J(csdid.nobs , 1 , - csdid.spcsdid[i].mn_wattgt);
arr2[csdid.spcsdid[i].index] = (csdid.spcsdid[i].wattgt )
rr1 = rr1 :+ arr2*(attw[i]:-(mn_all:/ sum(mn_wgt)))
}
//*
// Aggregates into 1
**class spaggte scalar aux_aggte
//aux_aggte.aggte = rr1
//aux_aggte.mn_aggte = mn_all
return(rr1:+mn_all)
if (length(csdid.wvar)>0) {
// over all toselect.
for(i=1;i<=ntoselect;i++) {
mn_wgt[,i] = csdid.spcsdid[i].mn_wattgt
mn_rif[,i] = csdid.spcsdid[i].mn_attgt
}
mn_all=mean(mn_rif',mn_wgt')
// Stays as is
real matrix wgtw, attw
wgtw = (mn_wgt ) :/sum(mn_wgt)
attw = (mn_rif ) :/sum(mn_wgt)
real matrix rr1
rr1=J(csdid.nobs,1,0)
//rr2=rr3 = J(10,1,0)
for(i=1;i<=ntoselect;i++) {
rr1[csdid.spcsdid[i].index] = rr1[csdid.spcsdid[i].index]:+ wgtw[i]:*(csdid.spcsdid[i].attgt )
arr2 = J(csdid.nobs , 1 , - csdid.spcsdid[i].mn_wattgt);
arr2[csdid.spcsdid[i].index] = (csdid.spcsdid[i].wattgt )
rr1 = rr1 :+ arr2*(attw[i]:-(mn_all:/ sum(mn_wgt)))
}
return(rr1:+mn_all)
}
else {
for(i=1;i<=ntoselect;i++) {
// mn_wgt[,i] = csdid.spcsdid[i].mn_wattgt
mn_rif[,i] = csdid.spcsdid[i].mn_attgt
}
mn_all=mean(mn_rif')
// Stays as is
real matrix wgtw, attw
attw = (mn_rif ) :/cols(mn_rif)
real matrix rr1
rr1=J(csdid.nobs,1,0)
//rr2=rr3 = J(10,1,0)
for(i=1;i<=ntoselect;i++) {
rr1[csdid.spcsdid[i].index] = rr1[csdid.spcsdid[i].index]:+ (1:/cols(mn_rif)):*(csdid.spcsdid[i].attgt )
}
return(rr1:+mn_all)
}
}
// Only for Weights (if needed)
real matrix csdid_estat::spwaggte(class csdid matrix csdid ,
Expand All @@ -205,7 +226,7 @@ void csdid_estat::atts_asym(class csdid scalar csdid){
error = 0
if (test_type==1) {
// ATTGT
erif_attgt()
erif_attgt(csdid)
//=select(csdid.frif,select_data(csdid)')
if (length(csdid.cvar)==0) {
bvcv_asym(erif)
Expand Down Expand Up @@ -284,14 +305,14 @@ void csdid_estat::make_table(){
}
// This fixes ERIF
void csdid_estat::fixrif(real matrix erif){
void csdid_estat::fixrif(real matrix arif){
real matrix mn_rif, rif2
real scalar cnmiss
cnmiss = colnonmissing(erif)
mn_rif= colsum(erif):/cnmiss
erif = erif:-mn_rif
erif = editmissing(erif,0)
erif = mn_rif:+erif:*(rows(erif):/cnmiss)
cnmiss = colnonmissing(arif)
mn_rif= colsum(arif):/cnmiss
arif = arif:-mn_rif
_editmissing(arif,0)
arif = mn_rif:+arif:*(rows(arif):/cnmiss)
}
// erif should have only the important ERIFs
// This Saves
Expand Down Expand Up @@ -537,7 +558,7 @@ void csdid_estat::pretrend(class csdid scalar csdid ){
if (sum(toselect)>0) {
if (csdid.sparse == 1) {
// gets all RIFs
erif_attgt()
erif_attgt(csdid)
if (length(csdid.cvar)==0) bvcv_asym(erif)
else bvcv_clus(erif,csdid.cvar)
}
Expand Down Expand Up @@ -570,7 +591,7 @@ void csdid_estat::pretrend2(class csdid scalar csdid ){
if (sum(toselect)>0) {
if (csdid.sparse == 1) {
// gets all RIFs
erif_attgt()
erif_attgt(csdid)
if (length(csdid.cvar)==0) bvcv_asym(erif)
else bvcv_clus(erif,csdid.cvar)
}
Expand Down

0 comments on commit 279e978

Please sign in to comment.