Início / coding / Perguntas / 78946168
Accepted
James Rider
Asked: 2024-09-04 05:47:36 +0800 CST

Alteração percentual total em uma coluna em R?

  • 772

Tenho um conjunto de dados de painel em R que documenta a mudança percentual na população de cidades ao longo de vários anos. Uma versão simplificada se parece com isso:

city_pop<-tibble(city=c("NYC", "NYC", "LA", "LA"), 
       percentpopchange=c("0.1", "0.2", "0.5", "0.3"), 
       year=c(1980, 1981, 1980, 1981))

Gostaria de calcular a mudança populacional total (em porcentagem) para NYC e LA em 1980 e 1981. Por exemplo, seria 0,32 para NYC. Meu conjunto de dados ideal com tal saída seria:

cum_city_pop<-tibble(city=c("NYC", "LA"), 
       total_percentpopchange=c("0.32", "0.95"))

Existe alguma maneira de fazer isso em R? Eu agradeceria muito qualquer ajuda.

2 respostas

  1. Best Answer

    Supondo que seu 0,93 para LA seja um erro de digitação:

    library(dplyr)

city_pop <- tibble(
  city = c("NYC", "NYC", "LA", "LA"), 
  percentpopchange = c("0.1", "0.2", "0.5", "0.3"), 
  year = c(1980, 1981, 1980, 1981)
)

city_pop %>% 
  group_by(city) %>% 
  summarize(
    total_percentpopchange = as.character(
      prod(1 + as.numeric(percentpopchange)) - 1
  ))
#> # A tibble: 2 × 2
#>   city  total_percentpopchange
#>   <chr> <chr>                 
#> 1 LA    0.95                  
#> 2 NYC   0.32

    Criado em 2024-09-03 com reprex v2.1.1

    • 3

  2. Se assumirmos que a população de NY era de 1.000 em 1980, então sua população muda da seguinte forma após 1980 e 1981:

    1980: $ 1.000 + $ 1.000 * 0,1 = $ 1.100

    1981: $ 1.100 + $ 1.100 * 0,2 = $ 1.320.

    Então, ao longo de 1980 e 1981, a população de NY aumentou por um fator de: ($1320-$1000)/$1000 = 0,32.

    Um pouco de código R que calcula isso está colado abaixo:

    city_pop<-tibble(city=c("NYC", "NYC", "LA", "LA"), 
                 percentpopchange=c("0.1", "0.2", "0.5", "0.3"), 
                 year=c(1980, 1981, 1980, 1981))
city_pop$percentpopchange <- as.numeric(city_pop$percentpopchange)

cum_city_pop_final <- tibble(city=unique(city_pop$city))

unique_cities <- unique(city_pop$city)
total_percent_pop_changes <- c()
for (i in 1:length(unique(city_pop$city))){
  pop <- 100
  for (j in 1:length(unique(city_pop$year))){
    pop <- pop + pop * city_pop$percentpopchange[which(city_pop$city==unique(city_pop$city)[i] & city_pop$year==unique(city_pop$year)[j])]
    if (j==length(unique(city_pop$year))){
      total_percent_pop_changes <- c(total_percent_pop_changes, (pop-100)/100)
    }
  }
}

cum_city_pop_final$total_percentpopchange <- total_percent_pop_changes

    A propósito, acho que a mudança populacional percentual total de Los Angeles entre 1980 e 1981 seria de 0,95, não de 0,93.

    • 0

relate perguntas