- R-br - DINF UFPR

Erro função sendmail - Windows
by Bruno Santos 09 Ago '11

09 Ago '11

Boa tarde, pessoal! Eu consegui utilizar normalmente a função sendmail do pacote sendmailR na minha máquina. > Sys.info() sysname release version nodename "Linux" "2.6.32-33-generic" "#71-Ubuntu SMP Wed Jul 20 17:30:40 UTC 2011" "bruno-laptop" machine login user "i686" "unknown" "bruno" Mas ao tentar usar o mesmo programa numa máquina Windows, ocorreu um erro que está depois do exemplo reproduzível. Seguem as informações da máquina Windows que eu fiz o teste. > Sys.info() sysname release version nodename machine login user "Windows" "7 x64" "build 7600" "RIGEL" "x86" "bramos" "bramos" Exemplo reproduzível, com argumento do email incorreto por segurança. ################################################################################# library(sendmailR) from <- sprintf("<sendmailR@%s>", Sys.info()[4]) to <- "<email@valido>" #Substituir por email válido subject <- "Hello from R" body <- list("It works!", mime_part(iris)) sendmail(from, to, subject, body, control=list(smtpServer="ASPMX.L.GOOGLE.COM")) ################################################################################# Erro em writeLines(msg, sock, sep = "\r\n") : argumento 'text' inválido Alguém tem alguma sugestão do que pode estar acontecendo? Obrigado, Bruno

5 5

Modelos não lineares....
by Thiago De paula protásio 09 Ago '11

09 Ago '11

Pessoal; Novamente recorro a ajuda de vcs...Estou tendo dificuldades no ajuste de um modelo não linear aos meus dados. Vejam mensagem de erro que está aparecendo. > a=6> b=0.16> c=-0.125 > plot(HT~DAP,data=dados)> curve(1.30+(x^2/a+b*x+c*x^2),add=T) > n <- nls(HT~(DAP^2/a+b*DAP+c*DAP^2)+1.3,data=dados,+ start=list(a=a, b=b, c=c))Erro em nlsModel(formula, mf, start, wts) : matriz gradiente singular nas estimativas iniciais dos parâmetros Em anexo segue a base de dados! Muito obrigado pela ajuda! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

mudança de nomes de colunas (variáveis) em uma lista
by Cassiano 09 Ago '11

09 Ago '11

Olá pessoal, sou iniciante no R e estou quebrando a cabeça em um problema que, acho, deve ter uma solução simples. Tenho uma lista de data.frames com os mesmos nomes nas colunas (nomes das variáveis). No total são 8 itens (data.frames). Quero mudar o nome de algumas variáveis e estou usando para isso, gsub(). Com um laço for consegui executar a tarefa. O código que usei foi for(i in 1:length(Indicador)) { names(Indicador[[i]]) <- gsub("Pontos", "Total.pontos.atend", names(Indicador[[i]])) names(Indicador[[i]]) <- gsub("Comercio", "Comercio", names(Indicador[[i]])) names(Indicador[[i]]) <- gsub("(.)*Folha", "Liberada.folha", names(Indicador[[i]])) names(Indicador[[i]]) <- gsub("(.)*Atendimento", "Ponto.atend", names(Indicador[[i]])) } onde, Indicador = list(), e as strings dentro do gsub são os nomes da variáveis originais e os nomes novos, respectivamente. Porém acredito que exista uma forma mais eficaz, provavelmente sem usar o laço for. Tentei de outra maneira (para o primeiro gsub do laço anterior) com > gsub("Pontos", "Total.pontos.atend", lapply(Indicador, names)) Esse último gsub fez as mudanças, porém gera uma saída como [1] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [2] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [3] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [4] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [5] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [6] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [7] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [8] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [9] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" Aqui, "Pontos" foi alterado para "Total.pontos.atend". Se isso estiver certo como faço para atribuir o novo nome para os nomes das variáveis dos data.frames em cada item da lista "Indicador". Abraços! Cassiano Novais Consultor economista

2 1

Barras de erro em barplot
by Marcelo Claro de Souza 08 Ago '11

08 Ago '11

Possuo uma planilha distribuída da seguinte forma: cidade pH desvpad..... Brasília 4,5 0,5 Itirapina 4,3 0,2 Como faço p plotar acrescentar o desvio padrão superior e inferios a minha figura partindo de barplot(pH~cidade,data=dados, .....) Muito obrigado.

