Você está aqui: start » cursos » rbelem » dia2

Diferenças

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

Link para esta página de comparações

cursos:rbelem:dia2 [2010/05/26 22:55] – criada paulojuscursos: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>
  

QR Code
QR Code cursos:rbelem:dia2 (generated for current page)