disciplinas:ce223:comandos2008

Não foi possível enviar o arquivo. Será algum problema com as permissões?

Você está aqui: start » disciplinas » ce223 » comandos2008

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- disciplinas:ce223:comandos2008 [2008/02/28 13:29]
paulojus
+++ disciplinas:ce223:comandos2008 [2008/03/17 16:53]
ehlers
@@ Linha 23: / Linha 23: @@
 seq(0,1, by=0.1)
 (0:10)/10
+*(0:10)
+seq(0,20,by=2)
 :0
 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//
@@ Linha 37: / Linha 40: @@
 x
 </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.
@@ Linha 43: / Linha 55: @@
 x
-pesos <- c(67, 83, 56, 91, 58, 47, 82, 75)
 pesos
+pesos[4]
 pesos[4] <- 81
 pesos
+pesos[pesos > 80]
 pesos[pesos > 80] <- 85
 pesos
 pesos[pesos > 80] <- NA
+pesos
+pesos[is.na(pesos)] <- 90
 pesos
 </code>
-Amostrando valores de um vetor.
+Identificando as posições dos elementos que satisfazem certo critério.
+<code R>
+pesos
+which(pesos == 56)
+which(pesos == 90)
+which(pesos < 70)
+</code>
+Amostrando valores de um vetor. Note uso do ''rep'' para definir amostra com reposição. ''args()'' mostra os argumentos da função.
 <code R>
 dados <- c(34, 28, 31, 32, 43, 40, 45, 39, 26, 29)
 sample(dados, 3)
 sample(dados, 3)
-sample(dados, 3, rep=TRUE) # com reposição
+sample(dados, 3, rep=TRUE)
-<code>
+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>
+* 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     ^
+Comentários sobre ordem de entrada dos dados, cliclagem das variáveis e definição das dimensões do array<code R>
+freqs <- scan()
+: 45
+: 28
+: 37
+: 16
+: 22
+: 15
+: 21
+: 34
+: 56
+: 33
+: 21
+: 30
+: 40
+: 50
+: 85
+: 45
+: 37
+: 29
+:
+freqs
+Af <- array(freqs, dim=c(3,2,3))
+Af
+</code>
+==== Semana 4 ====
+=== 17/03/2008 ===
+<code R>
+freqs = scan(file='http://leg.ufpr.br/~ehlers/CE223/fumo.dat')
+array(freqs, dim=c(3,2,3))
+nomes = list(c('PR','SC','RS'), c('M','F'), c('nao fuma','fuma pouco','fuma muito'))
+hf = array(freqs, dim=c(3,2,3), dimnames=nomes)
+hf
+m1 <- matrix(1:12, ncol = 3)
+m1
+dimnames(m1)
+dimnames(m1) <- list(c("L1", "L2", "L3", "L4"), c("C1", "C2", "C3"))
+m1
+m2 <- cbind(1:5, 6:10)
+m2
+m3 <- cbind(1:5, 6)
+m3
+margin.table(m1, margin = 1)
+apply(m1, 1, sum)
+rowSums(m1)
+margin.table(m1, margin = 2)
+apply(m1, 2, sum)
+colSums(m1)
+#operacoes com matrizes
+m4 <- matrix(1:6, nc = 3)
+m5 <- matrix(10 * (1:6), nc = 3)
+m4
+m5
+m4 + m5
+m4 * m5
+m5 - m4
+m5/m4
+m4 %*% m5
+t(m4)
+m6 = t(m4)%*% m5
+solve(m6)
+m6[3,3]=20
+solve(m6)
+mat <- matrix(c(1, 5, 2, 3, -2, 1, -1, 1, -1), nc = 3)
+vec <- c(10, 15, 7)
+solve(mat, vec)
+</code>

Diferenças

Navegação

Busca

Ferramentas

QR Code