3 2

Re: [R-br] Leitura arquivo rtf
by Roney Fraga Souza 08 Ago '11

08 Ago '11

Eder, Tentei abrir esse arquivo (ST-26-E.rtf) no Mac OS pelo TextEdit, Word, LibreOffice, Pages e Xcode, alguns apresentaram que o arquivo pode estar corrompido e outros não deram mensagem de erro, mas também não abrem. Só visualizei o conteúdo do documento pelo Komodo Edit, ou quando modifiquei a extensão para .txt, será que esse arquivo não esta corrompido? Att Roney

3 2

Re: [R-br] Barras de erro em barplot
by Ivan Bezerra Allaman 08 Ago '11

08 Ago '11

Ou, PH <- c(4.5, 4.3) error <- c(0.5,0.2) cidade <- c('Brasília','Itirapina') da <- data.frame(cidade,PH,error) library(ggplot2) limits <- aes(ymax=PH + error, ymin=PH - error) p <- ggplot(da, aes(cidade,PH)) p + geom_bar(fill='blue') + geom_errorbar(limits) Ultimamente estou ficando muito satisfeito com esta library. Para maiores detalhes, http://had.co.nz/ggplot2/ Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

função diff com variavel data tempo
by Vinicius Brito Rocha 08 Ago '11

08 Ago '11

Estou trabalhando em um código onde basicamente preciso calcular as diferenças de ordem 1 para uma variavel com o formato "2011-08-08 12:30:00" "2011-08-09 13:30:00" . . . "2011-08-15 09:30:00" . . . O problema surgiu quando quando variava o parâmetor lag da função e armazenar o primeiro valor das diferenças. guarda=0 for ( i in 1:10){ guarda[i]=diff(vetor_data_tempo,lag=i)[1] } pois dependendo da distância/tempo entra linhas a notação da resposta muda de minuto, para hora e por aí vai.. perciso fixar esta notação, alguém saberia como fazer? Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

1 0

Leitura arquivo rtf
by Eder David Borges da Silva 06 Ago '11

06 Ago '11

Pessoal, tenho arquivos neste formato rtf (primo distante do .doc e CIA).... Arquivo em: http://www.leg.ufpr.br/~eder/ST-26-E.rtf gostaria de saber se alguém já conseguiu carregar um documento com estas características, pois aparentemente é uma tabela, porem parece que cada célula é uma caixa de texto, muito estranho. Quanto utilizo o readLines vem uma serie de comandos de formatação: 1] "{\\rtf1\\ansi{\\fonttbl{\\f0\\fcharset0\\fprq0 Arial;}}{\\colortbl;\\red0\\green0\\blue0;\\red255\\green255\\blue255}" Att Éder

3 5

Re: [R-br] Teste de Normalidade Multivariado! Erro!
by Thiago De paula protásio 06 Ago '11

06 Ago '11

Ok Ivan... Obrigado pela dica....Eu encontrei o erro aqui......Estava usando a função matrix() e na verdade no meu caso a função é as.matrix(). C=as.matrix(carvão[,(9:14)]) mvShapiro.Test(C) #Teste de normalidade multivariado Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Re: [R-br] Teste de Normalidade Multivariado! Erro!
by Ivan Bezerra Allaman 06 Ago '11

06 Ago '11

A mensagem de erro já diz tudo!! X <- matrix(rnorm(22),ncol=2) mvShapiro.Test(X) X1<-matrix(rnorm(24),ncol=2) mvShapiro.Test(X1) Sem um exemplo reproduzível fica difícil verificar qual o real problema dos seus dados. Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

Teste de Normalidade Multivariado! Erro!
by Thiago De paula protásio 06 Ago '11

06 Ago '11

Prezados membros da R-Br; Estou tendo um problema com o teste de normalidade multivariado da função mvShapiro.Test() do pacote mvShapiroTest Estou trabalhando com 6 variáveis e 19 observações por variável....Talvez, vcs possam me ajudar a resolver esse erro. C=matrix(carvão[,(2:7)])> is.matrix(D)[1] TRUE > C=matrix(carvão[,(2:7)])> mvShapiro.Test(C)Erro em mvShapiro.Test(C) : Sample size must be between 12 and 5000. Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Join my network on LinkedIn
by Luis Iván Ortiz Valencia via LinkedIn 06 Ago '11

06 Ago '11

