function(df){
    df[,8] <- df[,1]
    df[,9] <- df[,1]
    l5 <- China35$`Hts Num`
    l10 <- China310$`Hts Num`
    l20 <- China320$`Hts Num`
    l25 <- China325$`Hts Num`
    listoflists <- list(l5,l10,l20,l25)
    c <- c(0,0,0,0)
    r <- c(0.05,0.1,0.2,0.25)
    r0 <- c(0.05,0.05,0.1,0.1)
    for (i in 1:3600){
        code <- df[i,1]
        for (j in 1:4){
            if (code %in% listoflists[[j]]){
                c[j] <- 1}}
        rates <- c*r
        rates <- rates[which(rates != 0)]
        df[i,9] <- as.numeric(mean(rates))
        rates0 <- c*r0
        rates0 <- rates0[which(rates0 != 0)]
        df[i,8] <- as.numeric(mean(rates0))}
        df
    }

> Data2 <- addrates_avg(Data1)
> Data2[500,1]
[1] "2005.10"
> Data2[500,8]
[1] "0.075"
> Data2[500,9]
[1] "0.15"

> df <- Data1
> df[,8] <- df[,1]
> df[,9] <- df[,1]
> l5 <- China35$`Hts Num`
> l10 <- China310$`Hts Num`
> l20 <- China320$`Hts Num`
> l25 <- China325$`Hts Num`
> listoflists <- list(l5,l10,l20,l25)
> c <- c(0,0,0,0)
> r <- c(0.05,0.1,0.2,0.25)
> r0 <- c(0.05,0.05,0.1,0.1)
> i = 500
> code <- df[i,1]
> code
[1] "2005.10"
> for (j in 1:4){
+     if (code %in% listoflists[[j]]){
+         c[j] <- 1}}
> c
[1] 0 1 0 0
> rates <- c*r
> rates <- rates[which(rates != 0)]
> rates
[1] 0.1
> df[i,9] <- as.numeric(mean(rates))
> rates0 <- c*r0
> rates0
[1] 0.00 0.05 0.00 0.00
> rates0 <- rates0[which(rates0 != 0)]
> df[i,8] <- as.numeric(mean(rates0))
> df[500,8]
> 0.05
> df[500,9]
> 0.1

> df <- Data1[c(400:505),]
> df <- addrates_avg(df)
> df$V8[which(df$X1 == "2005.10")]
[1] "0.075"

> df <- Data1[c(450:505),]
> df <- addrates_avg(df)
> df$V8[which(df$X1 == "2005.10")]
[1] "0.075"

> df <- Data1[c(475:505),]
> df <- addrates_avg(df)
> df$V8[which(df$X1 == "2005.10")]
[1] "0.0833333333333333"

> df <- Data1[c(490:505),]
> df <- addrates_avg(df)
> df$V8[which(df$X1 == "2005.10")]
[1] "0.0833333333333333"

> df <- Data1[c(495:505),]
> df <- addrates_avg(df)
> df$V8[which(df$X1 == "2005.10")]
[1] "0.0833333333333333"

> df <- Data1[c(499:502),]
> df <- addrates_avg(df)
> df$V8[which(df$X1 == "2005.10")]
[1] "0.05"