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