LinkedIn ------------ Luis Iván Ortiz Valencia requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Luis Iván Ortiz Valencia http://www.linkedin.com/e/-33g73l-gr0nt2c7-1e/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View invitation from Luis Iván Ortiz Valencia http://www.linkedin.com/e/-33g73l-gr0nt2c7-1e/vhj1xG42iBdUFK5d_HMEGL4iFMC51… ------------------------------------------ DID YOU KNOW LinkedIn can help you find the right service providers using recommendations from your trusted network? Using LinkedIn Services, you can take the risky guesswork out of selecting service providers by reading the recommendations of credible, trustworthy members of your network. http://www.linkedin.com/e/-33g73l-gr0nt2c7-1e/svp/inv-25/ -- (c) 2011, LinkedIn Corporation

1 0

Como armazenar o valor-p de um teste de normalidade multivariada?
by Fernando Barbosa 06 Ago '11

06 Ago '11

Gostaria de comparar as taxas de erro tipo I e o desempenho do poder de dois testes de normalidade multivariada do R : um é a função mvShapiroTest do pacote mvnormtest e o outro teste é a função mvnorm.etest do pacote energy. Essa comparação será feita por meio de simulação Monte Carlo. Gostaria de saber se existe como armazenar o valor-p sem ter que implementar essas funções uma vez que elas já estão disposta no R. Como fazer isso?Tem como ver como esta função foi implementada no R?Grato.Fernando

1 0

Ajuda com a função sink()
by Maria Papa 04 Ago '11

04 Ago '11

Olá pessoal, boa tarde, novamente Eu estou com um problema para usar a função sink(). Há anos eu uso da mesma forma e nunca tive problema. Eu estou usando esta função da seguinte forma: # Resultados sink("Resultado_linha1.txt") cat("\nTabela: pontos medidos(X, Y e Z), Estimativas boot para as coordenadas X, Y e Z e seus respectivos intervalos de 95% de confiaça\n") show(dados_final) cat("\nDesvios observados entre as coordenadas medidas de X e as suas estimativas boot\n") show(desvio_peca_X) cat("\nDesvios observados entre as coordenadas medidas de Y e as suas estimativas boot\n") show(desvio_peca_Y) cat("\nDesvios observados entre as coordenadas medidas de Z e as suas estimativas boot\n") show(desvio_peca_Z) sink() e tenho a seguinte mensagem de erro: > sink("Resultado_linha1.txt") Erro em file(file, ifelse(append, "a", "w")) : não é possível abrir a conexão Além disso: Mensagens de aviso perdidas: In file(file, ifelse(append, "a", "w")) : não foi possível abrir o arquivo 'Resultado_linha1.txt': Permission denied Não estou conseguindo entender o motivo da permissão negada. Obrigada Célia

2 1

Re: [R-br] Ajuda com a função sink()
by Ivan Bezerra Allaman 04 Ago '11

04 Ago '11

Talvez ajude, http://r.789695.n4.nabble.com/file-reading-problem-unique-to-windows-Err-in… Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

Re: [R-br] Ajuda para gráfico 3D de uma fç não paramétrica
by Ivan Bezerra Allaman 04 Ago '11

04 Ago '11

Não paramétrica ou não-linear? Nos forneça um CMR. Caso não saiba, siga o guia de postagem em: http://www.leg.ufpr.br/doku.php/software:rbr Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

Ajuda para gráfico 3D de uma fç não paramétrica
by Maria Papa 04 Ago '11

04 Ago '11

Olá Pessoal, boa tarde, Gostaria da ajuda de vcs plotar um gráfico em 3D de uma curva não paramétrica. Obrigada Célia

2 1

operando em arquivos na lista de objetos
by Luis Iván Ortiz Valencia 04 Ago '11

04 Ago '11

Prezados tenho na lista de objetos arquivos do tipo [4] "cuiaba2001" "cuiaba2002" "cuiaba2003" [7] "cuiaba2004" sabendo que todos têm a mesma estrutura de colunas. preciso juntar estes arquivos num so arquivo denominado cuiaba (juntando as linhas) a partir de uma função que só use a palavra do UF neste caso cuiaba. o rbind resolveria chamando os 4 arquivos, mas não me interessa essa solução. atte Ivan ...............................................

2 2

Re: [R-br] Converter informações de linhas duplicadas em colunas
by Roney Fraga Souza 04 Ago '11

04 Ago '11

