sábado, 31 de maio de 2014

Programa R: a dieta de Boldrini

Em uma postagem anterior, foi proposta a resolução de um problema de dietética usando o Maxima como instrumento. Eis uma proposta para resolução valendo-se do programa R. O programa R cada vez mais torna-se uma importante ferramenta, uma vez que é gratuito, de código aberto, e largamente usado para análises estatísticas com alto grau de confiabilidade.


# Vetores por alimento
alimentoI <- c(1,10,1,2,2)
alimentoII <- c(9,1,0,1,1)
alimentoIII <- c(2,2,5,1,2)
alimentoIV <- c(1,1,1,2,13)
alimentoV <- c(1,1,1,9,2)

# Matriz de coeficientes
A <- rbind(alimentoI,alimentoII,alimentoIII,alimentoIV,alimentoV)

# Vetor/matriz de termos independentes
b <- c(170,180,140,180,350)
B <- cbind(b)

# Transpondo A (cada linha deve totalizar a quantidade de um mesmo nutriente)
At <- t(A)

# Calculando Matriz de incógnitas (AX=B)
Vitaminas <- paste("Vitamina",LETTERS[1:5],sep=" ")
X <- data.frame(solve(At,B),row.names=Vitaminas)
X

# A fórmula para cálculo é apenas "solve(At,B)"
# O restante do código serve apenas
# para apresentar melhor os resultados,
# evitando-se confusão entre vitaminas
# de A a E, e alimentos de I a V

Baixar código para R

Nenhum comentário:

Postar um comentário

O saber só faz sentido se for coletivo. Construamos juntos!