Setembro 2011 - R-br

[NEWS] Mudanças nas versões do R e melhora de performance
by Fernando Mayer 16 Fev '12

16 Fev '12

Não sei se todos já estão sabendo, mas li isso hoje e achei importante repassar para a lista: Sobre as mudanças de releases do R: R 2.14 to be released on October 31; R 2.13 patch on September 13 http://blog.revolutionanalytics.com/2011/09/r-214-to-be-released-on-october… Brian Ripley on The R Development Process http://blog.revolutionanalytics.com/2011/08/brian-ripley-on-the-r-developme… Sobre a melhora de performance: With byte compiler, R 2.14 will be even faster http://blog.revolutionanalytics.com/2011/08/with-byte-compiler-r-214-will-b… []s, --- Fernando Mayer Dept. Ecologia e Zoologia - ECZ/CCB Universidade Federal de Santa Catarina - UFSC URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com

4 4

gerar dados
by Ricardo Pereira 04 Out '11

04 Out '11

Outra duvida. Estou tentando simular dados, tal como um exemplo em outro programa... Exemplo abaixo: proc iml; do i=1 to 1000; seed=456; var=150; w = RANNOR(seed)*SQRT(var); end; *tentei:* set.seed(456) var <- 150 for (i in 1:1000) { w [i] <- rnorm * sqrt(varg) } *ERRO:* Erro em W[j] <- 1 : objeto 'W' não encontrado

10 15

Mlp redes neurais
by Fernando Neto 03 Out '11

03 Out '11

Por favor, estou precisando usar a pacote que tenha a funcao mlp (multilayer perceptron) so q o pacote q achei como referencia esta dando Not Found. Alguem poderia me ajudar? Estou usando o R no ubuntu. Grato, Fernando

3 7

Função não reconhece variável extraída de um modelo!
by Ivan Bezerra Allaman 30 Set '11

30 Set '11

Boa tarde senhores! Gostaria de saber como posso fazer uma função ler uma variável a partir de um modelo, nls por exemplo. Vamos supor um CMR. da <- expand.grid( A=factor( LETTERS[1: 3]), x=1:10) da$y <- 10/(1+exp(-( da$x-5)/1) )+rnorm(da$ x,0,0.1) n0 <- nls(y~As[A]/ (1+exp(-( x-x0[A])/ S[A])), data=da, start=list(As= c(10,10,10) , x0=c(5,5,5), S=c(1,1,1))) fun <- function(modelo){ mean(all.vars(summary(modelo)$formula)[1]) } Percebam que eu quero que a função calcule a média de uma variável, cuja variável foi declarada num determinado modelo. A função acima não dá certo, pois a função entende o seguinte: mean("VARIÁVEL") Para dar certo, eu acho, a variável não poderia estar entre aspas. Enfim, alguém tem alguma idéia de como eu posso fazer isso? Desde já agradeço a todos. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 4

Duvida sobre if e ifelse no R
by Augusto Ribas 30 Set '11

30 Set '11

Ola a todos, tenho uma duvida aqui. Estou lendo um livro sobre programação em R e empaquei num problema aqui. Eu tenha a função: x ≤ 0 ∈ (0, 1] > 1 f(x) −x3 x2 √x ai resolvo no intervalo de -2 a 2 com o comando if: #começo do script rm(list=ls()) # input x.values <- seq(-2, 2, by = 0.1) # for each x calculate y n <- length(x.values) y.values <- rep(0, n) for (i in 1:n) { x<-x.values[i] # your expression for y goes here if (x>1) { y<-sqrt(x) } else { if (x<=0) { y<- -x^3 } else { y<- x^2 } } #end of y expression y.values[i] <- y } # output plot(x.values, y.values, type = "l") #até aqui tudo funcionou beleza. #mas dai fui tentar fazer com ifelse a mesma coisa y.values2 <- rep(0, n) resultado<-ifelse(x.values>1,y.values2<-sqrt(x.values), ifelse(x.values<=0,y.values2<- -x.values^3, y.values2<- x.values^2) ) #o vetor resultado fica correto resultado==y.values #mas o vetor y.values2 fica errado y.values2==y.values #o resultado ficou [1] 8.000000 6.859000 5.832000 4.913000 4.096000 3.375000 2.744000 2.197000 [9] 1.728000 1.331000 1.000000 0.729000 0.512000 0.343000 0.216000 0.125000 [17] 0.064000 0.027000 0.008000 0.001000 0.000000 0.010000 0.040000 0.090000 [25] 0.160000 0.250000 0.360000 0.490000 0.640000 0.810000 1.000000 1.048809 1.303840 1.341641 1.378405 [41] 1.414214 #mas os valores que ficam salvos são no vetor y.values2 são: > y.values2 [1] 4.00 3.61 3.24 2.89 2.56 2.25 1.96 1.69 1.44 1.21 1.00 0.81 0.64 0.49 0.36 [16] 0.25 0.16 0.09 0.04 0.01 0.00 0.01 0.04 0.09 0.16 0.25 0.36 0.49 0.64 0.81 [31] 1.00 1.21 1.44 1.69 1.96 2.25 2.56 2.89 3.24 3.61 4.00 #minha duvida então é pq ele so resolve tudo elevando ao quadrado no y.values2 #o que acontece, por exemplo, com o primeiro valor de resultado que vejo 8 qd mas fica salvo 4 quando vou olhar o y.values2. #alguém pode me iluminar sobre o que esta acontecendo? -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