Benilton, Muito obrigado! Levaria vários dias para solucionar esse problema. Os comandos que vc sugeriu, reshape e grep, abrem novos horizontes para tratamento de dados. Essas dias também são úteis para todos que trabalham com a POF, Pesquisa de Orçamentos Familiares do IBGE. Att Roney

1 0

operando em arquivos na lista de objetos
by Roney Fraga Souza 04 Ago '11

04 Ago '11

Ivan, No Rcmdr tem a função mergeRows, acredito que é o que vc deseja. Caso o código abaixo não seja útil é interessante olhar no menu do Rcmdr. library(Rcmdr) cba1e2 <- mergeRows(cuiaba2001, cuiaba2002, common.only=FALSE) cba1e2e3 <- mergeRows(cba1e2, cuiaba2003, common.only=FALSE) cba1e2e3e4 <- mergeRows(cba1e2e3, cuiaba2004, common.only=FALSE) Att Roney

1 0

Join my network on LinkedIn
by Lincoln Neves via LinkedIn 04 Ago '11

04 Ago '11

LinkedIn ------------ Lincoln Neves requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Lincoln Neves http://www.linkedin.com/e/-33g73l-gqxqlhyc-31/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View invitation from Lincoln Neves http://www.linkedin.com/e/-33g73l-gqxqlhyc-31/vhj1xG42iBdUFK5d_HMEGL4iFMC51… ------------------------------------------ DID YOU KNOW your LinkedIn profile helps you control your public image when people search for you? Setting your profile as public means your LinkedIn profile will come up when people enter your name in leading search engines. Take control of your image! http://www.linkedin.com/e/-33g73l-gqxqlhyc-31/ewp/inv-22/ -- (c) 2011, LinkedIn Corporation

1 0

Plot com dois eixos Y
by Marcelo Claro de Souza 04 Ago '11

04 Ago '11

Gostaria de plotar uma figura onde os meses ficassem no eixo X, a precipitação ficasse na forma de barras com a escala Y do lado esquerdo e a temperatura ficasse na forma de linha com a escala do lado esquerdo. Os dados seguem abaixo: mês prec tmed Jan 161.2 22.8 Feb 158.3 23.4 Mar 155.6 23.5 Apr 142.4 21.7 May 39.0 18.4 Jun 22.4 18.4 Jul 7.0 18.4 Aug 34.4 20.5 Sep 50.5 20.6 Oct 62.8 23.4 Nov 113.3 23.0 Dec 185.8 23.2Como devo proceder após usar barplot(mês~prec)? Muito obrigado.

3 2

Converter informações de linhas duplicadas em colunas
by Roney Fraga Souza 03 Ago '11

03 Ago '11

