#========================================================================================== # ce083 Estatistica Computacional I - Profs Walmes Zeviani & Elias Krainski (LEG/DEST/UFPR) # Aula 6 - Operacoes aritmeticas, Operacoes com matrizes, funcoes (29/09/2011) # Elias T Krainski #========================================================================================== ### operações aritmeticas 2+2 2+2*3+2^2 2 + 3*2^2 2^2 ### igual a 2**2 a <- c(2, 3, 4, 1, -2, 0) a+10 a*3 2^a a/3 3/a b <- c(-3, 0, 2) b/0 c <- 2 c+3 c(2, 3) ### funções prontas no R sqrt(9) sin(3) pi sin(pi) format(pi, dig=19) print(pi, dig=19) factorial(5) choose(5, 2) combn(5, 2) d <- c(2, 3, NA, 4) d a mean(a) mean(d) mean(d, na.rm=TRUE) var(d, na.rm=TRUE) sd(a) quantile(a) quantile(a, c(0.3)) quantile(a, c(0.3, .60, .9)) min(a) max(a) range(a) diff(c(2, 5)) diff(range(a)) a diff(a) summary(a) d is.na(d) !is.na(d) ### quantos NAs tem em d sum(is.na(d)) is.finite(2) is.finite(Inf) a 2/a is.finite(2/a) summary(a) ### operações com matrizes A <- matrix(c(2, 1, 3, 1), 2) B <- c(1, 3) D <- matrix(c(2,3, 1,2, -1,2), 2) E <- matrix(c(1,2, -1,-2), 2) A + E t(A) A * E A %*% E A %*% B A%*%D det(A) diag(A) sum(diag(A)) solve(A, c(2,3)) solve(A) ### inversa de A diag(2) solve(A, diag(2)) ### inversa de A iA <- solve(A) ### inversa de A A%*%iA diag(c(2, 3, 1)) ### criando funções meusumario <- function(x) { resultado <- c(Mínimo=min(x), Mediana=median(x), Média=mean(x), Máximo=max(x)) return(resultado) } meusumario(a) meusumario(1:100) ### funcao para cálculo do IMC imc <- function(altura, massa) massa/(altura^2) imc(1.8, 80) alt <- c(1.8, 1.75) imc(alt, 80) imc(1.8, c(70, 80, 90))