- R-br - DINF UFPR

Limitação de memória
by Edson Lira 15 Mar '12

15 Mar '12

Bom dia a todos, estou executando uma rotina que lê informações de doadores de sangue. Uso uma base sql, a rotina está me dando a seguinte mensagem: Erro: não é possível alocar vetor de tamanho 62.5 Mb Meu sistema tem a seguinte configuração: Dual Core 2,7 Ghz Memoria: 2 giga S.O de 32 bits O que faço? Edson Lira Estatístico Manaus-Amazonas

2 1

Manipulação de dados (strings)
by Estevam Caixeta 15 Mar '12

15 Mar '12

Boa noite, possuo um banco de dados com 3 variáveis: 1- Nome da Pessoa 2- Data de nascimento 3- Nome completo da mãe Uma melhor visão do banco pode ser dada assim: NOM_PESSOA DTA_NASC_PESSOA NOM_COMPLETO_MAE_PESSOA 1 *A*NTONIO *D*O *N*ASCIMENTO VIEIRA 17081968 GENESIA *M*ARIA *D*O NASCIMENTO 2 ANTONIA* M*ARGARETE *R*OCHA *D*E ARAUJO 05091975 MARIA *D*A *A*SSUNCAO ROCHA 3 JOAO *P*EDRO *A*RAUJO VIEIRA 06031998 ANTONIA *M*ARGARETE *R*OCHA *D*E ARAUJO 4 ERICA *D*E *A*RAUJO *V*IEIRA 03112001 ANTONIA *M*ARGARETE *R*OCHA *D*E ARAUJO o que eu gostaria de fazer é criar uma quarta variável com a seguinte característica: PrimeiroNomeÚltimoNome*IniciaisNomeDoMeio* DataNascPrimeiroNomeMaeUltimoNomeMae*IniciaisNomeDoMeioMa*e para poder identificar cada indivíduo, já que no meu banco real eu possuo muitos homônimos. Para exemplificar, utilizando os nomes acima fornecidos, eu gostaria de criar algo do tipo: 1 ANTONIOVIEIRA*DN*17081968GENESIANASCIMENTO*MD* 2 ANTONIAARAUJO*MR*D05091975MARIAROCHA*DA* 3 JOAOVIEIRA*PA*06031998ANTONIAARAUJO*MRD* 4 ERICAVIEIRA*DA*03112001ANTONIAARAUJO*MRD* Desde já agradeço. Estevam C.M. Teixeira PS: Coloquei a parte relacionada aos nomes do meio para tentar facilitar a visualização e compreensão do que eu gostaria de fazer

4 3

Stepwise em glm
by Vinicius Brito Rocha 14 Mar '12

14 Mar '12

Senhores, me passou uma dúvida pela cabeça. Stepwise funciona em modelog glm? 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.

3 2

Erro ao plotar 2 climatogramas no mesmo frame
by Marcelo Claro de Souza 14 Mar '12

14 Mar '12

Tentei plotar 2 climatogramas em um mesmo frame, mas não sei o pq que está sendo plotado uma única figura entre as colunas 1 e 2 ao invés de plotar uma figura em cada coluna. Estou usando o pacote "climatol". Segue abaixo o script e o conjunto de dados do próprio pacote. layout(matrix(1:2,nrow=1,ncol=2)) par(mar = c(2.5,5,1,1)) data(datcli) diagwl(datcli,est="Example station",alt=100,per="1961-90",mlab="en") Muito obrigado pela ajuda. Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

1 0

duvida sobre presunção de homocedasticidade
by Fernando Antonio de souza 14 Mar '12

14 Mar '12

Caros amigos Estou trabalhando com dados de desenvolcimento de fetos ovinos. Estou modelando o seu crescimento ao longo da gestação. Para tal tenho utilizado o modelo de gompertz da forma log(Y)= a-b*exp(-c*x) sendo x a idade gestacional em dias e Y o peso do feto. Os meus dados atendem o presuposto de normalidade, porém não o faz para o pressuposto de igualdade das variâncias. Tenho lido na literatura a possibilidade de modelar a variância nestes casos. A variância residual observada graficamente se eleva e em seguida rebaixa conforme as médias do modelo se elevam. Na própria literatura comenta que este tipo de comportamento pode ser modelado utilizando o método VarPower() na função glns do pacote nlme. Eu assim o fiz entretanto o novo modelo não obteve melhor ajuste que o anterior quando comparados pela função Anova (p-value e AIC,BIC). Minha dúvida é a seguinte sem a igualdade das variâncias eu não posso validar o modelo. Existe alguma forma de corrigir este problema? abaixo disponho do CMR library(nlme) dput(gestsimp) gompertz1<-nls(log10(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=8.0,b=10.40,c=0.0023),data=gestsimp) gompertz2<-gnls(log10(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=8.0,b=10.40,c=0.0023),weights=varPower(),data=gestsimp) summary(gompertz1) plot(log10(PESO.FETO)~Gest,data=gestsimp,main="Peso do feto em função da idade de cordeiros") curve(5.697742-8.360311*exp(-0.003531*x),add=T,col=2) caso o dput()não funcione segue a saída dele dput(gestsimp)structure(list(Animal = c(3L, 51L, 83L, 89L, 91L, 92L, 93L, 48L, 52L, 55L, 69L, 76L, 78L, 11L, 25L, 41L, 42L, 46L, 50L, 53L, 57L, 10L, 19L, 20L, 30L, 38L, 45L), Manejo = c(1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), Gest = c(90L, 90L, 90L, 90L, 90L, 90L, 90L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L, 140L, 140L, 140L, 140L), Fetos = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L ), PV = c(51.1, 46.4, 53.5, 58.2, 30.7, 43.2, 41.5, 42, 43.4, 48.1, 42.2, 47.5, 32, 64.1, 49.9, 40.5, 42, 40.4, 50, 49, 52, 60.1, 55.9, 57, 51.5, 45, 42.8), PCVZ = c(36.93, 36.135, 46.125, 48.02, 23.095, 33.75, 32.575, 34.99, 34.465, 39.578, 30.84, 38.045, 27.455, 54.797, 43.595, 32.12, 33.841, 33.16, 42.06, 40.605, 42.325, 50.436, 46.515, 48.845, 43.075, 36.55, 35.435), PESO.FETO = c(0.215, 0.555, 0.49, 0.53, 0.225, 0.465, 0.61, 1.465, 0.585, 1.755, 0.62, 1.615, 1.045, 0.985, 2.79, 2.42, 2.286, 2.41, 2.3, 3.1, 3.105, 5.38, 4.065, 4.258, 3.878, 3.49, 3.032), FETOGKG = c(60.5065, 50.8621, 67.9386, 32.7586, 41.7391, 46.1815, 47.7506, 51.0421, 66.0087, 63.5711, 50.8569, 60.5784, 49.8817, 62.1197, 60.1357, 0, 0, 62.2886, 0, 55.4768, 53.5301, 68.0826, 58.413, 59.6224, 61.1382, 50.2531, 62.5518 )), .Names = c("Animal", "Manejo", "Gest", "Fetos", "PV", "PCVZ", "PESO.FETO", "FETOGKG"), row.names = c(11L, 13L, 15L, 19L, 20L, 21L, 22L, 30L, 31L, 32L, 35L, 36L, 38L, 42L, 45L, 46L, 47L, 48L, 49L, 51L, 55L, 63L, 66L, 67L, 71L, 76L, 80L), class = "data.frame", na.action = structure(22L, .Names = "62", class = "omit"))