Caros, Estou com um problema que não consegui resolver, caso alguém saiba um caminho para solucioná-lo será uma grande ajuda. Vamos lá, no data.frame (inv): inv <- data.frame( seq=c(1,1,1,1,1,1,1,2,2,2,2,2,2), dv=c(9,9,9,9,9,9,9,7,7,7,7,7,7), dom=c(1,1,2,2,2,4,4,2,2,3,3,4,4), veiculo=c(301,101,201,101,301,401,101,101,301,301,101,401,201), qtde=c(1,2,2,2,1,1,4,1,1,3,1,2,3), ano=c(2002,2008,1993,2007,1993,2003,2000,2008,2007,2000,2003,2005,1994), moto_bic=rep(0,13), carro=rep(0,13), tt_vei=rep(0,13), ano_m=rep(0,13)) para facilitar a visualização, seq dv dom veiculo qtde ano moto_bic carro tt_vei ano_m 1 9 1 301 1 2002 0 0 0 0 1 9 1 101 2 2008 0 0 0 0 1 9 2 201 2 1993 0 0 0 0 1 9 2 101 2 2007 0 0 0 0 1 9 2 301 1 1993 0 0 0 0 1 9 4 401 1 2003 0 0 0 0 1 9 4 101 4 2000 0 0 0 0 2 7 2 101 1 2008 0 0 0 0 2 7 2 301 1 2007 0 0 0 0 2 7 3 301 3 2000 0 0 0 0 2 7 3 101 1 2003 0 0 0 0 2 7 4 401 2 2005 0 0 0 0 2 7 4 201 3 1994 0 0 0 0 temos informações de residências e veículos, as residencias são determinadas por combinações das variáveis (seq dv e dom), assim, a sequencia 1 9 1 indica uma residencia, e a combinação 1 9 2 outra residencia, por fim, temos a residencia de número 2 7 4. A variável veiculo representa os veículos de cada residência, sendo: 101 bicicleta 201 moto 301 carro 401 outros veículos Como em uma casa podem existir vários veículos, as informações que definem uma residencia (seq dv e dom) foram repetidas ao longa das linhas. Exemplo, na residência 2 7 2 tem uma bicicleta e um carro, cada um veículo em uma linha. A variável qtde indica quantos veículos de determinado tipo existe em uma residencia e a variável ano indica o ano de fabricação do veículo. Meu objetivo principal e dispor as informações dos veículos em novas colunas, de modo que elimine a duplicidade das residencias e ao mesmo tempo não perca informações sobre os veículos. O resultado final que busco seria: seq dv dom moto_bic carro tt_vei ano_m 1 9 1 2 1 3 2006 1 9 2 4 1 5 1998 1 9 4 4 0 5 2000 2 7 2 1 1 2 2008 2 7 3 1 3 4 2002 2 7 4 3 0 5 1998 Para chegar no resultado acima é necessário, na variável moto_bic: i) obter a quantidade de motos e de bicicletas da respectiva residencia na variável carro: ii) informar a quantidade de carros da respectiva residencia na variável tt_vei: iii) informar o total de veículos da residencia, valendo notar que como existe outros veículos, apenas a soma das variáveis moto_bic com carro não é suficiente na variável ano_m: iv) informar o ano médio dos veículos de cada residencia v) por fim, é necessário eliminar as residencias duplicadas. Obs.: o data.frame (inv) criado para exemplificar o problema é pequeno, mas chego a utilizar data.frames com mais de 1 milhão de linhas. Desde já grato Roney

2 1

Intervalos de Confiança
by Fernanda Balbino 03 Ago '11

03 Ago '11

Olá, boa tarde a todos, Depois de calcular o intervalo de confiança para uma amostra existe uma forma de verificar no R quantos pontos da amostra estão contidos no intervalo? T1 818,17542 956,86414 987,78540 1026,28223 1045,81201 1049,56323 1049,88818 1063,65503 1086,95178 1088,12463 1090,06958 1150,73083 1155,69067 1166,01892 1184,06067 1201,17676 1210,82520 1212,08447 1219,90637 1239,87903 1252,41418 1267,99792 1273,25818 1363,03088 1504,01331 IC (1098,511 ; 1194,63) Obrigada. Fernanda.

4 3

Re: [R-br] IDE R - Eclipse + StatET
by Roney Fraga Souza 03 Ago '11

03 Ago '11

Leandro, Acredito que esse post será útil para esclarecer suas dúvidas. http://r-br.2285057.n4.nabble.com/R-br-Eclipse-R-interface-grafica-para-R-t… Att Roney

2 1

IDE R - Eclipse + StatET
by Leandro Marino 03 Ago '11

03 Ago '11

Caros, Estou tentado a testar novas IDEs para o R. Já estou fazendo alguns testes com o RStudio no entanto, gostaria de testar junto ao Eclipse. Bom primeiramente qual a versão do Eclipse devo instalar.. São muitas opções e fiquei meio confuso! · Eclipse IDE for Java EE Developers, 212 MB · Eclipse Classic 3.7, 174 MB · Eclipse IDE for Java Developers, 122 MB · Actuate BIRT Designer Pro Promoted Download · Eclipse IDE for JavaScript Web Developers, 109 MB · Eclipse IDE for Java and Report Developers, 252 MB · Eclipse Modeling Tools, 271 MB · Eclipse for RCP and RAP Developers, 181 MB · Eclipse for Testers, 86 MB · Eclipse IDE for C/C++ Developers (includes Incubating components), 107 MB · Eclipse for Scout Developers, 174 MB · Eclipse IDE for Parallel Application Developers (includes Incubating components), 178 MB Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

1 0

Re: [R-br] Testes de normalidade multivariada do R
by Ivan Bezerra Allaman 03 Ago '11

03 Ago '11

Ou, library(sos) ???''multivariate normality test'' Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

inverter um data frame
by Mauricio Cardeal 03 Ago '11

03 Ago '11

