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.
|
cursos:rbelem:dia2 [2010/05/26 22:55] paulojus criada |
cursos:rbelem:dia2 [2010/05/29 11:12] (atual) paulojus |
||
|---|---|---|---|
| Linha 64: | Linha 64: | ||
| ## ex de reg e despacho de métodos | ## ex de reg e despacho de métodos | ||
| xydata <- data.frame(x = 1:20, | xydata <- data.frame(x = 1:20, | ||
| - | y = c(0.9,5.1,2.4,8.1,4.2,7.1,5.6,7.6,5.9,7.7,11.8,6.9,9.3,10.9,8.4,11.6,13.0,13.8,13.1,9.3) | + | y = c(0.9,5.1,2.4,8.1,4.2,7.1,5.6,7.6,5.9,7.7, |
| + | 11.8,6.9,9.3,10.9,8.4,11.6,13.0,13.8,13.1,9.3) | ||
| ) | ) | ||
| reg <- lm(y ~ x, data = xydata) | reg <- lm(y ~ x, data = xydata) | ||
| Linha 96: | Linha 97: | ||
| with(quest, boxplot(Idade ~ Toler)) | with(quest, boxplot(Idade ~ Toler)) | ||
| - | ## fazendo gráfico de função | + | ########################################################################## |
| + | ## fazendo gráfico de uma função | ||
| x <- seq(0, 60, length=500) | x <- seq(0, 60, length=500) | ||
| y <- 3 + 15 * exp(-x/10) | y <- 3 + 15 * exp(-x/10) | ||
| Linha 106: | Linha 108: | ||
| class(fx) | class(fx) | ||
| - | ## definind um parâmetro que pode ser alterado | + | ## definindo um parâmetro da função que pode ser alterado |
| fx <- function(x, alpha) 3 + 15 * exp(-x/alpha) | fx <- function(x, alpha) 3 + 15 * exp(-x/alpha) | ||
| - | plot(x, fx(x, alpha=10), ty="l") | + | plot(x, fx(x, alpha=10), ty="l", ylab="f(x)") |
| lines(x, fx(x, alpha=5), ty="l", col=2) | lines(x, fx(x, alpha=5), ty="l", col=2) | ||
| lines(x, fx(x, alpha=15), ty="l", col=4) | lines(x, fx(x, alpha=15), ty="l", col=4) | ||
| - | legend("topright", c("alpha=5", "alpha=10", "alpha=15"), col=c(2,1,4), lty=1) | + | legend("topright", c(expression(alpha==5), |
| - | title(expression(alpha)) | + | expression(alpha==10), |
| + | expression(alpha==15)), col=c(2,1,4), lty=1) | ||
| + | title(expression(paste("Função f(x) com diferentes valores de ",alpha))) | ||
| pdf("exponenciais.pdf") | pdf("exponenciais.pdf") | ||
| - | plot(x, fx(x, alpha=10), ty="l") | + | fx <- function(x, alpha) 3 + 15 * exp(-x/alpha) |
| + | plot(x, fx(x, alpha=10), ty="l", ylab="f(x)") | ||
| lines(x, fx(x, alpha=5), ty="l", col=2) | lines(x, fx(x, alpha=5), ty="l", col=2) | ||
| lines(x, fx(x, alpha=15), ty="l", col=4) | lines(x, fx(x, alpha=15), ty="l", col=4) | ||
| - | legend("topright", c("alpha=5", "alpha=10", "alpha=15"), col=c(2,1,4), lty=1) | + | |
| - | title(expression(alpha)) | + | legend("topright", c(expression(alpha==5), |
| + | expression(alpha==10), | ||
| + | expression(alpha==15)), col=c(2,1,4), lty=1) | ||
| + | title(expression(paste("Função f(x) com diferentes valores de ",alpha))) | ||
| dev.off() | dev.off() | ||
| Linha 126: | Linha 134: | ||
| fx <- function(x) 3 + 15 * exp(-x/10) | fx <- function(x) 3 + 15 * exp(-x/10) | ||
| integrate(fx, 10, 20) | integrate(fx, 10, 20) | ||
| - | integrate(fx) | + | integrate(fx) # tem, obrigatoriamente, que especificar limites de integração!!! |
| integrate(fx, 0, 10) | integrate(fx, 0, 10) | ||
| - | ## passando argumento adicional a função | + | ## passando argumento adicional para a função |
| fx <- function(x, alpha) 3 + 15 * exp(-x/alpha) | fx <- function(x, alpha) 3 + 15 * exp(-x/alpha) | ||
| integrate(fx, 10, 20, alpha=5) | integrate(fx, 10, 20, alpha=5) | ||
| Linha 135: | Linha 143: | ||
| ## nem tudo é integrável... | ## nem tudo é integrável... | ||
| integrate(fx, 10, Inf, alpha=5) | integrate(fx, 10, Inf, alpha=5) | ||
| - | Erro em integrate(fx, 10, Inf, alpha = 5) : | + | #Erro em integrate(fx, 10, Inf, alpha = 5) : |
| - | the integral is probably divergent | + | # the integral is probably divergent |
| ## definindo uma outra função | ## definindo uma outra função | ||
| Linha 142: | Linha 150: | ||
| ## .. é verificando que é uma f.d.p. | ## .. é verificando que é uma f.d.p. | ||
| integrate(fx, 0, Inf) | integrate(fx, 0, Inf) | ||
| - | 1 with absolute error < 0.00011 | + | #1 with absolute error < 0.00011 |
| ## a função só é definida para valores não negativos... | ## a função só é definida para valores não negativos... | ||
| integrate(fx, -Inf, Inf) | integrate(fx, -Inf, Inf) | ||
| - | Erro em integrate(fx, -Inf, Inf) : non-finite function value | + | #Erro em integrate(fx, -Inf, Inf) : non-finite function value |
| ## e a função fica melhor definida especificando seu domínio | ## e a função fica melhor definida especificando seu domínio | ||
| fx <- function(x) {ifelse(x>0, (1/10) * exp(-x/10), 0)} | fx <- function(x) {ifelse(x>0, (1/10) * exp(-x/10), 0)} | ||
| integrate(fx, -Inf, Inf) | integrate(fx, -Inf, Inf) | ||
| - | 1 with absolute error < 0.00011 | + | #1 with absolute error < 0.00011 |
| - | ## para o exemplo no materia use a função abs() para definir módulo | + | ## para o exemplo no material use a função abs() para definir módulo |
| abs(-5) | abs(-5) | ||
| </code> | </code> | ||
