Junho 2012 - R-br - DINF UFPR

Duvida envelope simulado para ZINB
by Humberto Hazin 07 Jun '12

07 Jun '12

Olá Pessoal, Estou precisando fazer uma análiser de resíduo para um ZINB e gostaria de aplicar um envelope simulado! existe algum pacote ou codico que faça? Humberto

2 2

[OFF-TOPIC] Free Webinar: Spatial Analysis in R
by ASANTOS 07 Jun '12

07 Jun '12

Bom tarde Pessoal, Para quem trabalha com estatística espacial ou tem interesse sobre o assunto, encontrei um webinar grátis interessante na próxima semana (13/06 as 14:00h), intitulado: Finding Meaning in Points, Areas and Surfaces: Spatial Analysis in R com o professor David Unwin, segue o link para inscrição: http://www.revolutionanalytics.com/news-events/free-webinars/2012/spatial-a… Abraço, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

1 0

Re: [R-br] Help me please! ^^
by Ivan Bezerra Allaman 07 Jun '12

07 Jun '12

Irei partir do princípios que as medidas repetidas são independentes. Se o interesse é comparar mais de duas médias podemos usar a técnica da análise de variância seguida de um teste de comparação múltipla. Segue o script. DIA1 <- c(7.306,7.296,7.326,7.334,7.334) DIA2 <- c(7.244,7.258,7.288,7.304,7.240) DIA3 <- c(6.998,7.044,7.044,7.076,7.086) DIA4 <- c(6.940,6.956,6.98,6.998,7.002) DIA5 <- c(7.116,7.130,7.142,7.148,7.150) pot <- c(DIA1,DIA2,DIA3,DIA4,DIA5) dias <- rep(c('DIA1','DIA2','DIA3','DIA4','DIA5'),rep(5,5)) dad <- data.frame(dias,pot) dad mod <- aov(pot ~ dias,data=dad) summary(mod) plot(mod)#Avaliando os pressupostos. library(agricolae) HSD.test(mod,'dias') (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Help me please! ^^
by Patricia alcantara nguyen 07 Jun '12

07 Jun '12

Olá, Tudo bem? Preciso fazer um teste de repetitividade da potência de uma balança, mas estou tendo dificuldades com o R. Uma balança de radiação foi usada para medir a potência de um equipamento de ultrassom (3 MHz) durante 5 dias. Em cada dia foram feitas 5 medições de uma determinada potência. Como implementar isso no R??? Quero verificar se a média da potência dos 5 dias são iguais. Potência Nominal - 7 W MEDIDA 1 MEDIDA 2 MEDIDA 3 MEDIDA 4 MEDIDA 5 DIA 1 7,306 7,296 7,326 7,334 7,334 DIA 2 7,244 7,258 7,288 7,304 7,240 DIA 3 6,998 7,044 7,044 7,076 7,086 DIA 4 6,940 6,956 6,98 6,998 7,002 DIA 5 7,116 7,130 7,142 7,148 7,150 Muito Obrigada! Patrícia Nguyen

2 1

plotar valores específicos com função image( )
by Vitor Aguiar 06 Jun '12

06 Jun '12

Olá, desculpe por estar postando a pergunta novamente, mas não tive resposta da primeira vez. Se alguém puder ajudar, agradeço muito. Tenho uma matrix com que estou fazendo um heatmap usando a função image( ). Gostaria de plotar não apenas os valores observados na matrix, mas também valores que seriam possíveis mas não são observados. Isso significa o seguinte: Tenho uma matrix 16x16. Valores possíveis, mas não observados, correspondem a "0" (zeros) onde a soma (índice da linha + índice da coluna) são iguais ou menores que 15. Exemplo: Se eu tenho um valor "0" na posição [14, 3], eu NÃO quero plotar, pois 14 + 3 = 17. Se eu tenho um valor "0" na posição [11, 3], eu quero plotar, pois 11 + 3 = 14. Eu queria que esses valores fossem plotados em cor diferente do meu heatmap, por exemplo, em azul. É possível fazer isso? Matrix: x = c(3045, 893, 692, 830, 617, 155, 246, 657, 105, 60, 18, 7, 7, 4, 2, 11234, 2985, 2242, 2471, 1575, 366, 503, 1283, 170, 79, 32, 6, 4, 1, 3, 19475, 4756, 3233, 3251, 1810, 409, 575, 1210, 139, 41, 11, 4, 2, 0, 0, 20830, 4739, 2990, 2531, 1346, 298, 325, 612, 60, 17, 1, 0, 1, 0, 0, 15304, 3196, 1885, 1440, 610, 117, 115, 185, 14, 2, 0, 0, 0, 0, 0, 8026, 1535, 806, 539, 223, 33, 37, 39, 0, 0, 0, 0, 0, 0, 0, 3300, 562, 286, 141, 45, 14, 5, 12, 0, 0, 0, 0, 0, 0, 0, 1067, 160, 65, 40, 14, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 277, 47, 6, 2, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 72, 6, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) xmat = matrix(x, ncol = 12) xmat = cbind(xmat, rep(0, 15), rep(0, 15), rep(0, 15), rep(0, 15)) xmat = rbind(xmat, rep(0, 16)) dimnames(xmat) = list(0:15, 0:15) image(log(xmat)) Muito Obrigado, Vitor

