Code for Saddaf Akhtar's email inquiry

Saddaf_lx_smooth.R

library(tidyverse)
library(DemoTools)
library(readxl)
IN <-read_excel("Data/Saddaf_New_file.xlsx")|> 
  mutate(age = parse_number(Age),
         lx = parse_number(lx))

lx2qx <- function(lx){
  -diff(c(lx,0)) / lx
}

IN |> 
  ggplot(aes(x=age,y=lx,color=Sex,group=Sex)) +
  geom_line() +
  facet_wrap(Time~State+Residence) +
  labs(title = "hard to see which might need smoothing")

IN |> 
  group_by(Time,State,Residence,Sex) |> 
  mutate(qx = lx2qx(lx)) |> 
  ungroup() |> 
  ggplot(aes(x=age,y=qx,color=Sex,group=Sex)) +
  geom_line() +
  facet_wrap(Time~State+Residence) +
  scale_y_log10() +
  labs(title = "seems like only Delhi rural women need smoothing")