Olá RMasters, Tenho um data frame que está assim: v1 <- c("id","A","B","C") v2 <- c(1,1,2,1) v3 <- c(0,0,1,1) dat <- as.data.frame(cbind(v1,v2,v3)) dat v1 v2 v3 1 id 1 2 2 A 1 0 3 B 0 1 4 C 1 1 Gostaria de deixá-lo assim: id A B C 1 1 0 1 2 0 1 1 Como faço para inverter a partir de dat ? Já tentei o reshape, mas não tive êxito. Vi uma resposta anterior na lista, mas também não consegui reproduzir. Obrigado, Mauricio UFBA

4 7

Testes de normalidade multivariada do R
by Fernando Barbosa 02 Ago '11

02 Ago '11

Olá a todos, recentemente questionaram meu trabalho em um evento, perguntando: Quais são os testes de normalidade multivariada que já possuem suas funções já implementadas no R? Sem sucesso, até então não consegui resposta para essa questão. Portanto, alguém saberiam me responder essa questão ou onde posso obter essa resposta. Grato pela atenção. Att. Fernando Costa

2 1

Re: [R-br] Sites para colocar arquivos para ser lido no R
by Ivan Bezerra Allaman 02 Ago '11

02 Ago '11

Opa!! agora sim! Valeuuuuuuuuuuuu Gustavo. Agora eu e os demais usuários temos duas opções para compartilhar arquivos de forma que estes possa ser lido pelo R. Obrigado a todos, estas dicas serão muito úteis a todos nós. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

3 2

Passar valores ao programa
by Carlos Mendonça 02 Ago '11

02 Ago '11

Caros, eu queria saber se existe no R alguma "opção" para se passar valores ao programa. Explicando melhor, os meus programas lêem e gravam vários arquivos. Os nomes desses arquivos contém as informações de ano, mes, etc. como no exemplo a seguir: tm0411 = read.csv2("C:\\2009\\Dados Abril 2011.csv") Eu precisava de alguma coisa do tipo informar no início do programa esses valores de ano e de mes e que substituisse no nome do arquivo automaticamente, mais ou menos assim: ano = 2011 mes = "Abril" ... ... ... ... tm0411 = read.csv2("C:\\*&ano*\\Dados *&mes &ano*.csv") É possível? Obrigado, Carlos Mendonça

4 9

extrair colunas como linhas
by Ricardo Pereira 02 Ago '11

02 Ago '11

Boa noite, Agora tenho o seguinte problema. Tenho um arquivo de dados como abaixo Z <- read.table("/test.txt", sep="\t", stringsAsFactors = FALSE) Z: QW B1 4 2 QW B1 4 2 QW B1 4 2 QW D2 2 1 QW D2 2 1 QW D2 2 1 QW D3 1 3 QW D3 1 3 QW D3 1 3 preciso escrever esse arquivo "test.txt" em um outro arquivo, ou seja, SALVAR como "final.txt". 1) Selecionar o elemento da 2a coluna e escrever na 1a linha; 2) Escrever todos os elementos da 3a que correspondem ao intem 1, em linha; 3) Escrever todos os elementos da 4a que correspondem ao intem 1, em linha; 4) Retorna ao intem 1 (ciclo) Obs. No arquivo test.txt tenho muitas linhas + 500. A 1a é descartada... B1 4 4 4 2 2 2 D2 2 2 2 1 1 1 D3 1 1 1 3 3 3 Muito obrigado caros amigos da lista. Ricardo Pereira

2 3

[Dúvida] Como pegar uma substring do fim para o inicio
by Pedro Rafael 02 Ago '11

02 Ago '11

Pessoal como pegar extrair uma substring que comece a contar do final da string? Por exemplo pegar os 4 ultimos caracteres da string. Pedro

4 6

Um convite especial
by Mauro Sznelwar 02 Ago '11

02 Ago '11

Olá R-br, Criei um perfil no Facebook com minhas fotos, vídeos e eventos e quero adicionar-lhe aos amigos para que você possa ver meu perfil. Primeiro, você precisa cadastrar-se no Facebook! Uma vez cadastrado, você também poderá criar seu próprio perfil. Obrigada, Mauro Para se cadastrar no Facebook, clique no link abaixo: http://www.facebook.com/p.php?i=1531843709&k=AQCzBSKCGLaLypSQBiUhcJ009eWjuf… Já possui uma conta? Adicionar este endereço de e-mail à sua conta: http://www.facebook.com/n/?merge_accounts.php&e=r-br%40listas.c3sl.ufpr.br&… ======================================= r-br(a)listas.c3sl.ufpr.br foi convidado a participar do Facebook por Mauro Sznelwar. Se você não deseja mais receber estes e-mails do Facebook, use o link abaixo para cancelar sua inscrição. http://www.facebook.com/o.php?k=5e824d&u=100002736790529&mid=4a1f49fG5af3b3… Saiba mais sobre esse e-mail: http://www.facebook.com/help/?faq=17151\nFacebook, Inc. P.O. Box 10005, Palo Alto, CA 94303