1 0

Krigagem em 3D
by Tiago Vieira 14 Mar '12

14 Mar '12

Bom dia a todos! Vocês sabem me informar se existe algum pacote do R que faça krigagem em 3D? Também aceito sugestões de literaturas, pois estou com dificuldades para encontrar. Desde já, agradeço pela atenção de todos! Att, Tiago Vieira.

2 1

Estatisticas de uso do R frente a outros programas
by Marcelo Laia 13 Mar '12

13 Mar '12

Prezado, Tempos atras foi postado aqui na lista um site ou paper com as estatisticas de uso do R frente a outros programas. Se nao me engano foi o Benilton quem postou. Mas, não estou achando a postagem. Alguem teria ela ai pronta ou teria os dados atuais? -- Marcelo

5 4

Abrir conexões com os nós de um sock cluster
by Lucas Barbosa 13 Mar '12

13 Mar '12

Boa noite pessoal, Eu estava tentando executar em paralelo uma função -- distancia() -- que baixa certas informações de páginas (html) e retorna alguns valores de meu interesse. Eu tentei fazer usando o snow em uma só máquina (sock cluster), mas todas as vezes que eu tento acessar a internet a partir de um "snow worker" resulta nesse erro: *Erro em checkForRemoteErrors(val) : 4 nodes produced errors; first error: não é possível abrir a conexão *Já deu tudo certo usando o multicore em uma máquina com linux, mas eu pretendo usar esse tipo de função em um cluster com várias máquinas e isso o multicore não resolve. Há algum tipo de opção que eu devo usar na makeCluster() para poder abrir esse tipo de conexão a partir dos nós do cluster?

3 5

Re: [R-br] Teste de Kruskal Wallis!
by Thiago de Paula Protásio 13 Mar '12

13 Mar '12

Muito obrigado pela dica Emmanuel Arnhold . Já peguei o livro que vc me sugeriu e de fato tem um exemplo muito bom sobre o teste. Att., -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

Re: [R-br] Teste de Kruskal Wallis!
by Thiago de Paula Protásio 13 Mar '12

13 Mar '12

Olá Fernando! Abaixo seguem os resultados. A linha de comando é relativamente simples para as funções dos dois pacotes. > require(agricolae)> kruskal(O, Clone, alpha = 0.05, group=F, main = "Clone") Study: Clone Kruskal-Wallis test's Ties or no Ties Value: 11.6483 degrees of freedom: 6 Pvalue chisq : 0.07029077 Clone, means of the ranks O replication G008 22.5 5 G084 16.4 5 G122 6.8 5 U030 14.4 5 U034 20.3 5 U060 26.8 5 U083 18.8 5 Comparison between treatments mean of the ranks Difference pvalue sig LCL UCL G008 - G084 6.1 0.301092 -5.7600575 17.96006 G008 - G122 15.7 0.011312 * 3.8399425 27.56006 G008 - U030 8.1 0.172796 -3.7600575 19.96006 G008 - U034 2.2 0.706832 -9.6600575 14.06006 U060 - G008 4.3 0.463864 -7.5600575 16.16006 G008 - U083 3.7 0.527988 -8.1600575 15.56006 G084 - G122 9.6 0.108468 -2.2600575 21.46006 G084 - U030 2.0 0.732354 -9.8600575 13.86006 U034 - G084 3.9 0.506096 -7.9600575 15.76006 U060 - G084 10.4 0.083256 . -1.4600575 22.26006 U083 - G084 2.4 0.681656 -9.4600575 14.26006 U030 - G122 7.6 0.199964 -4.2600575 19.46006 U034 - G122 13.5 0.027140 * 1.6399425 25.36006 U060 - G122 20.0 0.001776 ** 8.1399425 31.86006 U083 - G122 12.0 0.047532 * 0.1399425 23.86006 U034 - U030 5.9 0.316920 -5.9600575 17.76006 U060 - U030 12.4 0.041064 * 0.5399425 24.26006 U083 - U030 4.4 0.453640 -7.4600575 16.26006 U060 - U034 6.5 0.271122 -5.3600575 18.36006 U034 - U083 1.5 0.797474 -10.3600575 13.36006 U060 - U083 8.0 0.177984 -3.8600575 19.86006 > require(pgirmess)> kruskalmc(O~Clone, probs = 0.05, cont=NULL) Multiple comparison test after Kruskal-Wallis p.value: 0.05 Comparisons obs.dif critical.dif difference G008-G084 6.1 19.68897 FALSE G008-G122 15.7 19.68897 FALSE G008-U030 8.1 19.68897 FALSE G008-U034 2.2 19.68897 FALSE G008-U060 4.3 19.68897 FALSE G008-U083 3.7 19.68897 FALSE G084-G122 9.6 19.68897 FALSE G084-U030 2.0 19.68897 FALSE G084-U034 3.9 19.68897 FALSE G084-U060 10.4 19.68897 FALSE G084-U083 2.4 19.68897 FALSE G122-U030 7.6 19.68897 FALSE G122-U034 13.5 19.68897 FALSE G122-U060 20.0 19.68897 TRUE G122-U083 12.0 19.68897 FALSE U030-U034 5.9 19.68897 FALSE U030-U060 12.4 19.68897 FALSE U030-U083 4.4 19.68897 FALSE U034-U060 6.5 19.68897 FALSE U034-U083 1.5 19.68897 FALSE U060-U083 8.0 19.68897 FALSE Obrigado! -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

3 3

distância entre municípios
by jairan dionizio 13 Mar '12

13 Mar '12

