私は非常に単純な if else ループを作成して、別の変数を四分位数にビン化する新しい変数を作成しています。これは非常に単純な手順のように見えますが、ループはすべてのデータを中央値と 3 番目の四分位数 (四分位数の定義に違反しています) にグループ化します。
これが私のデータの構造です：
> str(tmp)
Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 435 obs. of 12 variables:
$ CD112FP : chr "01" "02" "03" "04" ...
$ State : chr "ALABAMA" "ALABAMA" "ALABAMA" "ALABAMA" ...
$ Year : num 2011 2011 2011 2011 2011 ...
$ Alignment : num 0 0 0 0 0 0 1 0 0 0 ...
$ State_Aligned : num 0 0 0 0 0 0 0 1 0 0 ...
$ PercentFunding : num 0.0658 0.29 0.6764 0.0174 0.047 ...
$ fips : chr "01" "01" "01" "01" ...
$ ssa : int 1 1 1 1 1 1 1 NA 3 3 ...
$ region : int 3 3 3 3 3 3 3 NA 4 4 ...
$ division : int 6 6 6 6 6 6 6 NA 8 8 ...
$ abb : chr "AL" "AL" "AL" "AL" ...
$ PercentFundingBinned: chr "0.0625-0.1799" "0.0625-0.1799" "0.0625-0.1799" "0.0625-0.1799" ...
これは私のデータの頭です：
head(tmp)
# A tibble: 6 x 12
CD112FP State Year Alignment State_Aligned PercentFunding fips ssa region division abb PercentFundingBinned
<chr> <chr> <dbl> <dbl> <dbl> <dbl> <chr> <int> <int> <int> <chr> <chr>
1 01 ALABAMA 2011 0 0 0.0658 01 1 3 6 AL 0.0625-0.1799
2 02 ALABAMA 2011 0 0 0.290 01 1 3 6 AL 0.0625-0.1799
3 03 ALABAMA 2011 0 0 0.676 01 1 3 6 AL 0.0625-0.1799
4 04 ALABAMA 2011 0 0 0.0174 01 1 3 6 AL 0.0625-0.1799
5 05 ALABAMA 2011 0 0 0.0470 01 1 3 6 AL 0.0625-0.1799
6 06 ALABAMA 2011 0 0 0.0440 01 1 3 6 AL 0.0625-0.1799
次の if else ループを使用しています。
tmp$PercentFundingBinned <- NULL
if (tmp$PercentFunding >= quantile(tmp$PercentFunding, 0.75)) {
tmp$PercentFundingBinned <- paste0(round(quantile(tmp$PercentFunding, 0.75), 4), "-",
round(max(tmp$PercentFundingBinned), 4))
} else if (tmp$PercentFunding >= median(tmp$PercentFunding)){
tmp$PercentFundingBinned <- paste0(round(median(tmp$PercentFunding),4), "-",
round(quantile(tmp$PercentFunding, 0.75),4))
} else if (tmp$PercentFunding >= quantile(tmp$PercentFunding, 0.25)){
tmp$PercentFundingBinned <- paste0(round(quantile(tmp$PercentFunding, 0.25),4), "-",
round(median(tmp$PercentFunding),4))
} else {
tmp$PercentFundingBinned <- paste0(round(min(tmp$PercentFunding),4), "-",
round(quantile(tmp$PercentFunding, 0.25),4))
}
次のカテゴリを返します。
unique(tmp$PercentFundingBinned)
[1] "0.0625-0.1799"
何をすべきか、どのように適合するかわかりません。これは本当に簡単な手順のようです。どんなアドバイスも役に立ちます、ありがとう！