# 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!