1 0

Re: [R-br] eliminar linha completa de amostra rejeitada
by Ricardo Pereira 01 Ago '11

01 Ago '11

Olá Benilton! Estou usando o RStudio. e retorna "i" como NA_integer_, tb cria C com apenas uma linha com todos os valores NA. Como pode ver sou novato no R... Obrigado pelo help. Abraço Ricardo

2 4

RStudio
by Gilson Sanchez 01 Ago '11

01 Ago '11

Pessoal estou iniciando a usar o RStudio, como eu faço para exportar os script do R para o RStudio. obrigado -- Eng. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

5 14

Trabalhando com shapefile dentro do spatstat
by Alexandre Santos 01 Ago '11

01 Ago '11

Boa noite pessoal, Estava tentanto usar um objeto SpatialPolygonsDataFrame contido em um arquivo shapefile dentro do pacote spatstat, bom fiz segundo o exemplo do material Handling shapefiles in the spatstat package do Baddeley, mas não consegui plotar o objeto owin criado, fiz: require(maptools) require(spatstat) x <- readShapeSpatial("mtaremanne92.shp") class(x) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp" cregions <- slot(x, "polygons") cregions <- lapply(cregions, function(x) {SpatialPolygons(list(x))}) cwindows <- lapply(cregions, as.owin) Checking 1 polygons...1. done. Checking 3 polygons...1, 2, 3. done. Checking for cross-intersection between 3 polygons...1, 2. done. plot(cwindows) Erro em xy.coords(x, y, xlabel, ylabel, log) : 'x' is a list, but does not have components 'x' and 'y' Alguem que já tenha utilizado shapfiles dentro do spatstat saberia me dar um help, Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

1 0

Dúvida - Operadores %any%, :: e :::
by Pedro Rafael 01 Ago '11

01 Ago '11

