感知压力数据
The total score of PSS is obtained by reversing the scores on the positive items and then summing across all the items, with a higher score indicating higher perceived stress.
PSS的总分是通过将阳性项目的得分反转，然后在所有项目上求和而获得的，得分越高，表示感知到的压力越高
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8369204/学习文章1
学习文章2-https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7082906/
# 先把depress的提取出来
library(haven)
library(dplyr)
# wave是2015的
depress <- read_sas("C:/Users/xxx/Desktop/CHNS/Master_PE_PA_201908/pstress_12.sas7bdat")
names(depress)depress_14 <- c("U551" , "U552" , "U553" ,"U554" , "U555" , "U556" , "U557" , "U558" , "U559" , "U560", "U561" ,"U562" , "U563" , "U564"
)
# 只去了两个
depress <- tidyr::drop_na(depress, all_of(depress_14))qaq= "U551"
for (qaq in depress_14) {indice = depress[[qaq]] == 9depress[[qaq]][indice] <- NA
}count(depress,U560) # 实际没有0# second round del 9
depress <- tidyr::drop_na(depress, all_of(depress_14))length(unique(depress$Idind))
# negative subscale (items 1,2,3,8,11,12 and 14) and
# # positive subscale (items 4,5,6,7,9,10 and 13)
# 0, 1, 2, 3, 4
# 4, 3, 2, 1, 0
nni <- c("U554","U555","U556","U557", "U559","U560", "U563")
depress$U554[1:20]
count(depress,U554 )
for (ni in nni) {depress[[ni]] <- dplyr::case_when(depress[[ni]] == 1 ~ 5,depress[[ni]] == 2 ~ 4,depress[[ni]] == 3 ~ 3,depress[[ni]] == 4 ~ 2,depress[[ni]] == 5 ~ 1,depress[[ni]] == 9 ~ NA,)
}
count(depress,U554)depress$d_score <- depress$U551+depress$U552+depress$U553+depress$U554+depress$U555+depress$U556+depress$U557+depress$U558+depress$U559+depress$U560+depress$U561+depress$U562+depress$U563+ depress$U564 - 14
summary(depress$d_score)quantiles <- quantile(depress$d_score, probs = c(0, 1/4, 2/4, 3/4, 1),na.rm = T )
depress$d_score_Cat <- cut(depress$d_score,breaks = quantiles, include.lowest = T,labels = c("q1","q2","q3","q4"))
count(depress, d_score_Cat) length(unique(depress$Idind))
table(depress$wave) # 只有2015年的