Bom dia a todos ! Estou querendo saber como faço para usar distancia entre municipíos na análise de agrupamento. Jairan Dionizio Graduando em Estatística Skype:jairan.dionizio Facebook: Jairan Dionizio Twitter: @jairandionizio Tel.:(84) 9423-5286 "A força não provém da capacidade física e sim de uma vontade indomável." (Mahatma Gandhi) Cuidados ao enviar e-mail Ao reenviar apague o meu e-mail e o meu nome, apague também os endereços dos amigos antes de reenviar qualquer e-mail, encaminhe como cópia oculta (Cco ou Bcc) aos SEUS destinatários.Agindo sempre assim dificultaremos a disseminação de vírus, spams, banners e outros.

2 2

Re: [R-br] Teste de Kruskal Wallis!
by Thiago de Paula Protásio 13 Mar '12

13 Mar '12

Prezados membros da R-Br; Como nas postagens anteriores estou tentando encontrar respostas sobre o teste t que a função kruskal() pacote agricolae executa após o teste não paramétrico em questão. Contudo, analisando as postagens antigas dos demais membros percebi que existe um outro pacote que faz o teste de Kruskal-Wallis: pgirmess função kruskalmc (). O problema é que essa função retornou um resultado diferente do obtido pela função kruskal() do agricolae. Alguém poderia me ajudar a esclarecer essa dúvida.? Muito obrigado!! Att., -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

2 1

Re: [R-br] Teste de Kruskal Wallis!
by Thiago de Paula Protásio 13 Mar '12

13 Mar '12

Obrigado Allaman pela dica! Olhei com atenção a documentação do pacote, mas não encontrei uma explicação sobre o teste t-Student que é aplicado após a função kruskal(). Não mencionam se é um teste t simples ou com algum tipo de proteção como Bonferroni. Muito obrigado! Att., -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

Re: [R-br] Ajuda com modelos mistos (medidas repetidas)
by Ivan Bezerra Allaman 13 Mar '12

13 Mar '12

continuando as perguntas, d) o intervalo de 48 h entre as medidas foi o suficiente para torná-las independentes? (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}

2 1

Re: [R-br] Teste de Kruskal Wallis!
by Ivan Bezerra Allaman 13 Mar '12

13 Mar '12

Coisas interessantes em http://tarwi.lamolina.edu.pe/~fmendiburu/ talvez encontre resposta para sua pergunta! (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}

1 0

Ajuda com modelos mistos (medidas repetidas)
by Olympio Neto 13 Mar '12

13 Mar '12

Pessoal, gostaria de uma ajuda de vocês. Estou com um problema de estudo fatorial com medidas repetidas (assim considero) com as seguintes características: 1- Cada corpo de prova foi submetido a 1 tratamento com 3 níveis distintos, sendo que cada grupo continha 18 unidades amostrais (fator 1), 54 corpos de prova ao total; 2- Dentro de cada grupo definido pelo nível do tratamento no fator 1, metade das unidades amostrais do grupo (9) foram submetidas ao tratamento 2 com 2 níveis (fator 2), de modo que o i-ésimo indivíduo está submetido ao j-ésimo nível do fator 1 e k-ésimo nível do fator 2, ou seja, são independentes com relação ao fator 2. 3-A variável resposta foi medida em dois tempos distintos, tempo 0 e tempo 48h, por isso considerei medidas repetidas. Para cada unidade amostral, portanto, eu considerei um modelo assim: fixando-se a unidade amostra i,que foi somente submetida a 1 tratamento do fator 1 e 1 tratamento do fator 2. Logo, j e k também estão fixos e o que varia dentro de cada indivíduo é o tempo (l) em que a medida foi obtida. O modelo abaixo Y_l=b0+b1Trat1_l+b2*Trat2_l+b3*Tempo_l+todas as interações+a_i+e_il, l=1,2, a_i ~N(0, var_1) (efeito aleatório) e e_il~Normal(0, var_2) Desta maneira, as pseudorrepresentações são os tempos, com tamanho 2 para cada unidade e isso vem causando sérios problemas, porque acredito que o tratamento 2 deveria estar aninhado no tratamento 1, mas quando tento algo, os graus de liberdade são, obviamente insuficientes e um erro ocorre. Estimado o modelo acima, com a função lme, (algo do tipo lme(y~Trat1*Trat2*Tempo, random=1|Id)), os resíduos não são normais. Existe alguma maneira de utilizar a transformação de box-cox para tentar normalizar os resíduos. Eu fiz um aproximação usando o modelo linear comum com estes dados e a função boxcox do pacote MASS, e depois usei a transformação indicada no modelo linear neste misto, variando os valores de lambda em torno do valor estimado e vendo qual se aproximava mais da normalidade. Ficou muito ruim, ainda. Então, lanço as minhas perguntas: a) Com a experiência de vocês me digam, por favor; estou usando um modelo adequado para que quero encontrar? b) Existe algum pacote no R que calcule a transformação de box-cox em objetos lme? c) Será que o fator 2, dada natureza do experimento, deve ser aninhada ao fator 1 como efeito fixo? Algo como Y~Trat1*Trat1/Trat2*Tempo? Sei que estas perguntas fogem ao escopo deste fórum, mas logo imaginei que vocês pudessem me ajudar. Desde já muito obrigado. Olympio

1 0

Teste de Kruskal Wallis!
by Thiago de Paula Protásio 12 Mar '12

12 Mar '12

Prezados membros da R-Br; Novamente recorro ao conhecimento de vcs. Estou tendo algumas dúvidas sobre a função kruskal() do pacote agricolae. Essa função executa o teste de kruskal e posteriormente um teste t. Recorri ao help() para entender esse teste t, mas não obtive sucesso. Gostaria de saber se alguém poderia confirmar ou não se esse teste t tem algum tipo de proteção como Bonferroni ou se é apenas um teste t simples (nesse caso não seria o mais apropriado eu acho). Alguém poderia indicar uma função que executa o teste de kruskal e posteriormente o teste t com a proteção de Bonferroni? Grato pela valiosa ajuda. -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

problemas com função nls
by Fernando Antonio de souza 11 Mar '12

11 Mar '12

caros amigos estou tentando rodar um modelo não linear de Gormpertz da forma Log(Y) = a-b*exp(-c*X) utiliando a função nls. A função é como segue: gompertz<-nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=3.4,b=22.99,c=0.018),data=gestsimp) dput(gestsimp) Acontece que toda vez que rodo, recebo o aviso de erro:" Erro em nls(log(PESO.FETO) ~ a - b * exp(-c * Gest), start = list(a = 3.4, : número de iterações excedeu o máximo de 50 Ja tentei elevar o máximo de interação utilizando o comando controle mas nunca converge. Alguem pode me esclarecer em que estou errando. os valore iniciais foram obtido sutilizando a função manipulate. um grande abraço

2 1

erro "initRefFields"
by D'Apolito 10 Mar '12

10 Mar '12

