Eu tenho uma pergunta simples. Eu gostaria de somar duas distribuições não paramétricas.

Aqui está um exemplo. Existem duas cidades que têm 10 casas. sabemos o consumo de energia para cada casa. (editado) Quero obter a distribuição de probabilidade da soma de uma casa aleatória escolhida de cada cidade.

A1 <- c(1,2,3,3,3,4,4,5,6,7) #10 houses' energy consumption for city A
B1 <- c(11,13,15,17,17,18,18,19,20,22) #10 houses' energy consumption for city B

Eu tenho uma distribuição de probabilidade de A1 e B1, como posso obter a distribuição de probabilidade de A1 + B1? Se eu apenas usarA1+B1 em R, dá12 15 18 20 20 22 22 24 26 29. No entanto, não acho que isso esteja certo. Porque não há ordem nas casas.

Quando mudo a ordem das casas, dá outros resultados.

# Original
A1 <- c(1,2,3,3,3,4,4,5,6,7)
B1 <- c(11,13,15,17,17,18,18,19,20,22)
#change order 1
A2 <- c(7,6,5,4,4,3,3,3,2,1) 
B2 <- c(22,20,19,18,18,17,17,15,13,11)
#change order 2
A3 <- c(3,3,3,4,4,5,6,7,1,2) 
B3 <- c(17,17,18,18,19,13,20,11,22,15)
sum1 <- A1+B1; sum1
sum2 <- A1+B2; sum2
sum3 <- A3+B3; sum3

As linhas vermelhas são soma1, soma2 e soma3. Não tenho certeza de como posso obter a distribuição da soma de duas distribuições.Por favor, dê-me algumas ideias.Obrigado!

(Se essas distribuições forem distribuições normais ou uniformes, eu poderia obter a soma da distribuição facilmente, mas elas não são normais e não há pedido)