4 5

Distância entre strings
by Pedro Rafael 29 Set '11

29 Set '11

Pessoal qual o objeto para compara duas strings e observar a distância entre elas? Exemplo: Gostaria que Pedro Rafael fosse igual à peDro Rafael, que também fosse igual à Rafael Pedro. Existe alguma forma de fazer isto? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

7 18

DIFFTIME
by Edson Lira 29 Set '11

29 Set '11

Gente, estou reenviando este email, tenho os dados abaixo: DIA SAIDA CHEGADA 1 09/09/2011 08:11 11:20 2 09/09/2011 08:17 08:25 3 09/09/2011 08:45 09:05 4 09/09/2011 09:15 09:20 5 09/09/2011 09:30 09:35 6 09/09/2011 09:50 10:15 7 09/09/2011 10:40 10:55 8 12/09/2011 08:25 11:05 9 12/09/2011 08:30 08:45 10 12/09/2011 09:25 09:40 Todas as variáveis tem a estrutura de fatores. Vejam que as datas (DIA), estão separadas da hora da saida (SAIDA) e da hora da chegada (CHEGADA). Preciso calcular a diferença em minutos entre a hora da chegada e a hora da saida (CHEGADA-SAIDA). O Augusto me respondeu sugerindo usar o help para: ?as.POSIXct ?strptime ?difftime Verifiquei, mas não consegui, sempre tenho dificuldades em trabalhar com data. O pacote chron tem conversões, calcula diferenças, mas não consegui implementar seus comandos para o meu problema acima. []'s. Edson Lira Estatístico Manaus-Amazonas

2 1

conversão coordenadas decimais em cartesianas
by Heloíse Pavanato 29 Set '11

29 Set '11

Olá, Alguém sabe se existe um pacote para converter coordenadas geográficas em cartesianas?? Obrigada, Heloise

2 2

Problema com matrizes e vetores
by Thiago Touya 28 Set '11

28 Set '11

Senhores, preciso estimar um vetor de parâmetros (Beta) usando as saídas do svd() para comparar com a estimação pelo método dos mínimos quadrados. O grande problema é que nao sei isolar Beta(3x1) na seguinte equação: H'XBeta = H'Y = lambda'Beta, onde lambda'(1x3) Eu tenho H'(1x3), X(3x3) e Y(3x1), mas nao estou visualizando uma maneira de isolar beta, pois H'X é um vetor. X <- matrix(c(3,0,-3,0,6,3,-3,3,6),ncol=3,nrow=3) dvs <- svd(X) lambda <- (dvs$d)^2 Y <- c(5,2,3) D <- matrix(c(lambdas[1,1],0,0,0,lambdas[2,1],0,0,0,lambdas[3,1]), ncol=3, nrow=3) C <- dvs$u Vlinha <- t(dvs$v) sabe-se que pela Decomposição em valores singulares X=CDV' (neste caso em que X é simétrica e nao temos nenhum autovalor igual a zero) Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

1 0

Dúvida simples
by Victor Eduardo 28 Set '11

28 Set '11

Gente, como posso colocar o mi linha no R com o comando expression ??? Abraços!

3 4