Olá, Alguém sabe como solucionar o erro: "Error in initFields(scales = scales) : could not find function "initRefFields" " ...sem ter que reinstalar o R 2.14.1? obrigado Carlos

2 7

Re: [R-br] problemas com função nls
by Ivan Bezerra Allaman 10 Mar '12

10 Mar '12

Provavelmente seus chutes estão errados! Acesse o site http://ridiculas.wordpress.com/ que provavelmente irás resolver o seu problema. Do contrário, nos forneça um CMR. (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 8

fatorial 5x3 + testemunha
by Marcelo Claro de Souza 10 Mar '12

10 Mar '12

Preciso rodar uma análise fatorial 5x3 + testemunha e estou usando i seguinte script segue em anexo o conjunto de dados para tornar o script reproduzível mamona=read.csv2(file.choose(),h=T) mamona.f1=TukeysHSD(aov(IVG~palha:test,data=mamona)) mamona.f2=TukeysHSD(aov(IVG~posicao:test,data=mamona)) mamona.f1.f2=TukeysHSD(aov(IVG~palha*posicao:test,data=mamona)) Vcs poderiam me informar se o erro está no script ou se meu conjunto de dados está organizado errado? Muito onrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Fwd: [R-pkgs] ggmap : ggplot2 and RgoogleMaps
by Benilton Carvalho 09 Mar '12

09 Mar '12

---------- Forwarded message ---------- From: David Kahle <david.kahle(a)gmail.com> Date: 5 March 2012 17:20 Subject: [R] [R-pkgs] ggmap : ggplot2 and RgoogleMaps To: r-packages(a)r-project.org, ggplot2 <ggplot2(a)googlegroups.com> Dear useRs - ggmap v1.2 is now available on CRAN. ggmap is a relatively new package which combines the power of ggplot2 with the spatial contextual information in Google Maps or OpenStreetMaps (via RgoogleMaps) to provide an easy, consistent and modular framework for spatial graphics. For examples on how to get started with ggmap, see ?ggmapplot for several examples of its use. Happy mapping! Cheers david kahle. ------------------------------------------------------------- ------------------------------------------------------------- David J. Kahle, Ph.D. Assistant Professor of Statistical Science Department of Statistical Science Baylor University One Bear Place #97140 Waco, Texas 76798 Office : (254) 710-6102 http://sites.google.com/site/davidkahle/

1 0

arvore de decisão
by Vinicius Brito Rocha 09 Mar '12

09 Mar '12

Senhores, estou com um problema que me parece de uma árvore de decisão, mas não sei se é popssível resolve-lo da forma que gostaria no R. É o seguinte tenho observações diárias de uma árvore, algo mais ou menos assim(extamente assim em uma planilha excel). Data Nível 1 Nível 2 Nível 3 Nível 4 Nível 5 Dia 1 - (Nível 1, Evento 2) - (Nível 2, Evento 1) - (Nível 3, evento 6) Dia 2 - (Nível 1, Evento 3) - (Nível 2, Evento 1) - (Nível 3, evento 1) Dia 3 - (Nível 1, Evento 3) - (Nível 2, Evento 1) - (Nível 3, evento 1) - (Nível 4, evento 1) . . . ou seja, Tenho 5 níveis (não preciso percorrer necessáriamente os 5 níveis) onde dentro de cada nível tenho uma gama de possibilidades, e ao final uma decisão é tomada. Minha pergunta é se existe alguma forma de ajustar uma árvore de decisão, onde eu force cada varíavel a classificar na orde que eu desejo, por exemplo: A oder de particionamento deve sempre seguir a ordem natural dos níveis, não posso começar a particionar pelo nível 3, sempre nível 1,2,3,4,5. 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

Dúvida: R para Windows 7 - 64bits
by Núbia Batista 09 Mar '12

09 Mar '12

Boa noite... Qual a versão do R mais apropriada para o Windows 7 - 64Bits? Instalei a última disponível no site, mas vários erros ocorrem durante a utilização! Agradeço desde já, Núbia Batista

6 7

Baixar figura
by Luís Gustavo 08 Mar '12

08 Mar '12

Pessoal, Estou com o seguinte problema, estou gerando um gráfico utilizando a API do Google Charts, mas agora quero salvar esta figura, como faço para baixá-la via R? http://chart.apis.google.com/chart?cht=bvs&chd=t:10,20,30,20,10&chs=170x100… Desde já agradeço. -- Luís Gustavo Silva e Silva

3 7

sobre a utilização de processamento em vários núcleos
by Antonio Silva 07 Mar '12

07 Mar '12

Olá, Rodo o R 64-bit no Ubuntu 10.04 em uma máquina com processador Core i7 e 8 Gb de RAM. Nesta semana começei a fazer uma análise de cluster, utilizando basicamente o pacote vegan, sobre uma matriz grande. O cálculo das ligações (hclust) levou três hora e meia para ser efetuado (!!!). Vi que apenas um núcleo do processador estava sendo utilizado. Começei então a procurar alternativas que permitissem o aceleramento dos cálculos e a utilização dos outros núcleos, que estão sub-utilizados. Li sobre o pacote multicore e outros, inclusive em mensagens nesta lista. Muitas coisas não ficaram claras para mim e, no final, não compreendi como proceder sua utilização no meu caso (ou mesmo se é possível). É possivel rodar uma análise de cluster fazendo utilização dos diversos núcleos? Agradeço qualquer orientação e/ou indicação de documentação. Abraços, Antônio -- Antônio Olinto Ávila da Silva Biólogo / Oceanógrafo Instituto de Pesca São Paulo, Brasil

4 7

Decision tree - espaço vertical dos nós
by Nathália Machado 06 Mar '12

06 Mar '12

Olá a todos! Estou trabalhando com a Árvore de Classificação (pacote *rpart*) e tenho uma dúvida sobre um dos argumentos da função "*plot.rpart*": Quando usamos "uniform=F", o espaço vertical dos nós será proporcional ao erro no ajuste, como explicado pelo manual do *rpart*. Porém, preciso construir uma escala que mostre isso para cada um dos nós. Então, como poderia ter acesso a esses valores de cada nó, para que possa construir a escala? Por exemplo: library(rpart) fit <- rpart(Price ~ Mileage + Type + Country, cu.summary) plot(fit, uniform=F) text(fit, use.n=TRUE) Cada nó tem um tamanho diferente de acordo com a explicação, certo? Com eu faço pra ver esse valores? Desde já agradeço! -- Nathália Machado PhD student - Ecology and Evolution Universidade Federal de Goiás Conservation Biogeography Lab Lab: wix.com/rdloyola/lab

1 0

[R-BR] Erro em rodar os variogramas empiricos com o pacote geoR
by Rodrigo Sant'Ana 06 Mar '12

