Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior Próxima revisão Ambos lados da revisão seguinte | ||
disciplinas:ce223:comandos2008 [2008/02/28 13:38] paulojus |
disciplinas:ce223:comandos2008 [2008/03/13 23:52] paulojus |
||
---|---|---|---|
Linha 23: | Linha 23: | ||
seq(0,1, by=0.1) | seq(0,1, by=0.1) | ||
(0:10)/10 | (0:10)/10 | ||
+ | |||
+ | 2*(0:10) | ||
+ | seq(0,20,by=2) | ||
10:0 | 10:0 | ||
seq(10,0, by=-1) | seq(10,0, by=-1) | ||
- | <code> | + | </code> |
Selecionando indivíduos pela sua posição (indexando valores de um vetor). Note que comandos a seguir somente mostram resultados sem alterar //x// | Selecionando indivíduos pela sua posição (indexando valores de um vetor). Note que comandos a seguir somente mostram resultados sem alterar //x// | ||
Linha 37: | Linha 40: | ||
x | x | ||
</code> | </code> | ||
+ | |||
+ | Selecionando valores do vetor segundo outros critérios. Nos comandos a seguir o objeto original não é alterado. | ||
+ | <code R> | ||
+ | pesos <- c(67, 83, 56, 91, 58, 47, 82, 75) | ||
+ | pesos[pesos > 80] | ||
+ | pesos[pesos < 50 | pesos > 80] | ||
+ | pesos[pesos > 50 & pesos < 80] | ||
+ | </code> | ||
+ | |||
Substituindo valores de um vetor. Note que comandos alteram valores do vetor. | Substituindo valores de um vetor. Note que comandos alteram valores do vetor. | ||
Linha 43: | Linha 55: | ||
x | x | ||
- | pesos <- c(67, 83, 56, 91, 58, 47, 82, 75) | ||
pesos | pesos | ||
+ | pesos[4] | ||
pesos[4] <- 81 | pesos[4] <- 81 | ||
pesos | pesos | ||
+ | pesos[pesos > 80] | ||
pesos[pesos > 80] <- 85 | pesos[pesos > 80] <- 85 | ||
pesos | pesos | ||
pesos[pesos > 80] <- NA | pesos[pesos > 80] <- NA | ||
+ | pesos | ||
+ | pesos[is.na(pesos)] <- 90 | ||
pesos | pesos | ||
+ | </code> | ||
+ | |||
+ | Identificando as posições dos elementos que satisfazem certo critério. | ||
+ | <code R> | ||
+ | pesos | ||
+ | which(pesos == 56) | ||
+ | which(pesos == 90) | ||
+ | which(pesos < 70) | ||
</code> | </code> | ||
Linha 60: | Linha 83: | ||
sample(dados, 3, rep=TRUE) | sample(dados, 3, rep=TRUE) | ||
args(sample) | args(sample) | ||
+ | </code> | ||
+ | |||
+ | Ordenando valores | ||
+ | <code R> | ||
+ | dados | ||
+ | sort(dados) | ||
+ | rev(sort(dados)) | ||
+ | sort(dados, decreasing=T) | ||
+ | </code> | ||
+ | |||
+ | Acrescentando elementos e concatenando dois ou mais vetores<code R> | ||
+ | x <- 0:10 | ||
+ | x | ||
+ | x <- c(x, 11) | ||
+ | x | ||
+ | |||
+ | x1 <- c(51:60, 101:110) | ||
+ | x1 | ||
+ | |||
+ | x2 <- c(x, c(15, 18, 21, 15, 30)) | ||
+ | x2 | ||
+ | |||
+ | x3 <- c(x2, x1) | ||
+ | x3 | ||
+ | |||
+ | c(1:10, seq(21, 50, by=2), c(100, 110, 150)) | ||
+ | </code> | ||
+ | |||
+ | ==== Semana 2 ==== | ||
+ | === 03/03/2008 === | ||
+ | |||
+ | Operadores aritméticas e funções para operações aritméticas, prioridade de operações, uso de parêntesis. | ||
+ | <code R> | ||
+ | x <- 0:10 | ||
+ | px <- choose(10, x) * (0.3)^x * (1-0.3)^(10-x) | ||
+ | px | ||
+ | </code> | ||
+ | |||
+ | Gráficos | ||
+ | <code R> | ||
+ | plot(x, px) # ver gráfico | ||
+ | plot(x, px, type="h") # ver gráfico | ||
+ | </code> | ||
+ | |||
+ | Funções relacionadas a distribuições de probabilidades | ||
+ | <code R> | ||
+ | px <- dbinom(x, 10, 0.3) | ||
+ | px | ||
+ | </code> | ||
+ | |||
+ | Gráfico da função acumulada | ||
+ | <code R> | ||
+ | Fx <- pbinom(x, 10, 0.3) | ||
+ | Fx | ||
+ | cumsum(px) # note que soma acumulada de dbinom() fornece mesmo resultado que pbinom() | ||
+ | plot(x, Fx, type="s") | ||
+ | </code> | ||
+ | |||
+ | Lei da reciclagem (//recycling rule//) -- válida para as operações aritméticas | ||
+ | <code R> | ||
+ | 2 * x | ||
+ | (0.3)^x | ||
+ | x1 <- c(2, 5, 7, 8) | ||
+ | x1 | ||
+ | x2 <- c(10, 20) | ||
+ | x2 | ||
+ | x1 + x2 | ||
+ | x2 + x1 | ||
+ | x3 <- (1:3)*10 | ||
+ | x3 | ||
+ | x1 + x3 | ||
+ | x3 + x1 | ||
+ | x2 + x3 | ||
+ | x3 + x2 | ||
+ | x1 * x2 | ||
+ | x1 * x2 | ||
+ | exp(x1) + log(x2) | ||
+ | </code> | ||
+ | |||
+ | Argumentos de funções: ordem dos argumentos, nomes dos argumentos, casamento parcial de nomes, argumentos com //default// | ||
+ | <code R> | ||
+ | args(dbinom) | ||
+ | dbinom(x, 10, 0.3) | ||
+ | dbinom(x=x, size=10, prob=0.3) | ||
+ | dbinom(prob=0.3, x=x, size=10) | ||
+ | dbinom(size=10, prob=0.3, x=x) | ||
+ | dbinom(x, s=10, p=0.3) | ||
+ | dbinom(x, p=0.3, s=10) | ||
+ | dbinom(x, p=0.3, s=10) | ||
+ | dbinom(x, 10, 0.3, log=F) | ||
+ | dbinom(x, 10, 0.3, log=T) | ||
+ | </code> | ||
+ | |||
+ | === 05/03/2008 === | ||
+ | |||
+ | Verificando máximos, mínimos e suas posições em um vetor.<code R> | ||
+ | pesos <- c(67, 83, 56, 91, 58, 47, 82, 75) | ||
+ | max(pesos) | ||
+ | min(pesos) | ||
+ | which(pesos == max(pesos)) | ||
+ | which.max(pesos) | ||
+ | which(pesos == min(pesos)) | ||
+ | which.min(pesos) | ||
+ | </code> | ||
+ | |||
+ | Testando pela ocorrência de ''NA''. Note o uso do caracter de negação ''!'' <code R> | ||
+ | dat <- c(43, 56, NA, 23, 48, 33, NA, 29, 33, 39) | ||
+ | is.na(dat) | ||
+ | !is.na(dat) | ||
+ | which(is.na(dat)) | ||
+ | which(!is.na(dat)) | ||
+ | dat1 <- dat[!is.na(dat)] | ||
+ | </code> | ||
+ | |||
+ | ==== Semana 3 ==== | ||
+ | === 10/03/2008 === | ||
+ | |||
+ | Outras funções de probabilidades<code R> | ||
+ | args(dchisq) | ||
+ | args(dt) | ||
+ | </code> | ||
+ | |||
+ | Algumas contantes e operadores<code R> | ||
+ | exp(1) | ||
+ | exp(3) | ||
+ | print(pi, dig=12) | ||
+ | print(pi, dig=12) | ||
+ | pi | ||
+ | options(digits=12) | ||
+ | pi | ||
+ | exp(1) | ||
+ | </code> | ||
+ | |||
+ | Examinar opções de ''options()''<code R> | ||
+ | options() | ||
+ | </code> | ||
+ | |||
+ | |||
+ | === 12/03/2008 === | ||
+ | |||
+ | Funções para verificar exitência de ''NA'''s ''NaN'''s e valores infinitos (''Inf'') <code R> | ||
+ | x <- c(5, 0, -2) | ||
+ | x/0 | ||
+ | is.nan(x/0) | ||
+ | is.finite(x/0) | ||
+ | !is.finite(x/0) | ||
+ | </code> | ||
+ | |||
+ | Explicações sobre como o R armazana objetos (RAM e/ou dispositivos como por exemplo o HD)<code R> | ||
+ | save.image() | ||
+ | q() | ||
+ | </code> | ||
+ | |||
+ | Listando e apagando objetos | ||
+ | <code R> | ||
+ | ls() | ||
+ | objects() | ||
+ | rm(x) | ||
+ | rm(list=ls()) | ||
+ | args(ls) | ||
+ | ls() | ||
+ | ls(all=T) | ||
+ | </code> | ||
+ | |||
+ | Tipos de objetos no R: matrizes | ||
+ | <code R> | ||
+ | x <- 1:24 | ||
+ | x | ||
+ | m <- matrix(x, nr=6) | ||
+ | m | ||
+ | m <- matrix(x, nr=6, byrow=T) | ||
+ | m | ||
+ | x <- 1:25 | ||
+ | m <- matrix(x, nr=6) | ||
+ | m | ||
+ | attributes(x) | ||
+ | attributes(m) | ||
+ | </code> | ||
+ | |||
+ | Tipos de objetos no R: arrays | ||
+ | <code R> | ||
+ | x <- 1:24 | ||
+ | a <- array(x, dim=c(4,3,2)) | ||
+ | a | ||
+ | </code> | ||
+ | |||
+ | Digitação e conversão de uma tabela de tripla entrada (dada no quadro durante a aula) em um objeto do tipo ''array'' | ||
+ | | ^ PR ^^ SC ^^ RS ^^ | ||
+ | | ^ Masculino ^ Feminino ^ Masculino ^ Feminino ^ Masculino ^ Feminino ^ | ||
+ | |Não Fuma ^ 45 ^ 16 ^ 21 ^ 33 ^ 40 ^ 45 ^ | ||
+ | |Fuma pouco ^ 28 ^ 22 ^ 34 ^ 21 ^ 50 ^ 37 ^ | ||
+ | |Fuma muito ^ 37 ^ 15 ^ 56 ^ 30 ^ 85 ^ 29 ^ | ||
+ | |||
+ | <code R> | ||
+ | freqs <- scan() | ||
+ | freqs | ||
+ | Af <- array(freqs, dim=c(3,2,3)) | ||
+ | Af | ||
</code> | </code> | ||