Pessoal qual é a finalidade dos operadores %any%, :: e ::: ? Procurei no Help mas só vi a informação que segue no link ( http://stat.ethz.ch/R-manual/R-patched/library/base/html/Syntax.html) Praticamente poderiam me dar um exemplo? Pedro Rafael

3 3

[Dúvida] Período das versões do R
by Pedro Rafael 01 Ago '11

01 Ago '11

Pessoal alguém sabe me informar o período de atualização das versões do R, ou seja, se existe um espaço temporal definido para cada versão do R? Pedro Rafael

4 3

Re: [R-br] Aos interessados em execucao paralela de codigos em R
by Junior Beleti 01 Ago '11

01 Ago '11

A execução "automática" com certeza é muito bem vinda. Já ouviu falar no rparallel? Pelo que li, ele seria uma solução para execução paralela "automática", apenas com a inclusão de algumas diretivas. Porém, não consegui nem instalá-lo. Abraço. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

Aos interessados em execucao paralela de codigos em R
by Benilton Carvalho 01 Ago '11

01 Ago '11

O codigo-fonte do R a ser lancado em Outubro sugere a inclusao de um pacote chamado 'parallel' e (pelo menos inicialmente) mantido pelo Brian Ripley. Ainda e' muito cedo para afirmar detalhes, mas a impressao que tive e' que funcionara' como o 'multicore' e um dos comandos providos sera' (finjam surpresa) 'parlapply'. Vantagens de um sobre o outro (ou sobre qualquer outro front-end de sua preferencia): impossivel dizer agora. Mas a relevancia da noticia e' que se o Ripley esta' usando seu tempo com o 'parallel', em breve podemos comecar a ver rotinas comuns fazendo uso automagicamente de recursos de computacao paralela. b

1 0

Re: [R-br] [OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 01 Ago '11

01 Ago '11

Walmes, Estou impressionado com a ferramenta para converter .tex para .odt. O negócio é fantástico mesmo. É claro que não sai como o arquivo pdf. No entanto, o arquivo .odt fica muito bom. As fórmulas matemáticas ficaram perfeitas. As saídas do R ficaram perfeitas. Só não descobri, porque os gráficos não saíram. Vou tentar descobri e posto de novo. Aconselho a todos fazer um teste para ver como fica. Valeu Walmes pela excelente dica de novo. Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

[OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 01 Ago '11

01 Ago '11

Bom dia senhores!! Alguém sabe se tem como alterar a extensão pdf dos gráficos gerados pelo sweave para outra extensão de interesse? Por exemplo, ao compilar um arquivo sweave, gostaria que os gráficos gerados fossem em *.tiff*. Quando há poucos gráficos no arquivo dá para fazer manualmente mesmo, no entanto, quando há muitos gráficos, mudar a extensão dos gráficos já se torna penoso. Abraços! Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

5 5

Re: [R-br] [OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 31 Jul '11

31 Jul '11

Eu acho engraçado como alguns anos de convívio é o suficiente para algumas pessoas já descobrirem exatamente onde queremos chegar. rsrsrsr. É isso aí mesmo Walmes. É exatamente um relatório de análises de dados. É exatamente para uma pessoa com usa Windows e precisam de figuras JPG e lalala. É claro que para nós que usamos linux, os gráficos gerados em pdf são espetaculares. Mais na nossa área, e acredito que mais especificamente na área (retrógrada zootécnica), as revistas exigem figuras em formatos, png, tiff, lalala. As suas sugestões foram valiosas, como sempre. Só não entendi ainda como funciona os comandos $convert foo.pdf foo.png Que 'foo' é esse? Utilizei estes comandos com o meu arquivo relatorio.pdf, e ele gerou cada página do relatório em png. Mas descobri o seguinte, por acaso, juntamente com o Fábio. Se eu colocar: options(device=png) A compilação do arquivo Rnw me gera gráficos em pdf e também em png. Mais uma vez obrigado. Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

Re: [R-br] [OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 31 Jul '11

31 Jul '11

Excelente Faria! Buscava exatamente isso! Agora posso compilar do vim mesmo. Abraço! Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

3 2

Re: [R-br] [OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 31 Jul '11

31 Jul '11

Agora sim!! É que é trivial para uns, para outros nem tanto. Obrigado de novo Benilton! Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

Re: [R-br] [OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 31 Jul '11

31 Jul '11

Aí você falou grego agora!!! Reformulando a pergunta: Estou usando o vim para compilar os arquivos. Mas vamos supor que eu esteja compilando do terminal, qual comando eu devo utilizar para compilar um arquivo sweave, cujo os gráficos gerados saiam com uma extensão PNG por exemplo? Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

Re: [R-br] Cluster + ANOVA
by Thiago De paula protásio 30 Jul '11

30 Jul '11

Pessoal; Muito legal o tópico!Vários membros do grupo levantaram a questão da validação do agrupamento.No caso de agrupamentos hierárquicos aglomerativos podemos utilizar a correlação cofenética...ou mesmo a MANOVA como foi mencionado. No caso da comparação múltipla multivariada.....quando temos vários tratamentos ou grupos formados...podemos utilizar a aproximação de Bonferroni ou T² de hotelling. Mas acho que no R não existam funções implementadas pra isso... Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

[OFF TOPIC] "TEOREMA DOURADO" - Jacob Bernoulli (1731)
by Marcelo Luiz de Laia 29 Jul '11

29 Jul '11

Pessoal, sou assiduo leitor do Nassif e me deparei com esse post dias atras. Fui pesquisar e vi que ha relacao. Mas, vou perguntar aqui, haja vista a grande quantidade de estatisticos na lista. Me desculpem o OFF TOPIC. http://www.advivo.com.br/node/553641 O post completo esta aqui: http://www.advivo.com.br/blog/luisnassif/vox-populi-haddad-estreia-com-3 -- Marcelo Brazil Linux user number 487797

1 0

Cluster + ANOVA
by Vinicius Brito Rocha 29 Jul '11

29 Jul '11

Srs. Estou lendo um livro de multivariada onde existe uma avaliação dos cluster produzidos pelo K-means através de ANOVA. Não entendi muito bem o que estava sendo feito. Num primeito momento me pareceu que para cada variável realizava-se uma ANOVA, e amarzenava-se as somas de quadrados individualmente. Alguém já viu isso? Ou melhor. Alguem já avaliou cluster via ANOVA? Att. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

5 10