06 Mar '12

Boa noite Listeiros da R-BR, estou com um problema no uso de algumas funções do pacote geoR, entretanto, acredito que seja pelo tamanho da base de dados. Infelizmente, não posso disponibilizar a base de dados para replicarmos aqui. Mas simulei algo parecido para que vocês tenham uma ideia. set.seed(125) x <- grf(10000, grid = "reg", xlim = c(-3.34, 8.03), ylim = c(-2.85, 6.32), cov.pars = c(122, 0.2), mean = 45) y <- variog(x, option = "cloud", estimator.type = "modulus", max.dist = 14.61) variog: computing omnidirectional variogram Erro: não é possível alocar vetor de tamanho 381.4 Mb Além disso: Mensagens de aviso perdidas: In if (class(geodata) %in% ls(pattern = glob2rx("Spatial*DataFrame"), : a condição tem comprimento > 1 e somente o primeiro elemento será usado Na base que estou utilizando para análise, não aparece essa mensagem de aviso. Só aparece o Erro de alocação de vetor de tamanho. Será que vocês podem me ajudar a resolver esse impasse? abraço a todos e muito obrigado! _______________ Rodrigo Sant'Ana Oceanógrafo

2 4

semivariograma cruzado
by DANIELA CUNHA DA SE 05 Mar '12

05 Mar '12

Olá gostaria de saber se alguém tem alguma rotina para fazer semivariograma cruzado? Att Daniela

2 1

Escala Gráfica para mapas
by Julianna Trindade 05 Mar '12

05 Mar '12

Pessoal, Gostaria de saber como plotar uma escala gráfica respeitando as dimensões corretas do meu shapefile. Ou seja, mostrar o tamanho da escala condizente com o tamanho em Km do meu mapa. Desde já obrigada, Aguardo resposta. -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

1 0

String com um espaço que não é espaço
by Lucas Barbosa 05 Mar '12

05 Mar '12

Bom dia pessoal, Estou escrevendo um código para calcular a distancia entre cidades brasileiras com o gmaps. Para construir os links eu preciso trocar os espaços por "+" nos nomes das cidades. Nos meus dados (uma lista de cidades) algumas strings contém um caracter que parece um espaço em branco, mas não se comporta como tal: > cidades[[47]][[3]] [1] "Águas Lindas de Goiás" > > cidade <- cidades[[47]][[3]] > > strsplit( cidade ,"[ \t\n\r\f\v]") [[1]] [1] "Águas Lindas" "de" "Goiás" > > gsub( "[ \t\n\r\f\v]" , "\\+", cidade) [1] "Águas Lindas+de+Goiás" > > loucura <- substr(cidade, 6, 6) > > gsub( loucura , "\\+", cidade) [1] "Águas+Lindas de Goiás" > > gsub( sprintf("[%s \t\n\r\f\v]", loucura) , "\\+", cidade) [1] "Águas+Lindas+de+Goiás" Eu consegui resolver o problema (trocar espaços por "+"), mas gostaria de escrever meu código sem extrair o caracter *"loucura"* com a *substr()*. Se não é *[ \t\n\r\f\v]*, o que é? Alguém tem uma dica? Abs., Lucas Cusinato

3 9

Apostila de R - Jackson Aquino.
by Pedro Rafael 04 Mar '12

04 Mar '12

Galera, segue o link de uma apostila em R em português explicando alguns tópicos importantes em R. Julianna o capítulo 11 é interessante que explica passo a passo como criar mapas estáticos em R. A apostila é bem simples pois é escrita para pessoas que não estão da área de computação/programação, contudo muito bem escrita e interessante. São poucas páginas e vale a pena uma leitura. LINK: http://www.lepem.ufc.br/jaa/RparaCS-260611.pdf -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

8 8

Dá pra Rodar o R nos Tablets ANDROID??
by Márcio Luís Moreira de Souza 04 Mar '12

04 Mar '12

Olá, amigos! O título da mensagem já diz tudo. Estou precisando de um aplicado ou algo do gênero que faça o R rodar em TABLETS ANDROID ou iPhone's mesmo... Alguém já ouviu algo por aí?? Tem alguma dica? Qualquer orientação é bem vinda... Abração!! -- Prof. Márcio Souza, M.Sc. (COPPE/UFRJ) *Pesquisador - UFJF e INFOTESTE BRASIL (JF)* *Professor - INSTITUTO VIANNA JÚNIOR (JF)* www.profmsouza.blogspot.com http://lattes.cnpq.br/4578008002785684

6 15

Analise Descritiva Datas periodicas
by Eder David Borges da Silva 04 Mar '12

04 Mar '12

Bom dia pessoal, Tenho um conjunto de dados, de fenologia, (fases da planta: florescimento, frutificação, etc), e periodicamente é dado um nota de 0 a 5 para percentual da fase, por exemplo entre 40 a 60 % de florescimento recebe nota 3, e assim por diante dependendo do percentual temos uma nota. Oque me interessa saber é ter primeiramente um descritivo, por exemplo de valor médio da fase, e IC para a fase, do tipo, com 95% de confiança o florescimento com nota 1 vai de 3 de janeiro a 10 de fevereiro. Ai surge minha duvidas como fazer estes cálculos com data no R, pois não é simplesmente pedir um summary por exemplo pois oque me interessa e o período. Um exemplo do conjunto de dados é simulados abaixo: ny <- 5 #numero de anos n <- 365*ny #dias inicio <- as.Date(c('01/01/2000'),'%d/%m/%Y') periodo <- seq(inicio,inicio+(n-1),by=1) y <- as.vector(sapply(1:ny,function(x)sort(ceiling(runif(n/ny,0,4))))) y1 <- as.vector(sapply(1:ny,function(x)sort(ceiling(runif(n/ny,0,4))))) dados <- data.frame(rep=rep(1:2,each=n),periodo,resp=c(y,y1)) head(dados) lattice::xyplot(resp~periodo,groups=rep,dados) Gostaria de um descritivo por nota, alguma sugestão a mais será muito bem vinda. Obrigado.

2 1

Problema ao instalar/carregar pacote mgarchBEKK
by Rodrigo Coster 04 Mar '12

04 Mar '12

Tentei carregar o pacote mgarchBEKK (baixei o .zip dele aqui: http://www.quantmod.com/download/mgarchBEKK/) mas aparece a seguinte mensagem de erro: > require(mgarchBEKK) Carregando pacotes exigidos: mgarchBEKK Failed with error: ‘package ‘mgarchBEKK’ does not have a NAMESPACE and should be re-installed’ Alguém sabe o que fazer (Obviamente já tentei reinstalar e o erro persistiu...)? Estou com o R 2.14.1, pelo que diz no site a ultima atualização foi pro R 2.12.0, mas acredito que isso não deveria interferir. Eu usei o mesmo .zip para instalar no pc da faculdade e la deu erro (não sei a versão do R de la). []'s

1 0

identificação de grupos em um cluster
by Antonio Silva 03 Mar '12

03 Mar '12

Olá Gostaria de saber se há como indicar no dendrograma, na altura do corte, o número referente ao grupo, como indicado em cutree. Observei que em alguns dendrogramas que os grupos selecionados de acordo com um dado o nível de corte não seguem necessariamente uma ordem ascendente (grupo 3 pode não ser o terceiro da esquerda para direita). Seria algo semelhante ao da figura https://www.crops.org/images/publications/cs/42/5/1584f4.jpeg Por exemplo exemplo: library(vegan) data(dune) dis <- vegdist(dune) cluc <- hclust(dis, "complete") plot(cluc) rect.hclust(cluc, 3) grp <- cutree(cluc, 3) grp Notem que o grupo com os objetos 17, 19, 11 e 18, posicionado a esquerda no dedrograma é o grupo de número 3, e o grupo com os objetos 1,2, ..., 7, representado no meio do dendrograma é o de número 1. Eu gostaria que, no nível de corte, aparecesse no dendrograma, ao lado das barras verticais de cada cluster, a indicação 3, 1, 2. Agradeço qualquer dica. Abraços Antônio -- Antônio Olinto Ávila da Silva Biólogo / Oceanógrafo Instituto de Pesca São Paulo, Brasil

3 4

Inserir valor de um objeto numérico dentro de uma expressão
by Daniel C Bezerra 03 Mar '12

03 Mar '12

Alguém saberia como inserir o valor de um objeto numérico dentro de uma expressão? Eu queria fazer um gráfico assim: hist(rnorm(n=100,mean=67,sd=37),main=expression("Media = 67"%+-%37)) mas se eu faço: media<-67 dp<-37 hist(rnorm(n=100,mean=67,sd=37),main=expression("Media ="media%+-%dp)) hist(rnorm(n=100,mean=67,sd=37),main=paste("Media =", media ,expression(""%+-%""),dp)) Não dá certo. Alguém saberia como resolver? Desde já mais uma vez agradeço, Daniel

3 4

RGL: uma dica de ouro para instalar dependências
by Jose Claudio Faria 03 Mar '12

03 Mar '12

Pessoal, O pacote RGL tem muitas dependências e ficar checando passo a passo as mensagens de erro do R para instalar todas as dependências é muito trabalhoso. Ao instalar o Debian (estava com saudades e decidi comparar com o LMDE no notebook) me deparei novamente com o problema e resolvi dar uma olhada na web. Achei uma dica de ouro: http://r.789695.n4.nabble.com/Pkg-rgl-installation-fails-because-of-x11-td8… Ben Bolker que compartilho: sudo apt-get build-dep r-cran-rgl ou sudo aptitude build-dep r-cran-rgl Facilita muito! Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 1

RES: Rosa dos Ventos
by Jodavid Ferreira 03 Mar '12

03 Mar '12

Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com

5 10

[OFF TOPIC] Julia
by Fernando Mayer 03 Mar '12

03 Mar '12

Achei interessante repassar para conhecimento. (Douglas Bates, na Rd) From: bates(a)stat.wisc.edu Douglas Bates To: r-devel(a)r-project.org R-devel Date: Thu, 1 Mar 2012 14:06:51 -0300 Subject: [Rd] Julia My purpose in mentioning the Julia language (julialang.org) here is not to start a flame war. I find it to be a very interesting development and others who read this list may want to read about it too. It is still very much early days for this language - about the same stage as R was in 1995 or 1996 when only a few people knew about it - but Julia holds much potential. There is a thread about "R and statistical programming" on groups.google.com/group/julia-dev. As always happens, there is a certain amount of grumbling of the "R IS SOOOO SLOOOOW" flavor but there is also some good discussion regarding features of R (well, S actually) that are central to the language. (Disclaimer: I am one of the participants discussing the importance of data frames and formulas in R.) If you want to know why Julia has attracted a lot of interest very recently (like in the last 10 days), as a language it uses multiple dispatch (like S4 methods) with methods being compiled on the fly using the LLVM (http://llvm.org) infrastructure. In some ways it achieves the Holy Grail of languages like R, Matlab, NumPy, ... in that it combines the speed of compiled languages with the flexibility of the high-level interpreted language. One of the developers, Jeff Bezanson, gave a seminar about the design of the language at Stanford yesterday, and the video is archived at http://www.stanford.edu/class/ee380/. You don't see John Chambers on camera but I am reasonably certain that a couple of the questions and comments came from him. --- Fernando Mayer Universidade Federal de Santa Catarina - UFSC Departamento de Ecologia e Zoologia - ECZ/CCB URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com

1 0

Formato SPSS
by Edson Lira 03 Mar '12

03 Mar '12

Boa noite pessoal, o Pacote gmodels tem a função CrossTable que mostra a saída com formatos SAS, SPSS. O detalhe é mostra somente para tabelas de dupla entrada. Tem algum pacote ou função que me mostro formato do SPSS para tabelas simples? Edson Lira Estatístico Manaus-Amazonas

2 1

Batch Mode Windows 7 64bits
by Leandro Marino 02 Mar '12

02 Mar '12

*Caros,* * * *estou tentando rodar a seguinte linha de comando no em um arquivo chamado Simulacao.bat* *"C:\Programas Leandro\R-2.14.0\bin\x64\" R BATCH --no-restore --no-save Simulacao.r log.txt* * * *Na pasta que o arquivo simulacao.bat se encontra tenho o Simulacao.r. Entretanto, o Rterm é aberto e recebo a tela, que está no link: https://lh3.googleusercontent.com/-NSLGLbsnlSQ/T1EOz0BHFyI/AAAAAAAABfE/FTFe… * * * *Alguém saberia me ajudar nisto?* * * *Aproveitando, quero depois passar alguns argumentos, como devo fazer?* * * *Procurei nos manuais e na lista mas não consigo entender meu erro.* * * *Att,* *Leandro Marino* * *

5 11

Batch mode Windows 7 64bits
by Leandro Marino 02 Mar '12

02 Mar '12

*Caros,* * * *estou tentando rodar a seguinte linha de comando no em um arquivo chamado Simulacao.bat* *"C:\Programas Leandro\R-2.14.0\bin\x64\" R BATCH --no-restore --no-save Simulacao.r log.txt* * * *Na pasta que o arquivo simulacao.bat se encontra tenho o Simulacao.r. Entretanto, o Rterm é aberto e recebo a tela, que está no link: https://lh3.googleusercontent.com/-NSLGLbsnlSQ/T1EOz0BHFyI/AAAAAAAABfE/FTFe… * * * *Alguém saberia me ajudar nisto?* * * *Aproveitando, quero depois passar alguns argumentos, como devo fazer?* * * *Procurei nos manuais e na lista mas não consigo entender meu erro.* * * *Att,* *Leandro Marino* * * * * Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

1 0

manipulação de colunas
by Rodrigo Plei 02 Mar '12

02 Mar '12

Prezados, Como faço para colocar em uma única colunas os dados que estão em diversas colunas com o mesmo número de linhas? Por exemplo, tenho um arquivo com uma tablela de dados assim col1 col 2 col 3 col 4 154 122 258 155 128 137 259 144 111 128 233 143 Precio que fique assim: coluna 154 128 111 122 137 128 258 259 233 155 144 143 Agradeço desde logo. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

3 4

[Dúvida] Atribuindo projeção no mapa.
by Pedro Rafael 02 Mar '12

02 Mar '12

plot(mapa, col = colors[findInterval(dados$ano_2008, brks,all.inside=TRUE)], axes = TRUE) Senhores, gostaria de inserir o sistema de referência espacial (4291) na hora de adicionar o vetor de shapfile ao R. PS: O .shp não possui sistema de referência atribuido. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 1

Duvida sobre como modificar dados em um data.frame
by Augusto Ribas 02 Mar '12

02 Mar '12

Ola pessoal. Estou com uma duvida simples aqui. Eu estava analisando dados aqui e vi um grafico estranho, dai vi que os dados estavam digitados errados. Mas o erro é simples. Existem 2 fatores, e os niveis de um fator estão trocados, mas somente dentro de um nivel do segundo fator, o resto esta tudo correto. Basicamente existe um vetor e eu tenho que trocar A por B e vive versa. Ai eu consegui arrumar imendando um monte de comando mas ficou uma coisa enorme e especifica para algo que parece tão simples. ai pensei se alguém tem uma solução mais economica e geral que a que eu fiz? ### #Exemplo de dados ### resposta<-c(c(rnorm(10,5),rnorm(10,8)),c(rnorm(10,8),rnorm(10,5))) fator1<-c(rep(c("A","B"),each=10),rep(c("A","B"),each=10)) fator2<-rep(c("1","2"),each=20) dados<-data.frame(resposta,fator1,fator2) ### #Grafico Errado ### library(lattice) bwplot(resposta~fator1|fator2,data=dados) ### #Minha Solução ### dados[which(dados$fator2%in%"2"),"fator1"]<-sapply(dados[which(dados$fator2%in%"2"),"fator1"],function(x){ if(x=="B"){ replace(x, x == "B", "A") } else { replace(x, x == "A", "B") } }) ### #Grafico Certo ### bwplot(resposta~fator1|fator2,data=dados) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

4 4

R + GPU (Era: Dá pra Rodar o R nos Tablets ANDROID??)
by Fernando Mayer 01 Mar '12

01 Mar '12

Jogando lenha na fogueira (: gputools promete otimizações incríveis usando a tecnologia CULA da nVidia http://brainarray.mbni.med.umich.edu/Brainarray/Rgpgpu/ Ainda não testei. Alguém? []s, --- Fernando Mayer Universidade Federal de Santa Catarina - UFSC Departamento de Ecologia e Zoologia - ECZ/CCB URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2012/3/1 Benilton Carvalho <beniltoncarvalho(a)gmail.com>: > programar para usar a GPU para fazer processamentos usualmente feitos em CPU. b > > 2012/3/1 Rodrigo Sant'Ana <rodrigo.gringo(a)gmail.com>: >> Concordo, em termos de custo beneficio final, ainda vale mais a pena >> manter um Notebook, ou até mesmo >> netbook, e concordo tbm, e foi exatamente isso que quis dizer no email >> anterior, que para "produção" não vale a pena ainda. >> Mas penso que como ferramenta secundária/auxiliar ao seu PC de >> produção possa trazer algum beneficio por agora, só isso. >> >> Mas você está certo. Já havia entendido o que você quis dizer, só >> estava ressaltando que esta >> é uma área de inovação para atuação do R e até mesmo como uma área de >> desenvolvimento em ascensão. >> >> Quando você diz em programar para GPU, você diz para optimizarmos >> funções para utilizar melhor as funções gráficas da placa, ou >> para utilizar o processador da placa de video como um processador fisico? >> >> >> >> >> _______________ >> Rodrigo Sant'Ana >> Oceanógrafo >> >> >> >> Em 1 de março de 2012 18:36, Benilton Carvalho >> <beniltoncarvalho(a)gmail.com> escreveu: >>> Rodrigo, >>> >>> o que eu quero dizer eh que, se vc for usar um dispositivo movel com >>> teclado externo e um processador que nao seja compativel com i386 ou >>> x86_64, entao (hoje!) eh melhor vc usar um laptop ou, se vc estiver >>> mesmo preocupado com tamanho/peso, use um netbook que custa 1/3 do >>> valor e vai fazer exatamente o q tablet faria em termos de R... >>> >>> Se voce estiver nos EUA, o Asus Transformer Prime vai te custar, pelo >>> menos, USD 499 (tablet) + USD 149 (dock) = USD 648.... E, por esse >>> valor, vc compra algum Toshiba com processador i5 e 6GB RAM.... Vai >>> ter ate' quem prefira comprar um Dell com processador i3 com os 499 e >>> um netbook por alguns 200 dolares... Ai' vai do custo-beneficio de >>> cada um... >>> >>> Se for para eu produzir algo de util sem um teclado externo, eu nao >>> consigo ficar mais de 10min usando um tablet. Se for para usar um >>> tablet com teclado externo, prefiro um (net/note)book. (preferencia >>> pessoal: mais aplicativos disponiveis, melhor suporte de >>> software/hardware e por ai' vai) >>> >>> O mesmo aplica-se (novamente, *pelo menos por enquanto*) ao caso do >>> Ubuntu com seu desktop para Android: se vc nao tiver um teclado >>> externo e um monitor, continua "inutil" para nossas aplicacoes. >>> >>> Concordo que, cada vez mais, teremos `a disposicao equipamentos mais >>> portateis. Entretanto, para "producao" ainda ficaremos dependentes de >>> teclado e bateria por mais alguns (poucos, espero) anos... >>> >>> Enquanto isso, todo o exercicio de rodar R (ou equivalentes de outros >>> setores, etc) sob Android ou iOS vale (como ja' dito pelo Leonard) >>> pela prova do conceito e inovacoes que ainda virao disso... Mas, ainda >>> nao para "producao"... not yet.. >>> >>> Comentario provocativo: os comentarios acima tambem valem para os >>> "beneficios" da programacao para GPU vis-a-vis prog. para CPU... >>> >>> b >>> >>> >>> 2012/3/1 Rodrigo Sant'Ana <rodrigo.gringo(a)gmail.com>: >>>> Pô Benilton concordo com você em alguns pontos. Mas acho que podemos >>>> rodar algumas rotinas em tablets sim. Esse Asus vem com um dock de >>>> teclado bastante interessante, e talvez, para algumas matérias básicas >>>> e scripts simples, como mesmo calcular uma média, desvio padrão, >>>> tabela de frequencia, histogramas, teste t, e outros tipos de >>>> inferências mais comuns e simples porém didáticas e necessárias para >>>> uma sala de aula, talvez não no nível de aprofundamento de um curso de >>>> graduação mais direcionado, mas em matérias inicias, ou em mini-cursos >>>> com bases de dados sucintas utilizadas para exemplificação, seja >>>> interessante. >>>> >>>> Entendo que com grandes banco de dados não seja o caminha, mas >>>> utilizar um tablet para acesso remoto via ssh, onde o trabalho pesado >>>> será executado por um computador mais potente nos dê mobilidade, é um >>>> tecnologia interessante. Acho que a explanação do Leonardo é válida, >>>> talvez teremos celulares como computadores pessoais, realizando >>>> trabalho em conjunto com computadores mais potentes.. >>>> >>>> Acho bastante válido, mas entendo sua colocação Benilton, com certeza, >>>> no dia a dia de analises mais densas, não seja tão interessante. >>>> >>>> Abração >>>> _______________ >>>> Rodrigo Sant'Ana >>>> Oceanógrafo >>>> >>>> >>>> >>>> Em 1 de março de 2012 15:38, Leandro Marino >>>> <leandromarino(a)leandromarino.com.br> escreveu: >>>>> Para rodar o R via ssh você terá ter ter uma conexão ativa a internet o que >>>>> muitas vezes não é bom. >>>>> >>>>> Alguns softwares, como por exemplo, o QuantumGis já estão com versões para >>>>> Android. Acho que é um caminho sim. >>>>> >>>>> Cabe salientar que a Canonical (desenvolvedora do Ubuntu) já anunciou versão >>>>> do Ubuntu para Celulares e Tablets. Inclusive é possível ter um telefone da >>>>> motorola que tenha webtop como Atrix rodando o ubuntu. >>>>> >>>>> Com o desenvolvimento do Windows 8 suportando processadores ARM, e linux, >>>>> acredito que seja uma futura convergencia. Acredito que em breve os >>>>> celulares serão nossos computadores pessoais... :) >>>>> >>>>> >>>>> >>>>> Atenciosamente, >>>>> Leandro Marino >>>>> http://www.leandromarino.com.br (Fotógrafo) >>>>> http://est.leandromarino.com.br/Blog (Estatístico) >>>>> Cel.: + 55 21 9845-7707 >>>>> Cel.: + 55 21 8777-7907 >>>>> >>>>> >>>>> >>>>> Em 1 de março de 2012 15:22, Rodrigo Sant'Ana <rodrigo.gringo(a)gmail.com> >>>>> escreveu: >>>>> >>>>>> Depende Leonard, >>>>>> >>>>>> existem tablets hoje em dia bastante robustos, um exemplo é o Asus >>>>>> Transformer Prime (sem querer fazer propaganda. Que vem com um >>>>>> processador Quad-core da nVidia, num é a solução, nem se compara a um >>>>>> PC. Mas ter o R rodando nele, para rotinas rápidas, talvez seja >>>>>> funcional para quem trabalha viajando muito! >>>>>> >>>>>> Só quis colocar a ideia, pois lembrei destas novas gerações de >>>>>> tablets, e até mesmo, de uma necessidade que eu compartilho, pois acho >>>>>> interessante ter esta mobilidade no uso do R, em um tablet. Para aulas >>>>>> em geral, para demonstrações de códigos em palestras, para uso em >>>>>> qualquer ocasião. >>>>>> >>>>>> Outra coisa, achei bem interessante essa parte do acesso ao R via ssh, >>>>>> talvez seja um tópico interessante para iniciar a discussão aqui na >>>>>> lista, com demonstrações e tudo mais! Pois eu desconheço completamente >>>>>> e acho super interessante e importante. >>>>>> >>>>>> Abraço a todos! >>>>>> >>>>>> _______________ >>>>>> Rodrigo Sant'Ana >>>>>> Oceanógrafo >>>>>> >>>>>> >>>>>> >>>>>> Em 29 de fevereiro de 2012 13:06, Benilton Carvalho >>>>>> <beniltoncarvalho(a)gmail.com> escreveu: >>>>>> > Concordo contigo, Leonard... mas como o OP nao perguntou minha >>>>>> > opiniao, ne'? rsrsrs >>>>>> > _______________________________________________ >>>>>> > R-br mailing list >>>>>> > R-br(a)listas.c3sl.ufpr.br >>>>>> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>>> > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>>> > código mínimo reproduzível. >>>>>> _______________________________________________ >>>>>> R-br mailing list >>>>>> R-br(a)listas.c3sl.ufpr.br >>>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>>> código mínimo reproduzível. >>>>> >>>>> >>>>> >>>>> _______________________________________________ >>>>> R-br mailing list >>>>> R-br(a)listas.c3sl.ufpr.br >>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código >>>>> mínimo reproduzível. >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)listas.c3sl.ufpr.br >>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)listas.c3sl.ufpr.br >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

3 3

Agrupamento de Tocher
by Melissa Oda 01 Mar '12

01 Mar '12

Alguém conhece um função no R para o agrupamento de Tocher? Obrigada Melissa

1 0

[Dúvida] Uso de linguagem funcional para estatística.
by Pedro Rafael 29 Fev '12

29 Fev '12

Alguem aqui domina alguma linguagem de programação funcional? Existe alguma vantagem de quem é da estatística ou matemática também aprender profundamente alguma linguagem funcional como é o caso de Haskell? Atualmente tudo que preciso consigo em R e C mas ultimamente venho lendo sobre linguagem funcional e cálculo lambda e me parece bem interessante para quem trabalha matemática.... -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

7 22

ebooks com 50% desconto
by Benilton Carvalho 29 Fev '12

29 Fev '12

Celebrando a Strata Conference deste ano, a O'Reilly esta' oferecendo 50% de desconto nos livros relacionados a dados (e analise). Divirtam-se: http://oreilly.com/data/index.html b

2 1