2 1

Re: [R-br] Alteração em função de pacote (novo tópico)
by Junior Beleti 06 Jun '12

06 Jun '12

Olá Benilton.. Com as informações que você apresentou sobre as execuções, realizei a execução do script novamente. Passando como parâmetro primeiramente (num.simu=1000), a qual ficou executando por mais de 18 horas e não conseguiu chegar ao final. Na segunda execução, utilizei (num.simu=100) e obtive o seguinte resultado, quanto ao tempo: usuário sistema decorrido 60869.136 76.144 61005.444 A chamada a função foi a seguinte, modificando apenas o valor do parâmetro: system.time(preditos.simu <- volta.cokri(md.cov.ck ,num.simu=100,int.conf=0.95)) Gostaria de saber como você executou para ter esses tempos tão pequenos. Att, Junior Beleti.

2 1

Escrever diretorio formato windows em arquivo
by Leandro Marino 06 Jun '12

06 Jun '12

Caros, estou preparando uns arquivos para uma simulação é preciso que o R guarde um vetor da seguinte forma x <- 't:\pasta' Isto depois será usado por um outro programa. No entanto, recebo a mensagem de erro Erro: '\p' é uma seqüência de escape não reconhecida na cadeia de caracteres começando com "T:\p" Alguma sugestão? Abraços, Leandro

5 10

list para data.frame
by Leandro Marino 05 Jun '12

05 Jun '12

Caros, estou com um grande dúvida, tenho a seguinte lista x <- list() x$s1 <- list() x$s2 <- list() x$s1$tp1 <- data.frame(v1=1:10,v2=3:12,s=1,tp=1) x$s1$tp2 <- data.frame(v1=1:10,v2=3:12,s=1,tp=2) x$s1$tp3 <- data.frame(v1=1:10,v2=3:12,s=1,tp=3) x$s2$tp1 <- data.frame(v1=1:10,v2=3:12,s=2,tp=1) x$s2$tp2 <- data.frame(v1=1:10,v2=3:12,s=2,tp=2) x$s2$tp3 <- data.frame(v1=1:10,v2=3:12,s=2,tp=3) quero juntar de tal forma que : s1 <- rbind(x[[1]][[1]],x[[1]][[2]],x[[1]][[3]]) s2 <- rbind(x[[1]][[1]],x[[1]][[2]],x[[1]][[3]]) no entanto, s varia de 1 a 300. isto é possivel?

2 2

Ajuda em modelo não linear
by Leonardo Guimarães de Oliveira 05 Jun '12

05 Jun '12

Como faço para encontrar o ponto na curva obtida na CMR abaixo onde a curva muda de direção. x=structure(list(Vt = c(0, 0.34, 3.68, 10.38, 21.94, 32.08, 48.54, 79.69, 94.82, 107.99), t = c(0L, 3L, 6L, 9L, 12L, 16L, 24L, 36L, 48L, 72L)), .Names = c("Vt", "t"), class = "data.frame", row.names = c(NA, -10L)) m4=nls(Vt~(107.99)/(1+exp(2-4*k*(t-g))), data=x, start=c(k=0.05,g=0.1)) m4 f4=function(x){(107.99)/(1+exp(2-4*0.02853*(x-8.71643)))} curve(f4, 1,100) Leonardo G. Oliveira Médico Veterinário Doutorando em Ciência Animal - EVZ/UFG

2 1

pergunta simples - diagonal de matrix
by Vitor Aguiar 04 Jun '12

04 Jun '12

Boa noite, se eu tenho uma matrix, por exemplo, 7x7 eu posso facilmente extrair a diagonal usando diag(mat). Mas e se eu quiser algo como "uma segunda diagonal", uma diagonal composta de mat[1,2] mat[2, 3] mat[3, 4] mat[4, 5] mat[5, 6] mat[6, 7]? Tem alguma forma simples de fazer? obrigado, Vitor

2 1