// 清除之前的数据和变量
clear all
 
// 输入数据
input str10 race str15 age_group owns_cat does_not_own_cat
"White" "youth" 65 84
"White" "middle-aged" 92 107
"White" "elderly" 87 120
"Black" "youth" 4 14
"Black" "middle-aged" 3 17
"Black" "elderly" 4 15
"Other" "youth" 7 13
"Other" "middle-aged" 4 19
"Other" "elderly" 4 12
end
 
// 计算每个组别的拥有猫和不拥有猫的百分比
gen total = owns_cat + does_not_own_cat
gen owns_cat_percentage = (owns_cat / total) * 100
gen does_not_own_cat_percentage = (does_not_own_cat / total) * 100
 
// 重构数据为长格式
gen cat_status = "owns_cat"
gen percentage = owns_cat_percentage
tempfile temp
save `temp', replace
 
gen cat_status = "does_not_own_cat"
gen percentage = does_not_own_cat_percentage
append using `temp'
 
// 将 cat_status 转换为标签变量
replace cat_status = "Owns Cats" if cat_status == "owns_cat"
replace cat_status = "Does Not Own Cats" if cat_status == "does_not_own_cat"
label define cat_status_lbl 1 "Owns Cats" 2 "Does Not Own Cats"
encode cat_status, gen(cat_status_encoded)
label values cat_status_encoded cat_status_lbl
 
// 添加对应的数量
gen population = .
replace population = owns_cat if cat_status == "Owns Cats"
replace population = does_not_own_cat if cat_status == "Does Not Own Cats"
 
// 绘制百分比的堆叠柱状图
graph bar percentage, over(race) over(age_group) by(cat_status_encoded, legend(col(1))) ///
    stack title("Cat Ownership Percentage by Race and Age Group") ///
    ytitle("Percentage") xtitle("Race and Age Group") ///
    blabel(bar, format(%4.1f) position(center) size(small))
 

Ly8g5riF6Zmk5LmL5YmN55qE5pWw5o2u5ZKM5Y+Y6YePCmNsZWFyIGFsbAoKLy8g6L6T5YWl5pWw5o2uCmlucHV0IHN0cjEwIHJhY2Ugc3RyMTUgYWdlX2dyb3VwIG93bnNfY2F0IGRvZXNfbm90X293bl9jYXQKIldoaXRlIiAieW91dGgiIDY1IDg0CiJXaGl0ZSIgIm1pZGRsZS1hZ2VkIiA5MiAxMDcKIldoaXRlIiAiZWxkZXJseSIgODcgMTIwCiJCbGFjayIgInlvdXRoIiA0IDE0CiJCbGFjayIgIm1pZGRsZS1hZ2VkIiAzIDE3CiJCbGFjayIgImVsZGVybHkiIDQgMTUKIk90aGVyIiAieW91dGgiIDcgMTMKIk90aGVyIiAibWlkZGxlLWFnZWQiIDQgMTkKIk90aGVyIiAiZWxkZXJseSIgNCAxMgplbmQKCi8vIOiuoeeul+avj+S4que7hOWIq+eahOaLpeacieeMq+WSjOS4jeaLpeacieeMq+eahOeZvuWIhuavlApnZW4gdG90YWwgPSBvd25zX2NhdCArIGRvZXNfbm90X293bl9jYXQKZ2VuIG93bnNfY2F0X3BlcmNlbnRhZ2UgPSAob3duc19jYXQgLyB0b3RhbCkgKiAxMDAKZ2VuIGRvZXNfbm90X293bl9jYXRfcGVyY2VudGFnZSA9IChkb2VzX25vdF9vd25fY2F0IC8gdG90YWwpICogMTAwCgovLyDph43mnoTmlbDmja7kuLrplb/moLzlvI8KZ2VuIGNhdF9zdGF0dXMgPSAib3duc19jYXQiCmdlbiBwZXJjZW50YWdlID0gb3duc19jYXRfcGVyY2VudGFnZQp0ZW1wZmlsZSB0ZW1wCnNhdmUgYHRlbXAnLCByZXBsYWNlCgpnZW4gY2F0X3N0YXR1cyA9ICJkb2VzX25vdF9vd25fY2F0IgpnZW4gcGVyY2VudGFnZSA9IGRvZXNfbm90X293bl9jYXRfcGVyY2VudGFnZQphcHBlbmQgdXNpbmcgYHRlbXAnCgovLyDlsIYgY2F0X3N0YXR1cyDovazmjaLkuLrmoIfnrb7lj5jph48KcmVwbGFjZSBjYXRfc3RhdHVzID0gIk93bnMgQ2F0cyIgaWYgY2F0X3N0YXR1cyA9PSAib3duc19jYXQiCnJlcGxhY2UgY2F0X3N0YXR1cyA9ICJEb2VzIE5vdCBPd24gQ2F0cyIgaWYgY2F0X3N0YXR1cyA9PSAiZG9lc19ub3Rfb3duX2NhdCIKbGFiZWwgZGVmaW5lIGNhdF9zdGF0dXNfbGJsIDEgIk93bnMgQ2F0cyIgMiAiRG9lcyBOb3QgT3duIENhdHMiCmVuY29kZSBjYXRfc3RhdHVzLCBnZW4oY2F0X3N0YXR1c19lbmNvZGVkKQpsYWJlbCB2YWx1ZXMgY2F0X3N0YXR1c19lbmNvZGVkIGNhdF9zdGF0dXNfbGJsCgovLyDmt7vliqDlr7nlupTnmoTmlbDph48KZ2VuIHBvcHVsYXRpb24gPSAuCnJlcGxhY2UgcG9wdWxhdGlvbiA9IG93bnNfY2F0IGlmIGNhdF9zdGF0dXMgPT0gIk93bnMgQ2F0cyIKcmVwbGFjZSBwb3B1bGF0aW9uID0gZG9lc19ub3Rfb3duX2NhdCBpZiBjYXRfc3RhdHVzID09ICJEb2VzIE5vdCBPd24gQ2F0cyIKCi8vIOe7mOWItueZvuWIhuavlOeahOWghuWPoOafseeKtuWbvgpncmFwaCBiYXIgcGVyY2VudGFnZSwgb3ZlcihyYWNlKSBvdmVyKGFnZV9ncm91cCkgYnkoY2F0X3N0YXR1c19lbmNvZGVkLCBsZWdlbmQoY29sKDEpKSkgLy8vCiAgICBzdGFjayB0aXRsZSgiQ2F0IE93bmVyc2hpcCBQZXJjZW50YWdlIGJ5IFJhY2UgYW5kIEFnZSBHcm91cCIpIC8vLwogICAgeXRpdGxlKCJQZXJjZW50YWdlIikgeHRpdGxlKCJSYWNlIGFuZCBBZ2UgR3JvdXAiKSAvLy8KICAgIGJsYWJlbChiYXIsIGZvcm1hdCglNC4xZikgcG9zaXRpb24oY2VudGVyKSBzaXplKHNtYWxsKSkK