- R-br - DINF UFPR

Daniele deixou uma mensagem para você...
by Badoo 01 Out '11

01 Out '11

Daniele deixou uma mensagem para você... Só você pode ler o conteúdo desta mensagem e ver quem a enviou. Delete a qualquer momento ou responda imediatamente com o sistema de troca de mensagens. Para descobrir o que diz a mensagem, siga este link. http://us1.badoo.com/01132477499/in/Ng0pKMMcVZc/?lang_id=61&m=63&mid=4e8702… Mais gente que espera pacientemente por você: Ivan (Ilhéus, Brasil) Lyly Nogueira (Ilhéus, Brasil) Wagner (Curitiba, Brasil) http://us1.badoo.com/01132477499/in/Ng0pKMMcVZc/?lang_id=61&m=63&mid=4e8702… Se o link desta mensagem não funcionar, copie e cole-o em seu navegador. Este email faz parte de nosso sistema de entrega para a mensagem enviada por Daniele. Se recebeu este email por engano, por favor ignore-o. A mensagem será deletada em breve. Divirta-se! A Equipe Badoo Você recebeu este email, porque um membro do Badoo deixou uma mensagem para você no Badoo. Esta é uma mensagem de somente envio. Respostas a esta mensagem não são monitoradas ou respondidas. Se você não quer receber mais mensagens do Badoo, por favor nos notifique: http://us1.badoo.com/impersonation.phtml?lang_id=61&mail_code=63&email=r-br…

1 0

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 27 Set '11

27 Set '11

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

3 4

Conexão R com Oracle
by nani.design＠bol.com.br 27 Set '11

27 Set '11

2 1

Ponto de corte em cluster
by Marcelo Cardoso mello 27 Set '11

27 Set '11

Como obter no R um ponto de corte em um cluster (UPGMA ... etc). Algum critério como pseudo test F ou t ou algum outro? Obrigado. M

2 1