- R-br - DINF UFPR

Re: [R-br] Digest R-br, volume 8, assunto 33
by Danilo Scorzoni Ré 26 Out '11

26 Out '11

Obrigado pelas respostas, vocês são muitos solícitos. Já tinha usado a função glht() para dados com distribuição normal, usarei então junto com a função glm(). Um abraço à todos, Danilo Scorzoni Ré 2011/10/26 <r-br-request(a)listas.c3sl.ufpr.br> > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Gráficos no R ou Gnu Plot (Fernando Colugnati) > 2. Comparações Múltiplas em GLM (Danilo Scorzoni Ré) > 3. Re: Comparações Múltiplas em GLM (Edson Lira) > 4. Re: Comparações Múltiplas em GLM (Benilton Carvalho) > 5. Re: Comparações Múltiplas em GLM (Walmes Zeviani) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 25 Oct 2011 23:09:25 -0200 > From: Fernando Colugnati <fernando(a)ipti.org.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Gráficos no R ou Gnu Plot > Message-ID: > <CAOLzYt28ZB+K4we38999v_62mXJkK08inXdptTYX1SFzere+aA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Não devemos perpetuar práticas mesmo, e na estatística temos muitas > práticas > ruins sendo perpetuadas que não dizem respeito a como editar ou fazer um > gráfico para que mais pessoas utilizem. Temos a horrível prática do uso do > p-valor como oráculo, de testes de hipótese que testam hipóteses erradas, > modelos enviesados, etc, etc...com isso perco o sono, fico angustiado > etc...(e revistas de áreas usuárias da estatística pedem essas coisas, > assim > como pedem arquivos em xls...isso sim é o fim da picada). Com os gráficos > em > excel eu já realmente relaxei e passei a ter mais tempo pra estas outras > questões! Há pouco tempo um cara morreu de cancer tentando mudar o > mercado...e não conseguiu um grande sucesso neste quesito (mas em outros, > temos de concordar)! > Abs > > > Em 25 de outubro de 2011 21:29, Eder David Borges da Silva < > eder(a)leg.ufpr.br > > escreveu: > > > Pessoal, > > Este tópico é interessante é assim como o Pedro sei como é difícil > > utilizar um software que só você usar (R) dentro de uma empresa que só > > utilizar pacote office (E na maioria das vezes não desfrutando das > > poucas facilidades que eles te ajudam), tenho seguido a seguinte > > filosofia, não gostou do gráfico faça você!! com seu pacote preferido. > > Em minha análise eu creio que as pessoas de modo geral perder 50% do > > seu tempo "enfrufruzando" tabelas, gráficos e mapas de um relatório, > > acho isso um desperdício de vida, mas fazer oque né.... > > ate mais > > > > > > > > Em 25 de outubro de 2011 19:38, Paulo Justiniano > > <paulojus(a)leg.ufpr.br> escreveu: > > > perpetuar a pratica nao ajuda muito.... > > > > > > > > > On Tue, 25 Oct 2011, Pedro Rafael wrote: > > > > > >> Não Benilton é usar gráficos mais atraentes para ser colocado em um > > >> Power Point, não gráficos cegos e apagados para não dar sono no > > >> pessoal rsrsrsrs. Eu particularmente acho frescura. Sofro quando passo > > >> os gráficos feito no R porque tem algum fulano que quer editar tais > > >> gráficos. Criei várias funções que leem arquivos DBF que é formato que > > >> eu tenho aqui para calcular indicadores para estados, municípios, > > >> regionais de saúde mas os gráficos o pessoal reclama porque não sabem > > >> de R. Aí eu fico no dilema do que usar. Mas o Lema não é esse que você > > >> citou! Fazer um gráfico com um efeito em 3D a depender do gráfico não > > >> complica a cabeça de ninguém. > > >> > > >> 2011/10/25 beniltoncarvalho [via R-br] > > >> <ml-node+s2285057n3936773h40(a)n4.nabble.com>: > > >>> > > >>> Bem interessante a atitude da instituição: "para falar com leigos, > > >>> precisamos usar gráficos que complicam o entendimento da nossa > > mensagem". > > >>> > > >>> > > >>> -- > > >>> Successful people ask better questions, and as a result, they get > > better > > >>> answers. (Tony Robbins) > > >>> > > >>> _______________________________________________ > > >>> R-br mailing list > > >>> [hidden email] > > >>> 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. > > >>> > > >>> ________________________________ > > >>> If you reply to this email, your message will be added to the > > discussion > > >>> below: > > >>> > > >>> > > > http://r-br.2285057.n4.nabble.com/R-br-Graficos-no-R-ou-Gnu-Plot-tp3936555p… > > >>> To unsubscribe from R-br, click here. > > >> > > >> > > >> > > >> -- > > >> Saudações, > > >> Pedro Rafael Diniz Marinho. > > >> Estatístico - Secretaria de Estado da Saúde - PB. > > >> _______________________________________________ > > >> 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. > > > > > > -- > Fernando A.B. Colugnati > Pesquisador Associado > > Instituto de Pesquisas em Tecnologia e Inovação - IPTI > Tel. 55 11 8704-9812 > www.ipti.org.br > fernando(a)ipti.org.br > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111025/5a531e18/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 26 Oct 2011 10:11:06 -0200 > From: Danilo Scorzoni Ré <danilo.scorzoni(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Comparações Múltiplas em GLM > Message-ID: > <CA+bJVreC=uKJ_XJs4495NcK_tUuJ8FyQ+sWG5ruKhD77H+Mo3g(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia a todos, > > É a primeira vez que envio para essa lista e, ainda por cima, é um assunto > que eu já vi comentarem por aqui, mas não tenho mais o e-mail (e não > encontrei nenhum local na web que dê pra acessar as mensagens já > enviadas...), mas meu problema é o seguinte. > > Estou analisando um conjunto de dados de um delineamento casualizado em > blocos em duas áreas distintas. A primeira questão é que esse experimento > está instalado com o mesmo delineamento em duas áreas distintas e gostaria > de avaliar o efeito da área. Segundo o livro do Pimentel-Gomes, eu preciso > checar a homogeneidade de variâncias antes de proceder a análise usual, > usei > então um teste de Bartlett e um de Levenne. Em ambos, a hipótese nula não > foi rejeitada, ou seja, não tenho evidência de que as variâncias sejam > heterogêneas. Porém, dada a natureza dos dados, não acho que seja isso > (proporções de variâncias entre a menor e a maior estão acima de 7). Dessa > forma, procedi a análise através do GLM com resposta com distribuição Gama. > Consegui um bom ajuste com boa análise de resíduos. A minha dúvida agora é: > como fazer comparações múltiplas usando o GLM? Sei que o Tukey pressupõe > modelos paramétricos, mas não me recordo se ele precisa pressupor > distribuição normal. > > Desde já agradeço a atenção de todos e desculpe pelo texto extenso. > > Um abraço, > > -- > Danilo Scorzoni Ré > Engenheiro Florestal > Mestrando em Ciência Florestal > FCA / UNESP - Botucatu > (14) 8180-2494 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/54381e6a/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 26 Oct 2011 05:33:17 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Comparações Múltiplas em GLM > Message-ID: > <1319632397.14122.YahooMailNeo(a)web160707.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Danilo, caso o seus dados não atenda as pressuposições,tem um teste não > paramétrico que faz teste de médias. Veja o exemplo abaixo do pacote > agricolae. > > Não esqueça de instalar o pacote agricolae > > library(agricolae) > > ?kruskal > > data(corn) > attach(corn) > str(corn) > comparison<-kruskal(observation,method,group=TRUE, main="corn") > detach(corn) > > No caso da homogeneidade da variância no pacote MASS tem a transformação de > BOXCOX que é muito usada. > > #Examples > ?boxcox > require(MASS) > boxcox(Volume ~ log(Height) + log(Girth), data = trees, > lambda = seq(-0.25, 0.25, length = 10)) > > boxcox(Days+1 ~ Eth*Sex*Age*Lrn, data = quine, > lambda = seq(-0.05, 0.45, len = 20)) > > > []'s. > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: Danilo Scorzoni Ré <danilo.scorzoni(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Quarta-feira, 26 de Outubro de 2011 8:11 > Assunto: [R-br] Comparações Múltiplas em GLM > > > Bom dia a todos, > > É a primeira vez que envio para essa lista e, ainda por cima, é um assunto > que eu já vi comentarem por aqui, mas não tenho mais o e-mail (e não > encontrei nenhum local na web que dê pra acessar as mensagens já > enviadas...), mas meu problema é o seguinte. > > Estou analisando um conjunto de dados de um delineamento casualizado em > blocos em duas áreas distintas. A primeira questão é que esse experimento > está instalado com o mesmo delineamento em duas áreas distintas e gostaria > de avaliar o efeito da área. Segundo o livro do Pimentel-Gomes, eu preciso > checar a homogeneidade de variâncias antes de proceder a análise usual, usei > então um teste de Bartlett e um de Levenne. Em ambos, a hipótese nula não > foi rejeitada, ou seja, não tenho evidência de que as variâncias sejam > heterogêneas. Porém, dada a natureza dos dados, não acho que seja isso > (proporções de variâncias entre a menor e a maior estão acima de 7). Dessa > forma, procedi a análise através do GLM com resposta com distribuição Gama. > Consegui um bom ajuste com boa análise de resíduos. A minha dúvida agora é: > como fazer comparações múltiplas usando o GLM? Sei que o Tukey pressupõe > modelos paramétricos, mas não me recordo se > ele precisa pressupor distribuição normal. > > Desde já agradeço a atenção de todos e desculpe pelo texto extenso. > > Um abraço, > > -- > Danilo Scorzoni Ré > Engenheiro Florestal > Mestrando em Ciência Florestal > FCA / UNESP - Botucatu > (14) 8180-2494 > > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/dac3b546/atta… > > > > ------------------------------ > > Message: 4 > Date: Wed, 26 Oct 2011 13:38:20 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] Comparações Múltiplas em GLM > Message-ID: > <CAO-arWNG_7VnipqftsexdORqGrG0ZWjtt_FnKcqEgiwXZA5O5A(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Hipoteses de distribuicao para casos assim referem-se `a distribuicao > dos **estimadores**, nao dos dados... E minha (fraca) memoria nao se > recorda de nenhuma hipotese de distribuicao no teste de Tukey, mas > posso estar enganado.... > > > ------------------------------ > > Message: 5 > Date: Wed, 26 Oct 2011 10:45:07 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Comparações Múltiplas em GLM > Message-ID: > <CAFU=EkYkW5EsmUR8aKG5KTsQHa+eVLsW4wkNmc_9LkyScO2taw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Tukey supõe distribuição normal para os estimadores (como Benilton > assinalou) pois é baseado da distribuição da amplitude total studentizada. > Assintóticamente os estimadores nos GLMs terão distribuição normal. Então é > possível usar o procedimento. Detalhe é que como nos GLMs a variância não é > constante e é em geral dependente da média, as estimativas terão diferentes > erros padrões, e por isso assim também serão os contrastes entre elas. A > multcomp::glht() oferece recursos para isso, mas não aplica teste de Tukey, > faz os contrastes de Tukey (2 a 2) e corrige p-valor por diversos métodos > (opções da função). > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111026/13fcd070/atta… > > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 8, assunto 33 > **************************************** > -- Danilo Scorzoni Ré Engenheiro Florestal Mestrando em Ciência Florestal FCA / UNESP - Botucatu (14) 8180-2494

1 0

Comparações Múltiplas em GLM
by Danilo Scorzoni Ré 26 Out '11

26 Out '11

Bom dia a todos, É a primeira vez que envio para essa lista e, ainda por cima, é um assunto que eu já vi comentarem por aqui, mas não tenho mais o e-mail (e não encontrei nenhum local na web que dê pra acessar as mensagens já enviadas...), mas meu problema é o seguinte. Estou analisando um conjunto de dados de um delineamento casualizado em blocos em duas áreas distintas. A primeira questão é que esse experimento está instalado com o mesmo delineamento em duas áreas distintas e gostaria de avaliar o efeito da área. Segundo o livro do Pimentel-Gomes, eu preciso checar a homogeneidade de variâncias antes de proceder a análise usual, usei então um teste de Bartlett e um de Levenne. Em ambos, a hipótese nula não foi rejeitada, ou seja, não tenho evidência de que as variâncias sejam heterogêneas. Porém, dada a natureza dos dados, não acho que seja isso (proporções de variâncias entre a menor e a maior estão acima de 7). Dessa forma, procedi a análise através do GLM com resposta com distribuição Gama. Consegui um bom ajuste com boa análise de resíduos. A minha dúvida agora é: como fazer comparações múltiplas usando o GLM? Sei que o Tukey pressupõe modelos paramétricos, mas não me recordo se ele precisa pressupor distribuição normal. Desde já agradeço a atenção de todos e desculpe pelo texto extenso. Um abraço, -- Danilo Scorzoni Ré Engenheiro Florestal Mestrando em Ciência Florestal FCA / UNESP - Botucatu (14) 8180-2494

4 3

[OFF] Morre John McCarthy, pioneiro da Inteligência Artificial e pai do Lisp.
by Pedro Rafael 25 Out '11

25 Out '11

Morre John McCarthy, pioneiro da Inteligência Artificial e pai do Lisp. O matemático e pesquisador da área de ciências da computação John McCarthy, responsável pela criação do termo “inteligência artificial” e inventor da linguaguem de programação Lisp, morreu no último domingo (23) em Palo Alto, na região da Califórnia conhecida como Vale do Silício. http://g1.globo.com/tecnologia/noticia/2011/10/morre-john-mccarthy-pioneiro… -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Leitura arquivos .DB
by Eder David Borges da Silva 25 Out '11

25 Out '11

Bom dia pessoal, Gostaria de pedir ajudar a vocês na leitura de arquivos .DB (database), tentei funções do pacote RODBC, require(RODBC) con0 <- odbcConnectDbase('http://www.leg.ufpr.br/~eder/E2D00007.DB') con1 <- odbcConnectDbase('http://www.leg.ufpr.br/~eder/M00777600m000.DB') porem estou em windows 64 Bit, preciso ler neste sistema operacional, e da msn de erro pelo driver ser 32 bit e outras coisa s mais, desta forma alguém tem uma sugestão para leitura destes arquivos? Alternativa a essa? OBS: os arquivos estão no caminho indicado acima. Grato pela atenção

6 27

order
by Leonard de Assis 25 Out '11

25 Out '11

Boa noite a todos, Estava escrevendo um email gigantesco, detalhando um problema que eu estava enfrentando em uma simulação aqui, mas no meio do mesmo eu descobri um jeito de fazer o que pretendo, então sobrou pra mim o passo 2 da mesma: Conseguir um meio menos "marreta" de fazer. em meu caso, Menos "marreta" implica em algo mais rápido que o que já tenho, pois apesar de ser em um dataframe pequeno, tenho que fazer este passo algumas milhares de vezes em meu algorítmo. então vamos lá: eu preciso adicionar em um dataframe, uma coluna informando a ordem de cada um dos eventos que ocorreram. Eu atualmente faço da seguinte forma: G[order(-G$A,-G$B),]$Pos <- 1:length(G$V) do jeito que está, funciona, minha intenção é algo mais eficiente, se possível. -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

1 0

Regressão polinomial. Uma discussão teórica!
by Ivan Bezerra Allaman 24 Out '11

24 Out '11

Bom tarde senhores!! Me deparei com uma situação que em outrora não havia questionado tal situação, seguindo então o que a literatura indicava. A situação ao qual me refiro, é que ao desdobrar a soma de quadrados de regressão em termos lineares, quadráticos e cúbicos, os três termos foram significativos. FATORES GL SQ QM F value p value Silag:ordered(Tempo):FosForm11.L 1 16.26 16.26 33.97 0.0000 Silag:ordered(Tempo):FosForm11.Q 1 0.22 0.22 0.46 0.5000 Silag:ordered(Tempo):FosForm11.C 1 10.51 10.51 21.96 0.0000 Residuals 108 51.68 0.48 A literatura nos recomenda escolher o termo significativo de maior grau. Sinceramente, ao pararmos para pensar, porque não escolher o termo com a maior SQ, no caso de significância é claro dos dois termos? Afinal, o SQ nos informa o quanto da variação se deve aquele fator não? Esse pensamento, podemos ter no caso de um cálculo a mão certo! Mais hoje em dia, em que temos o valor do "p value", não seria mais prudente escolhermos o termo com o maior "p value"? que consequentemente, é claro, terá o maior SQ? Gostaria de saber a opinião dos colegas a respeito disso. Abraço a todos! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

6 7

Instalando versão mais nova do R no Ubuntu
by Pedro Rafael 24 Out '11

24 Out '11

Pessoal estou tentando instalar a versão mais recente do R (2.13.2) no Ubuntu 11.10 (64bist) e adicionei o repositório deb http://cran.stat.ucla.edu/bin/linux/ubuntu oneiric/ e em seguida dei um update. O seguinte erro apareceu na hora de atualizar o repositório. *W: Erro GPG: http://cran.stat.ucla.edu oneiric/ Release: As assinaturas a seguir não puderam ser verificadas devido à chave pública não estar disponível: NO_PUBKEY 51716619E084DAB9* Após esse erro tentei fazer como o site do R manda, fui atrás de validar a chave e fiz: *gpg --keyserver keyserver.ubuntu.com --recv-key E084DAB9* Após rodar o comando acima a falha continua. O que pode ser? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 17

ByVal, ByRef
by Leonard de Assis 23 Out '11

23 Out '11

Já existe no R alguma forma de se passar parâmetros uma função explicitando se estes são por valor ou por referência? -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

2 1

Erro ao instalar o pacote qtbase
by Luís Gustavo 23 Out '11

23 Out '11

Olá pessoal, estou tentando instalar o pacote "qtbase", seguindo os comandos no site http://qtinterfaces.r-forge.r-project.org/. Porém, quando chego na etapa de "R CMD INSTALL qtbase", vai instalando tudo ok, e já no final ocorre o erro logo abaixo. Lembrando que estou usando linux em uma maquina virtual. Como posso resolver este problema? Obrigado. Install the project... -- Install configuration: "" -- Up-to-date: /home/luis/qtbase/src/./qtbase.so make[1]: Saindo do diretório `/home/luis/qtbase/src-build' installing to /home/luis/R/i686-pc-linux-gnu-library/2.11/qtbase/libs ** R ** demo ** inst ** preparing package for lazy loading ** help *** installing help indices ** building package indices ... ** testing if installed package can be loaded Error : .onLoad failed in loadNamespace() for 'qtbase', details: call: FUN(c("operator!=", "operator&", "operator*", "operator+", "operator-", error: não foi possível encontrar a função "list2env" ERROR: loading failed * removing ‘/home/luis/R/i686-pc-linux-gnu-library/2.11/qtbase’ -- Luís Gustavo Silva e Silva

2 2

plotar erro padrão usando "arrows"
by Rodrigo Plei 21 Out '11

21 Out '11

Prezados, Como plotar o erro padrão em um gráfico de médias usando a função "arrows"? Agradeço desde logo a atenção. 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 6

Dúvida em Cluster!
by Elton Gean Araújo Araújo 21 Out '11

21 Out '11

Boa tarde prezados colegas da lista. É a primeira vez que trabalho com analise de cluster e me deparei com um problema. Eu quero agrupar os países segundo os indices de satisfação (Sat), mais quando rodo o gráfico, o eixo x não sai com os nomes dos países, mais sim a identificação só (1,2,3...). Alguém poderia me dizer o que está faltando para aparecer os nomes dos países? Segue o script que estou usando e os dados em anexo. hc <- hclust(dist(model.matrix(~-1 + Sat, dados)) , method= "single") plot(hc,hang=-0.1). Desde já, agradeço qualquer ajuda. Abraços, Elton.

3 3

substitui loop por funcoes do R
by Samuel Carvalho 21 Out '11

21 Out '11

Caros membros do grupo, Alguém poderia por gentileza me auxiliar com o código a seguir? No caso a idéia é substituir o loop que calcula a diferença entre as observações dentro de uma mesma parcela sendo que a primeira observação de cada parcela deve ter valor zerado. Para aqueles mais familiarizados estou calculando a taxa de incremento para sucessivas medições. Segue o código ###RCMR dados<-data.frame(parcela=rep(1:5,each=3), y=runif(15)) dados dados$taxa<-0 for(i in 1:(nrow(dados)-1)) { if(dados$parcela[i] == dados$parcela[i+1]) { dados$taxa[i+1] <- dados$y[i+1]-dados$y[i] } } dados #A seguir tento reproduzir o código acima, porém feito incorretamente dados$taxa<-0 dados$taxa[which(diff(dados$parcela)==0)] <- diff(dados$y) #feito de maneira incorreta dados Abraços Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 2

previsão de valores futuros no pacote nnet
by Alexandro (Yahoo) 20 Out '11

20 Out '11

#Boa tarde R users, supondo o CMR de uma rede neural do pacote nnet. Gostaria de saber: dado o vetor(x1,x2) >>> previsao=c(15,3) qual seria o resultado para a variável y. Porém o comando: #predict(rede,previsao=previsao, type = "class") #resulta num vetor em que não sei qual é a resposta certa ! y=c(0 ,1 ,0 ,0 ,1 ,1 ,1 ,0 ,0 ,0) x1=c(12,19,14,13,14,15,16,17,17,18) x2=c(3,7,8,6,2,9,5,1,8,7) amostra<-c(sample(1:5,5),sample(6:10,5));amostra dados=data.frame(yfactor=factor(y),x1,x2) require(nnet) rede<- nnet(yfactor~x1+x2, data = dados, subset = amostra, size = 2, rang = 0.1,decay = 5e-4, maxit = 200) matriz_confusao=table(y, predict(rede,dados, type = "class")) taxa_acertos= 100*( matriz_confusao[1,1]+ matriz_confusao[2,2] ) / sum(matriz_confusao) cat(" A taxa de acertos é de",taxa_acertos,"%","\n") previsao=c(15,3) predict(rede,previsao=previsao, type = "class") #Assim, qual meu erro na previsão de y dado o vetor previsão ? #Desde já, agradeço-lhes ! Felicidades ! #Alexandro

2 1

Incremento de Workpace
by Edson Lira 20 Out '11

20 Out '11

Bom dia a todos. Estou tentando fazer uma associação: CrossTable(b1$coinf_2,b1$Grupos_Final,fisher = T,prop.t=T,prop.r=T,prop.c=T,format=c("SPSS")) Que me dá a mensagem abaixo: Erro em fisher.test(t, alternative = "two.sided") : FEXACT error 7. LDSTP is too small for this problem. Try increasing the size of the workspace. Alguém poderia me dar alguma sugestão ou rotina para resolver esse problema? []'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Re: [R-br] Digest R-br, volume 8, assunto 22
by Cassiano 20 Out '11

20 Out '11

Re: Barra de erros usando Binomial Negativa (Walmes Zeviani) Walmes, Acho que é exatamente isso que estava procurando. Mas eu não entendi o output que o comando me ofereceu. Ou pelo menos, não sei usá-lo para construir meu gráfico. > confint(m2) Waiting for profiling to be done... 2.5 % 97.5 % (Intercept) 0.7804852 1.2279976 predadorsem -1.0985071 -0.3178279 Eu tenho as categorias com e sem predador. Ok, o Intercept é com o predador e a linha seguinte, sem o predador. Qual é o valor da minha barra superior e inferior da barra com predador, p.ex.? Como ploto isso em um barplot2? Obrigado! Cassiano Em 20 de outubro de 2011 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. ajuda ks.test (Rodrigo Plei) > 2. Re: ajuda ks.test (Gustavo Henrique M. A. Rocha) > 3. Re: ajuda ks.test (Rodrigo Plei) > 4. Re: ajuda ks.test (Gustavo Henrique M. A. Rocha) > 5. Barra de erros usando Binomial Negativa (Cassiano) > 6. Re: ajuda ks.test (Benilton Carvalho) > 7. Re: ajuda ks.test (Edson Lira) > 8. Re: ajuda ks.test (Rodrigo Plei) > 9. Re: Barra de erros usando Binomial Negativa (Walmes Zeviani) > 10. recuperar coeficientes da função discriminante quadrática > (Guilherme Moraes Ferraudo) > 11. Re: recuperar coeficientes da função discriminante quadrática > (Benilton Carvalho) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 19 Oct 2011 14:04:34 -0200 > From: Rodrigo Plei <rodrigo.plei(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] ajuda ks.test > Message-ID: > <CACs_BFNAu-TBO+q1eUQWG1e2yRGYnCgUmt6M7vuASiqBNAkgBQ(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Pessoal, > > Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote > "stats"), > para comparar frequências de comprimentos de peixes de duas áreas de pesca > distintas, mas algo está emperrando a análise: > > x <- data # frequência de comprimentos de peixes da área x > y <- data2 # frequência de comprimentos de peixes na área y > ks.test(x, y) > Erro em ks.test(x, y) : > 'y' deve ser numérico ou uma cadeia nomeando uma função válida > > O formado dos dados é: > > area_x area_y > 23 29 > 25 30 > 23.8 26 > 19 28 > > Será que estou entrando os dados de maneira errada? > > Agradeço qualquer ajuda 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/96e46770/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 19 Oct 2011 14:12:28 -0200 > From: "Gustavo Henrique M. A. Rocha" <ghmarbr(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] ajuda ks.test > Message-ID: > <CADsPq2q4g6W-abJ1gXWBbnUCqbd+MrzVkRUWC6iSag39v3k=2g(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Rodrigo, > > possivelmente seus dados não estão definidos como numéricos. > > Use a função is.numeric para verificar se os dados são numéricos ou não. > > Fazendo um teste com os dados que você mostrou consegui rodar o teste > normalmente. > > x = c(23,25,23.8,19) > y = c(29,30,26,28) > ks.test(x,y) > > Abraços! > > > Gustavo Henrique Mitraud Assis Rocha > Doutorando em Estatística pela Universidade Federal de Minas Gerais > Mestre em Estatística pela Universidade Federal de Minas Gerais > Bacharel em Estatística pela Universidade Federal de Minas Gerais > > > > Em 19 de outubro de 2011 14:04, Rodrigo Plei <rodrigo.plei(a)gmail.com> > escreveu: > > Pessoal, > > Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote > "stats"), > > para comparar frequências de comprimentos de peixes de duas áreas de > pesca > > distintas, mas algo está emperrando a análise: > > x <- data # frequência de comprimentos de peixes da área x > > y <- data2 # frequência de comprimentos de peixes na área y > > ks.test(x, y) > > Erro em ks.test(x, y) : > > 'y' deve ser numérico ou uma cadeia nomeando uma função válida > > O formado dos dados é: > > area_x area_y > > 23 29 > > 25 30 > > 23.8 26 > > 19 28 > > > > Será que estou entrando os dados de maneira errada? > > Agradeço qualquer ajuda 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 > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ > > 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. > > > > > ------------------------------ > > Message: 3 > Date: Wed, 19 Oct 2011 14:24:04 -0200 > From: Rodrigo Plei <rodrigo.plei(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] ajuda ks.test > Message-ID: > <CACs_BFN4zjFaY264+5GVONMhq1ZUC7kwyxzKaZqPtsTQNMYYGg(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Parece que é isso mesmo, Gustavo. > > Mas o que eu não entendo é que o comando "is.numeric" não está funcionando > para meus dados: > > > x1 <-is.numeric(x) > > x1 > [1] FALSE > > Alguma idéia? > > Obrigado desde logo. > > R. > > Em 19 de outubro de 2011 14:12, Gustavo Henrique M. A. Rocha < > ghmarbr(a)gmail.com> escreveu: > > > Rodrigo, > > > > possivelmente seus dados não estão definidos como numéricos. > > > > Use a função is.numeric para verificar se os dados são numéricos ou não. > > > > Fazendo um teste com os dados que você mostrou consegui rodar o teste > > normalmente. > > > > x = c(23,25,23.8,19) > > y = c(29,30,26,28) > > ks.test(x,y) > > > > Abraços! > > > > > > Gustavo Henrique Mitraud Assis Rocha > > Doutorando em Estatística pela Universidade Federal de Minas Gerais > > Mestre em Estatística pela Universidade Federal de Minas Gerais > > Bacharel em Estatística pela Universidade Federal de Minas Gerais > > > > > > > > Em 19 de outubro de 2011 14:04, Rodrigo Plei <rodrigo.plei(a)gmail.com> > > escreveu: > > > Pessoal, > > > Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote > > "stats"), > > > para comparar frequências de comprimentos de peixes de duas áreas de > > pesca > > > distintas, mas algo está emperrando a análise: > > > x <- data # frequência de comprimentos de peixes da área x > > > y <- data2 # frequência de comprimentos de peixes na área y > > > ks.test(x, y) > > > Erro em ks.test(x, y) : > > > 'y' deve ser numérico ou uma cadeia nomeando uma função válida > > > O formado dos dados é: > > > area_x area_y > > > 23 29 > > > 25 30 > > > 23.8 26 > > > 19 28 > > > > > > Será que estou entrando os dados de maneira errada? > > > Agradeço qualquer ajuda 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 > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ > > > 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. > > > > > > -- > =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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/f4fffa48/atta… > > > > ------------------------------ > > Message: 4 > Date: Wed, 19 Oct 2011 14:26:13 -0200 > From: "Gustavo Henrique M. A. Rocha" <ghmarbr(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] ajuda ks.test > Message-ID: > <CADsPq2qQ1uAoUiE5uFD6262e_CFWBYXTcKzF4G0eE80rpw-nsA(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Rodrigo, > > no caso o resultado diz que sua variável x não é numérica. > > Como não sei a estrutura de seus dados, tente usar > ks.test(as.numeric(x),as.numeric(y)) > > > Gustavo Henrique Mitraud Assis Rocha > Doutorando em Estatística pela Universidade Federal de Minas Gerais > Mestre em Estatística pela Universidade Federal de Minas Gerais > Bacharel em Estatística pela Universidade Federal de Minas Gerais > > > > Em 19 de outubro de 2011 14:24, Rodrigo Plei <rodrigo.plei(a)gmail.com> > escreveu: > > Parece que é isso mesmo, Gustavo. > > Mas o que eu não entendo é que o comando "is.numeric" não está > funcionando > > para meus dados: > >> x1 <-is.numeric(x) > >> x1 > > [1] FALSE > > Alguma idéia? > > Obrigado desde logo. > > R. > > Em 19 de outubro de 2011 14:12, Gustavo Henrique M. A. Rocha > > <ghmarbr(a)gmail.com> escreveu: > >> > >> Rodrigo, > >> > >> possivelmente seus dados não estão definidos como numéricos. > >> > >> Use a função is.numeric para verificar se os dados são numéricos ou não. > >> > >> Fazendo um teste com os dados que você mostrou consegui rodar o teste > >> normalmente. > >> > >> x = c(23,25,23.8,19) > >> y = c(29,30,26,28) > >> ks.test(x,y) > >> > >> Abraços! > >> > >> > >> Gustavo Henrique Mitraud Assis Rocha > >> Doutorando em Estatística pela Universidade Federal de Minas Gerais > >> Mestre em Estatística pela Universidade Federal de Minas Gerais > >> Bacharel em Estatística pela Universidade Federal de Minas Gerais > >> > >> > >> > >> Em 19 de outubro de 2011 14:04, Rodrigo Plei <rodrigo.plei(a)gmail.com> > >> escreveu: > >> > Pessoal, > >> > Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote > >> > "stats"), > >> > para comparar frequências de comprimentos de peixes de duas áreas de > >> > pesca > >> > distintas, mas algo está emperrando a análise: > >> > x <- data # frequência de comprimentos de peixes da área x > >> > y <- data2 # frequência de comprimentos de peixes na área y > >> > ks.test(x, y) > >> > Erro em ks.test(x, y) : > >> > 'y' deve ser numérico ou uma cadeia nomeando uma função válida > >> > O formado dos dados é: > >> > area_x area_y > >> > 23 29 > >> > 25 30 > >> > 23.8 26 > >> > 19 28 > >> > > >> > Será que estou entrando os dados de maneira errada? > >> > Agradeço qualquer ajuda 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 > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > _______________________________________________ > >> > 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. > > > > > > > > -- > > =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 > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ > > 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. > > > > > ------------------------------ > > Message: 5 > Date: Wed, 19 Oct 2011 14:27:28 -0200 > From: Cassiano <cassianosr(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Barra de erros usando Binomial Negativa > Message-ID: > <CAF78VjJQ8vpJkBFnXWzp02Mn8EMwe69sNFfRfA+HjosuHsVeWg(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Olá a todos, > > Tenho dados de contagem de indivíduos (no y) e a distribuição mais adequada > foi a Bionmial Negativa. > Como meu x são categorias, vou demostrar isso com um barplot. > Porém não sei como gerar as barras de erro assimétricas, uma vez que meus > dados não são Normais. > > Como eu devo calcular o desvio padrão nestes casos? > > Obrigado! > Cassiano > > -- > ======================================= > MSc Cassiano S. Rosa > Universidade Federal de Viçosa -UFV > Dept. Biologia Animal > Pós-graduação em Entomologia > Lab. Termitologia > > Curriculo lattes: > http://lattes.cnpq.br/3030825329622014 > ======================================= > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/0739ccfe/atta… > > > > ------------------------------ > > Message: 6 > Date: Wed, 19 Oct 2011 17:29:42 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] ajuda ks.test > Message-ID: > <CAO-arWPcFfQ=jD8Duvn1UQB0sLHTHn0LTF6kh2-TSb0yuTixzA(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > talvez vc queira "as.numeric" (assumindo q seja algo q possa ser > convertido para numeric). > > mas, o melhor de tudo mesmo seria um exemplo reproduzivel... tipo > (voltando ao seu exemplo inicial) > > dput(data) > dput(data2) > > veja tbm o topico para algumas sugestoes: > > > http://r-br.2285057.n4.nabble.com/R-br-Produzindo-um-Codigo-Minimo-Reproduz… > > http://goo.gl/wlNWC > > b > > 2011/10/19 Rodrigo Plei <rodrigo.plei(a)gmail.com>: > > Parece que é isso mesmo, Gustavo. > > Mas o que eu não entendo é que o comando "is.numeric" não está > funcionando > > para meus dados: > >> x1 <-is.numeric(x) > >> x1 > > [1] FALSE > > Alguma idéia? > > Obrigado desde logo. > > R. > > Em 19 de outubro de 2011 14:12, Gustavo Henrique M. A. Rocha > > <ghmarbr(a)gmail.com> escreveu: > >> > >> Rodrigo, > >> > >> possivelmente seus dados não estão definidos como numéricos. > >> > >> Use a função is.numeric para verificar se os dados são numéricos ou não. > >> > >> Fazendo um teste com os dados que você mostrou consegui rodar o teste > >> normalmente. > >> > >> x = c(23,25,23.8,19) > >> y = c(29,30,26,28) > >> ks.test(x,y) > >> > >> Abraços! > >> > >> > >> Gustavo Henrique Mitraud Assis Rocha > >> Doutorando em Estatística pela Universidade Federal de Minas Gerais > >> Mestre em Estatística pela Universidade Federal de Minas Gerais > >> Bacharel em Estatística pela Universidade Federal de Minas Gerais > >> > >> > >> > >> Em 19 de outubro de 2011 14:04, Rodrigo Plei <rodrigo.plei(a)gmail.com> > >> escreveu: > >> > Pessoal, > >> > Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote > >> > "stats"), > >> > para comparar frequências de comprimentos de peixes de duas áreas de > >> > pesca > >> > distintas, mas algo está emperrando a análise: > >> > x <- data # frequência de comprimentos de peixes da área x > >> > y <- data2 # frequência de comprimentos de peixes na área y > >> > ks.test(x, y) > >> > Erro em ks.test(x, y) : > >> > 'y' deve ser numérico ou uma cadeia nomeando uma função válida > >> > O formado dos dados é: > >> > area_x area_y > >> > 23 29 > >> > 25 30 > >> > 23.8 26 > >> > 19 28 > >> > > >> > Será que estou entrando os dados de maneira errada? > >> > Agradeço qualquer ajuda 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 > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > _______________________________________________ > >> > 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. > > > > > > > > -- > > =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 > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ > > 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. > > > > > > -- > Successful people ask better questions, and as a result, they get > better answers. (Tony Robbins) > > > ------------------------------ > > Message: 7 > Date: Wed, 19 Oct 2011 10:01:01 -0700 (PDT) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] ajuda ks.test > Message-ID: > <1319043661.15905.YahooMailNeo(a)web160713.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Se já resolveu, ótimo, se não é bom você usar a função str(data) que te > dará qual a estrutura de todas as variáveis de seu banco de dados. > > Depois disso, pode-se usar as.numeric, caso alguma variável não seja > numérica. Mas antes cheque o seu banco, pode ser que ao digitar tenha havido > alguns erros: por exemplo, pode ser que tenha um ponto no lugar de uma > vírgula, ou outro problema similar. > > > []'s. > > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Quarta-feira, 19 de Outubro de 2011 12:29 > Assunto: Re: [R-br] ajuda ks.test > > talvez vc queira "as.numeric" (assumindo q seja algo q possa ser > convertido para numeric). > > mas, o melhor de tudo mesmo seria um exemplo reproduzivel... tipo > (voltando ao seu exemplo inicial) > > dput(data) > dput(data2) > > veja tbm o topico para algumas sugestoes: > > > http://r-br.2285057.n4.nabble.com/R-br-Produzindo-um-Codigo-Minimo-Reproduz… > > http://goo.gl/wlNWC > > b > > 2011/10/19 Rodrigo Plei <rodrigo.plei(a)gmail.com>: > > Parece que é isso mesmo, Gustavo. > > Mas o que eu não entendo é que o comando "is.numeric" não está > funcionando > > para meus dados: > >> x1 <-is.numeric(x) > >> x1 > > [1] FALSE > > Alguma idéia? > > Obrigado desde logo. > > R. > > Em 19 de outubro de 2011 14:12, Gustavo Henrique M. A. Rocha > > <ghmarbr(a)gmail.com> escreveu: > >> > >> Rodrigo, > >> > >> possivelmente seus dados não estão definidos como numéricos. > >> > >> Use a função is.numeric para verificar se os dados são numéricos ou não. > >> > >> Fazendo um teste com os dados que você mostrou consegui rodar o teste > >> normalmente. > >> > >> x = c(23,25,23.8,19) > >> y = c(29,30,26,28) > >> ks.test(x,y) > >> > >> Abraços! > >> > >> > >> Gustavo Henrique Mitraud Assis Rocha > >> Doutorando em Estatística pela Universidade Federal de Minas Gerais > >> Mestre em Estatística pela Universidade Federal de Minas Gerais > >> Bacharel em Estatística pela Universidade Federal de Minas Gerais > >> > >> > >> > >> Em 19 de outubro de 2011 14:04, Rodrigo Plei <rodrigo.plei(a)gmail.com> > >> escreveu: > >> > Pessoal, > >> > Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote > >> > "stats"), > >> > para comparar frequências de comprimentos de peixes de duas áreas de > >> > pesca > >> > distintas, mas algo está emperrando a análise: > >> > x <- data # frequência de comprimentos de peixes da área x > >> > y <- data2 # frequência de comprimentos de peixes na área y > >> > ks.test(x, y) > >> > Erro em ks.test(x, y) : > >> > 'y' deve ser numérico ou uma cadeia nomeando uma função válida > >> > O formado dos dados é: > >> > area_x area_y > >> > 23 29 > >> > 25 30 > >> > 23.8 26 > >> > 19 28 > >> > > >> > Será que estou entrando os dados de maneira errada? > >> > Agradeço qualquer ajuda 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 > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > > >> > _______________________________________________ > >> > 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. > > > > > > > > -- > > =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 > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ > > 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. > > > > > > -- > Successful people ask better questions, and as a result, they get > better answers. (Tony Robbins) > _______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/f9a23689/atta… > > > > ------------------------------ > > Message: 8 > Date: Wed, 19 Oct 2011 15:33:45 -0200 > From: Rodrigo Plei <rodrigo.plei(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> > Subject: Re: [R-br] ajuda ks.test > Message-ID: > <CACs_BFOF_qeBoNNO9K8KJOxN01D95JfCCuKpGH244PFw_3+GMA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Prezados todos, > > Consegui resolver o problema. > > Bastou dizer ao R que o separador decimal dos dados originais erma vírgulas > e não pontos. > > 1 abraço e obrigado pela ajuda! > > 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/0488bd6e/atta… > > > > ------------------------------ > > Message: 9 > Date: Wed, 19 Oct 2011 15:41:16 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Barra de erros usando Binomial Negativa > Message-ID: > <CAFU=EkY7mwOG-wJrVUL=md5nSAMdndrhGt_7GQR5pvjrbdr_3w(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Pode-se usar um intervalo baseado em perfil de verossimilhança. Para > objetos > de classe lm, glm e nls, tem-se a função confint() que faz isso. Veja se > serve. > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/a62ea65f/atta… > > > > ------------------------------ > > Message: 10 > Date: Wed, 19 Oct 2011 18:49:49 -0200 > From: Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] recuperar coeficientes da função discriminante > quadrática > Message-ID: > <CA++u5Oh9FOsU=ezJeqJ4NKEMWtxxGPB910VhckLfDNATrZS+7A(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa noite, > > por favor, preciso dos coeficientes da função discriminante quadrática. > Estou usando o exemplo da função "qda" para tentar enteder a função. > > #código R > tr <- sample(1:50, 25) > train <- rbind(iris3[tr,,1], iris3[tr,,2], iris3[tr,,3]) > test <- rbind(iris3[-tr,,1], iris3[-tr,,2], iris3[-tr,,3]) > cl <- factor(c(rep("s",25), rep("c",25), rep("v",25))) > z <- qda(train, cl) > predict(z,test)$class > > > Desde já, muito obrigado. > > Att > > -- > Guilherme Moraes Ferraudo > Jaboticabal/Campinas - SP > http://lattes.cnpq.br/2096118558794430 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111019/c0c3cc72/atta… > > > > ------------------------------ > > Message: 11 > Date: Wed, 19 Oct 2011 22:05:09 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] recuperar coeficientes da função discriminante > quadrática > Message-ID: > <CAO-arWNp8_h=WEDx9dsjFeOMd-4=dwSnjfyuBeN4zGGvbuFE6Q(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > nao existe tal conceito em qda. > > vc tem as medias, matrizes de variancia-covariacia e prioris de cada > grupo... uma vez com elas, vc determina a posteriori de cada > observacao e classifica-a no grupo que der a probabilidade maxima. > > b > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 8, assunto 22 > **************************************** > -- ======================================= MSc Cassiano S. Rosa Universidade Federal de Viçosa -UFV Dept. Biologia Animal Pós-graduação em Entomologia Lab. Termitologia Curriculo lattes: http://lattes.cnpq.br/3030825329622014 =======================================

1 0

recuperar coeficientes da função discriminante quadrática
by Guilherme Moraes Ferraudo 19 Out '11

19 Out '11

Boa noite, por favor, preciso dos coeficientes da função discriminante quadrática. Estou usando o exemplo da função "qda" para tentar enteder a função. #código R tr <- sample(1:50, 25) train <- rbind(iris3[tr,,1], iris3[tr,,2], iris3[tr,,3]) test <- rbind(iris3[-tr,,1], iris3[-tr,,2], iris3[-tr,,3]) cl <- factor(c(rep("s",25), rep("c",25), rep("v",25))) z <- qda(train, cl) predict(z,test)$class Desde já, muito obrigado. Att -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430

2 1

Barra de erros usando Binomial Negativa
by Cassiano 19 Out '11

19 Out '11

Olá a todos, Tenho dados de contagem de indivíduos (no y) e a distribuição mais adequada foi a Bionmial Negativa. Como meu x são categorias, vou demostrar isso com um barplot. Porém não sei como gerar as barras de erro assimétricas, uma vez que meus dados não são Normais. Como eu devo calcular o desvio padrão nestes casos? Obrigado! Cassiano -- ======================================= MSc Cassiano S. Rosa Universidade Federal de Viçosa -UFV Dept. Biologia Animal Pós-graduação em Entomologia Lab. Termitologia Curriculo lattes: http://lattes.cnpq.br/3030825329622014 =======================================

2 1

ajuda ks.test
by Rodrigo Plei 19 Out '11

19 Out '11

Pessoal, Estou tentando utilizar o kolgomorov-smirnov test (ks test, pacote "stats"), para comparar frequências de comprimentos de peixes de duas áreas de pesca distintas, mas algo está emperrando a análise: x <- data # frequência de comprimentos de peixes da área x y <- data2 # frequência de comprimentos de peixes na área y ks.test(x, y) Erro em ks.test(x, y) : 'y' deve ser numérico ou uma cadeia nomeando uma função válida O formado dos dados é: area_x area_y 23 29 25 30 23.8 26 19 28 Será que estou entrando os dados de maneira errada? Agradeço qualquer ajuda 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

4 6

especificar coluna na leitura de arquivo TXT
by Mauricio Morgado Oliveira 18 Out '11

18 Out '11

Caros colegas, tem como especificar no *read.csv* a leitura de cada variável conforme o inicio e fim da coluna? Assim como no sas p.e.: caracteristica 4 nivel1 8 nivel2 18 medida 24-34 ep 41-50 1 1 1 -0.06899554 0.13888397 2 1 1 -0.08143789 0.17738412 1 1 2 0.05868521 0.13935269 pois tenho que ler um arquivo gerado por outro programa e vem formatado por diferentes números espaços entre colunas... Desde já agradeço Att, Mauricio

2 2

Re: [R-br] Identidade de modelos de regresão via lm
by Ivan Bezerra Allaman 18 Out '11

18 Out '11

Via nls que não altera as estimativas dos parâmetros você pode fazer: mg=nls(r ~ a[t] + b[t]*x,data=d,start=list(a=c(1,1),b=c(1,1))) ma=nls(r ~ a + b[t]*x,data=d,start=list(a=1,b=c(1,1))) anova(mg,ma) mb=nls(r ~ a[t] + b*x,data=d,start=list(a=c(1,1),b=1)) anova(mg,mb) (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

Identidade de modelos de regresão via lm
by Marcelo Cardoso mello 18 Out '11

18 Out '11

Como posso comparar os dois modelo (ma e mb) na CMR abaixo? Tipo um teste de identidade de modelos em lm? Como testar se os tratamentos 'a' e 'b' tem (estatisticamente) a mesma curva de regressão? Obrigado por qualquer dica. a<- 'a'; b<- 'b' t<-c(a,a,a,a,a,a,a,a,b,b,b,b,b,b,b,b) x<-c(0,0,10,10,30,30,45,45,0,0,10,10,30,30,45,45) r<-c(2,3,4,5,7,9,11,13,5,4,7,8,8,9,10,11) d<-data.frame(t,x,r);d da<-subset(d, t=='a'); da db<-subset(d, t=='b'); db m<-lm(r~x, data=d) ma<-lm(r~x, data=da) mb<-lm(r~x, data=db) m;ma;mb deviance(m)/df.residual(m) deviance(ma)/df.residual(ma) deviance(mb)/df.residual(mb)

2 1

Erro em Script SEM
by Joao Ricardo Nickenig Vissoci 18 Out '11

18 Out '11

Pessoal, boa noite. Sou recente no fórum e também no estudo com o R, mas estou aprendendo. Bom, queria tirar uma dúvida, estou trabalhando em um script para analise de Modelos de Equações Estruturais, e sempre que tento rodar o script que fiz encontro o seguinte erro: Iteration limit exceeded. O Script é esse, desculpem-me pelo Inglês, mas a discussão do grupo é toda em inglês. # MODEL II # Install the package 'sem' # carry data library library (sem) # Change R's directory to that in which the data base .CSV file is # Load database teste<- read.csv("BaseSEM.csv",header=T) teste # sem package model specification mod.wh.1 <- specify.model()# Type these values that specify the model's relations (just use de Ctrl+R over each relation). Support->SP,NA,1 Support->SM,sm1,NA Rejection->RP,rp1,NA Rejection->RM,NA,1 Overprot->OP,op1,NA Overprot->OM,NA,1 MotI->MIAO,NA,1 MotI->MIEE,moti1,NA MotI->MIC,moti2,NA MotE->MERE,NA,1 MotE->MEIn,mote1,NA MotE->MEId,mote2,NA Support->MotI,NA,1 Rejection->MotI,rej1,NA Overprot->MotI,over1,NA Support->MotE,sup,NA Rejection->MotE,rej2,NA Overprot->MotE,NA,1 MotI->SCI,misci,NA MotE->SCI,mesci,NA Support<->Support,ss1,NA Rejection<->Rejection,rr1,NA Overprot<->Overprot,oo1,NA OP<->OP,opcov,NA OM<->OM,omcov,NA SP<->SP,spcov,NA SM<->SM,smcov,NA RP<->RP,rpcov,NA RM<->RM,rmcov,NA MotI<->MotI,micov,NA MotE<->MotE,mecov,NA MIAO<->MIAO,miaocov,NA MIEE<->MIEE,mieecov,NA MIC<->MIC,miccov,NA MEIn<->MEIn,meincov,NA MEId<->MEId,meidcov,NA MERE<->MERE,merecov,NA SCI<->SCI,scicov,NA Support<->Rejection,suprej,NA Support<->Overprot,supoover,NA Rejection<->Overprot,rejover,NA MotI<->MotE,mots1,NA # Insert de covariance matrix cov.teste <- cov(teste, y = NULL, use = "everything", method = c("pearson", "kendall", "spearman")) cov.teste # Estimate de model (Here is where we have been finding difficulties) sem.wh.1 <- sem(mod.wh.1, cov.teste, N=126) summary(sem.wh.1) Obrigado!!! João Ricardo -- *Joao Ricardo N. Vissoci* Psicólogo CRP 08/12469 Prof. Ms. Faculdade Ingá Doutorando em Psicologia Social - PUCsp Grupo Pro-Esporte UEM/CNPq Tel. 4499298078/4399368588 joaovissoci@gmail.com/jrvissoci@ig.com.br proesporteuem.blogspot.com.br

2 2

Associação conjunta
by Edson Lira 18 Out '11

18 Out '11

Boa noite a todos, tenho um banco de dados com algumas variávei: AR(sim, não), AIN(sim, não, talvez), T_HIV(menos de 250, entre de 250 e 300, acima de 300), HISTO_A(normal, nia, nia I). Gostaria de avaliar a associação conjunta da variável AR com o as demais variáveis conjuntamente. Seria um modelo log-linear? Alguém poderia me dar alguma sugestão ou material para eu dar uma olhada? []'s Edson Lira Estatístico Manaus-Amazonas

3 3

Como obter a estrutura de um data.frame (linux)
by Jose Claudio Faria 18 Out '11

18 Out '11

Pessoal, Alguém poderia me informar uma função (no Linux) para exportar (por exemplo para a área de transferência) a estrutura de um data.frame. No Windows existe uma função (trExport) que acompanha o pacote TinnR que faz isso e não estou encontrando (ainda) nada similar no Linux. O objetivo é poder compartilhar dados reais (data.frame) pela lista de discussão sem a necessidade de conexões, read.table, etc. O que preciso é algo assim como o resultado para o data.frame BOD que fica no pacote datasets: structure(list(Time = c(1, 2, 3, 4, 5), demand = c(8,3, 10,3, 19, 16, 15.6, 19.8), .Names = c("Time", "demand"), row.names = c(NA, -6L), class = "data.frame") Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

Média e desvio de muitas variáveis de um fatorial!
by Ivan Bezerra Allaman 18 Out '11

18 Out '11

Bom dia senhores! Tenho um experimento fatorial 14x4, com 12 variáveis, e gostaria de saber se tem como colocar de uma só vez em um único quadro a média e o desvio padrão dos tratamentos e de todas as variáveis. Fiz algo como: medias <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,mean,na.rm=TRUE)) desvios <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,sd,na.rm=TRUE)) descrit <- list() for(i in 1:12){ descrit[[i]] <- cbind(medias[,i],desvios[,i]) } Mais sai como lista não sendo o ideal que procuro. Desde já agradeço pela atenção! \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 3

Dois plots no mesmo frame
by Marcelo Claro de Souza 18 Out '11

18 Out '11

Gostaria de plotar duas regressões lineares no mesmo frame. Para tal plotei os pontos utilizando os comandos abaixo: x1=data$SLA1 y1=data$LLS 1 x2=data$SLA2 y2=data$LLS2 plot(x1, y1, type='p', xlim=range(x1,x2), ylim=range(y1, y2), xlab='x', ylab='y') points(x2, y2, type='p', col="red", xlab='x', ylab='y') reg1=lm(y1~x1) reg2=lm(y2~x2) Mas ao tentar aplicar as linhas de tendência utilizando abline(reg1) e abline(reg2) não obtive resultado algum. Como devo proceder para adicionar as linhas das regressões a essas figuras? Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Strings
by Rodrigo Sant'Ana 18 Out '11

18 Out '11

Caros Companheiros da R-BR, estou com um problema, que provavelmente, seja muito simples de resolver. Mas confesso, já quebrei a cabeça para resolver e num consegui chegar ao resultado que desejo. Abaixo, tento replicar meu problema para que vocês possam me ajudar! Muito obrigado pela atenção de todos! # Eis uma réplica simplificada da base de dados! dados <- data.frame(municipio = c("A", "B", "C", "D", "E"), pesc.emb = c("PESCADOR JOAO","EMBARCACAO ZE","PESCADOR FELIPE","EMBARCACAO AMARGURA II","EMBARCACAO STAR")) dados municipio pesc.emb 1 A PESCADOR JOAO 2 B EMBARCACAO ZE 3 C PESCADOR FELIPE 4 D EMBARCACAO AMARGURA II 5 E EMBACACAO STAR # Preciso criar uma coluna classificando a informação da linha em "PESCADOR" ou "EMBARCACAO", ou seja, preciso extrair a string "PESCADOR" ou "EMBARCACAO" da coluna pesc.emb. Já tentei utilizar todas as funções do pacote stringr, porém não obtive exito, pois preciso ambos na mesma coluna. library(stringr) dados$classify <- str_extract(dados$pesc.emb, "PESCADOR") dados$classify2 <- str_extract(dados$pesc.emb, "EMBARCACAO") dados municipio pesc.emb classify classify2 1 A PESCADOR JOAO PESCADOR <NA> 2 B EMBARCACAO ZE <NA> EMBARCACAO 3 C PESCADOR FELIPE PESCADOR <NA> 4 D EMBARCACAO AMARGURA II <NA> EMBARCACAO 5 E EMBARCACAO STAR <NA> EMBARCACAO sei que pode parecer simples, mas estou enfrentando dificuldades em solucionar este problema. Fico super agradecido pela atenção de todos, abs _______________ Rodrigo Sant'Ana Oceanógrafo

6 10

Diagrama de Vernn (Package eVenn)
by Íris Cardoso 17 Out '11

17 Out '11

Olá, pessoal Primeiramente, sou nova na lista e gostaria de dizer que é muito bom encontrar um grupo com uma troca tão grande e rica de informações. E, como toda mensagem tem um "interesse" por trás, gostaria de esclarecer uma dúvida: Apanhei do R um bom tempo para conseguir plotar um diagrama de Venn (package eVenn). Sempre acusava um erro do tipo ausência de "TRUE/FALSE" nos dados. Optei por alterar o único ponto diferente dos dados dos exemplos, os nomes das variáveis. Depois de alterar os nomes para ficarem iguais aos nomes da base do exemplo "Liste_1_.xxx." (onde xxx = soma de 1's na variável), consegui plotar o diagrama. Porém ele plotou a 1ª coluna (id das observações) como uma variável do diagrama. Excluí a referida coluna e obtive o Diagrama do jeito que queria. Gostaria de saber se alguém já conseguiu plotar usando o "nome verdadeiro" das variáveis ou se essa alteração é mesmo necessária. Achei estranho pois, se o nome da variável fosse mesmo necessário nesse formato, ele não deveria plotar a coluna de id. Desde já agradeço a atenção. P.S.: A coluna de id estava com o mesmo nome do exemplo "row.names" e o conteúdo estava configurado como character. Íris Cardoso Estatística e Metre em Logística

2 2

metamodelagem
by Edmar Alino 17 Out '11

17 Out '11

Srs, Estou tentando criar metamodelos para otimizacao e me recomendaram os seguintes pacotes: - Redes Neurais (rbf e mlp): AMORE - GaSP: Dicekriging. Como sou iniciante no R, qualquer sugestao/recomendacao/orientacao pertinente a esses pacotes ou outros pacotes mais apropriados serao muito bemvindos. Grato desde ja, Edmar

2 1

Ajuda
by Paloma Palmieri 17 Out '11

17 Out '11

Prezados, Boa tarde! Estou com o seguinte problema: Tenho um programa em MATLAB que gera um arquivo .prn que é utilizado por um arquivo compilado em C (que é chamado em um outro arquivo de extensão .bat). O que acontece é que eu tenho que rodar os programas individualmente e isso não é muito eficiente. Neste caso, gostaria de ver se alguém da lista sabe como eu chamo e rodo um executável (ou .bat) em C dentro do MATLAB. Desde já agradeço. Abraços, Paloma Palmieri Alves paloma.palmieri(a)gmail.com Cel.: (21) 8811-5230 (21) 7961-8718

2 1

ordenação de data para plot
by julio cesar araujo 15 Out '11

15 Out '11

Boa noite pessoal Tenho uma série temporal de retornos de ativos e gostaria de plotar um gráfico simples com data no eixo X e os valores no eixo Y.....como existem feriados no Brasil, essa série não é de segunda a sexta. Possuo +- 2000 datas, e quando vou plotar plot(data, valor) os dados do eixo X não ficam em ordem cronologica. O que fazer para que fiquem as datas no eixo X em ordem cronologica e os valores no eixo Y em linha? Atc

2 2

Re: [R-br] ordenação de data para plot
by Ivan Bezerra Allaman 15 Out '11

15 Out '11

É provável que seus dados não esteja formatado da maneira que o R entende, ou seja, YY/MM/DD, se eu não me engano. Caso queira formatar as datas da maneira que mais lhe convém, verifique as funções, ?as.Date Para uma ajuda mais refinada nos forneça um CRM. (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

Re: [R-br] Selecionar colunas em uma matriz
by Thiago De paula protásio 14 Out '11

14 Out '11

Muito obrigado Luis....Era exatamente 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

Selecionar colunas em uma matriz
by Thiago De paula protásio 14 Out '11

14 Out '11

Pessoal; Estou com uma dúvida e gostaria da ajuda de vcs....Preciso selecionar duas colunas específicas em um matriz para construir um gráfico tipo biplot().Quero selecionar as colunas 1 e 3. Alguém poderei me ajudar? Muito obrigado a todos. 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

2 1

Ajuda com mistura de distribuições
by Alexandre Castagna 14 Out '11

14 Out '11

Caros, Possuo dados distribuídos em um padrão que aparenta ser resultado da soma de duas gaussianas e gostaria de tentar reconstruir essas duas gaussianas. Existe algum pacote no R que auxilie neste trabalho? Obrigado, -- Alexandre Castagna Mourão e Lima Laboratório de Radioecologia e Mudanças Globais - UERJ Rua São Francisco Xavier, 524, Maracanã Pav. Haroldo Lisboa da Cunha, Sub-solo

2 2

Desdobramento de contrastes GLM via glht()
by ASANTOS 14 Out '11

14 Out '11

Boa dia pessoal, Gostaria de saber se existe alguma forma de se atribuírem índices ao teste realizado com a função glht() do pacote multcomp, exemplificando: require(multcomp) y1<-c(rpois(100, lambda = 5),rpois(100, lambda = 10),rpois(100, lambda = 5),rpois(100, lambda = 7))##Geração da variável resposta trat<-c(rep(1,100),rep(2,100),rep(3,100),rep(4,100))##Geração dos tratamentos repeticao<-rep(1:100,4)##Imaginado um DIC colocando as repetições trat<-as.factor(trat) m.1<-glm(y1~trat+repeticao,family="poisson")##Ajustando GLM Poisson anova.m.1<-anova(m.1,test="Chi") anova.m.1##Verificando a significancia dos tratamentos e se não há efeito das repetições summary(glht(m.1,linfct=mcp(trat="Tukey")))##Desdobrando os tratamentos Como a função utiliza o Tukey acho bem complicado como no caso simulado, decidir entre tratamentos que possuem índice ambíguo, por exemplo, alguém conhece alguma forma de calcular a decisão da distribuição dos índices, Obrigado, Alexandre **Obs: Descupe por re-posts mais acredito que haja algum problema entre o yahoo mail e o c3sl, pois ou minhas msg não chegam ou chegam com um dia de atraso.

2 2

Pacote AER não instala
by ASANTOS 14 Out '11

14 Out '11

Boa noite pessoal, Estou tentando instalar o pacote AER, mas não consigo me dá a mensagem de erro: Attaching package: 'zoo' The following object(s) are masked from 'package:base': as.Date ** help *** installing help indices ** building package indices ... ** testing if installed package can be loaded Error : object 'nobs' not found whilst loading namespace 'AER' ERROR: loading failed * removing '/home/asantos/R/i686-pc-linux-gnu-library/2.12/AER' Warning in install.packages : installation of package 'AER' had non-zero exit status The downloaded packages are in '/tmp/Rtmpo2u5y0/downloaded_packages' Não sei se pode ser isso mais tenho instalado dois pacotes com o nome zoo, Obrigado, Alexandre

4 10

Pacote AER não instala
by ASANTOS 14 Out '11

14 Out '11

Estou tentando fazer pelo terminal mais ele me diz que já estou com a última versão, mas na verdade continuo com a 2.12, sendo: asantos@asantos-Aspire-5315:~$ sudo apt-get install r-base Lendo listas de pacotes... Pronto Construindo árvore de dependências Lendo informação de estado... Pronto r-base já é a versão mais nova. Os seguintes pacotes foram instalados automaticamente e já não são necessários: qgis-plugin-grass-common linux-headers-2.6.38-10-generic linux-headers-2.6.38-10 Use 'apt-get autoremove' para removê-los. 0 pacotes atualizados, 0 pacotes novos instalados, 0 a serem removidos e 1 não atualizados. On 10/14/2011 08:52 AM, Benilton Carvalho wrote: > O RStudio e' apenas uma janelinha com frico-frico que usa o R que vc > tinha anteriormente instalado no computador. Ou seja, o R e' > independente do RStudio e vc pode instala-lo usando o metodo padrao > (vai no website, baixa o binario, instala). > _______________________________________________ > 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. >

1 0

[Dúvida]Escala de um mapa
by Pedro Rafael 14 Out '11

14 Out '11

Como verificar a escala de um mapa no R? Estou com um arquivo shp e não sei quantos km equivale 1 cm no mapa. Existe algum comando para saber isso rapidamente? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Abrir URL usando o R
by Lucas Ferreira 14 Out '11

14 Out '11

Boa tarde.Gostaria de saber como entrar em um URL específico e extrair um determinado texto dentro do corpo da página.consegui abrir a pagina usando o pacote scrape, mas não consegui fazer uma procura por caracteres dentro da pagina. alguém pode ajudar? ------------------------------------------------------------------ -Lucas Ferreira da Silva Mestrando em Biotecnologia - UFESLaboratório de Biotecnologia Aplicada ao AgronegócioNúcleo de Biotecnologia/CCS/UFESAv. Marechal Campos, 146829040-090, Vitória, ES, BrasilTel/Fax. 55.27. 3335-7348 ----------------------------------------------------------------- ------------------------------------------------------------------

5 5

Circunferência centrada no centroide de raio especificado
by Pedro Rafael 14 Out '11

14 Out '11

Pessoal como desenhar uma circunferência centrada no centroide de um município de uma malha importada no R com raio específico? Gostaria de de especificar a quilometragem e o R desenhar uma circunferência com o raio com a quilometragem que eu atribui. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Cbind e Loop
by Lrg Finanças 13 Out '11

13 Out '11

Tenho vários objetos (de 1 coluna e várias linhas) com nomes sequenciais: A.1, A.2, A.3 ..., B.1, B.2,B.3, ..., C.1, ... Pretendo juntar os arquivos em colunas. Ex.: todos os arquivos com nome A. Tentei. for (i in 1:7){ cbind(get(paste("A",i,sep=".")))} Mas não funciona. Como posso juntar esses objetos? Grato, Luciano.

3 3

[OFF] Morre Denis Ritchie
by Pedro Rafael 13 Out '11

13 Out '11

Pessoal vi agora que o grande Denis Ritchie morreu e até agora não vi uma reportagem na TV falando do falecimento desse grande cientista. Para quem não o conhece ele é o responsável para criação da linguagem C e do sistema Unix. Quem nunca estudou uma linguagem de programação e criou um programa “hello, world”? Programadores também se lembram de Ritchie por seu famoso programa “hello, world”, que é muito utilizado em códigos de programação como um exemplo de um programa de computador bastante simples. http://g1.globo.com/tecnologia/noticia/2011/10/morre-dennis-ritchie-pai-do-… -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Desdobramento de contrastes GLM via glht()
by ASANTOS 13 Out '11

13 Out '11

Boa noite pessoal, Gostaria de saber se existe alguma forma de se atribuírem índices ao teste realizado com a função glht() do pacote multcomp, exemplificado: require(multcomp) y1<-c(rpois(100, lambda = 5),rpois(100, lambda = 10),rpois(100, lambda = 5),rpois(100, lambda = 7))##Geração da variável resposta trat<-c(rep(1,100),rep(2,100),rep(3,100),rep(4,100))##Geração dos tratamentos repeticao<-rep(1:100,4)##Imaginado um DIC colocando as repetições trat<-as.factor(trat) m.1<-glm(y1~trat+repeticao,family="poisson")##Ajustando GLM Poisson anova.m.1<-anova(m.1,test="Chi") anova.m.1##Verificando a significancia dos tratamentos e se não há efeito das repetições summary(glht(m.1,linfct=mcp(trat="Tukey")))##Desdobrando os tratamentos Como a função utiliza o Tukey acho bem complicado como no caso simulado, decidir entre tratamentos que possuem índice duplo, por exemplo, alguém conhece alguma forma de calcular a decisão de distribuição dos índices, Obrigado, Alexandre

1 0

comando diagnosis
by Alexandro (Yahoo) 13 Out '11

13 Out '11

#Boa tarde a todos, numa regressão logística múltipla, após a escolha de um ponto de corte, > >#criei a tabela de classificação conforme CMR abaixo, e usei o >#comando dignosis do pacote DiagnosisMed >#para calcular a sensibilidade e especificidade ! > > >require(DiagnosisMed) > > >tab<-matrix(c(64,112,17,406),nrow=2,ncol=2,byrow=T, >dimnames=list(Classificado=c(" 1 "," 0 "), >Observado=c(" 1 "," 0 ")));tab > > >diagnosis(tab)[] > > > >#Porém os valores estão trocados, o certo é >#Sensitividade =64 /(17+64) =79,01% >#Especificidade=406/(112+406)=78,378% >#Tentei usar o comando: > > > > >diagnosis(t(tab))[] > > > > >#Mas não deu certo. Sendo assim: >#1- o que ocorreu ? >#2- Tem alguma função nos pacotes epiR ou ROCR que >#calcula a sensibilidade e especificidade ? >#3- para achar o ponto de corte "ideal" tenho que fazer várias tabelas de classificação >#e plotar a sens. x espec. >#No "satanic", tem-se um comando do proc logistic que >#mostra todas as sensibilidade e especificidade para diferentes pontos de cortes. >#Tem algum comando similar no R para regressão logística múltipla ? >#no caso univariado tem-se: >y=c(0,1,0,0,1,1,1,0,0,0) >x=c(12,19,14,13,14,15,16,17,17,18) >require(DiagnosisMed);arearoc<-DiagnosisMed::ROC(y,x,Full=TRUE) > > > > >########### MUITO OBRIGADO ########### >#Alex

2 1

Erro em simulação
by Greice Laureano 13 Out '11

13 Out '11

Bom dia, Estou realizando um estudo de simulação e esta ocorrendo o seguinte erro: "Erro: não é possível alocar vetor de tamanho 2.7 Mb" Existe alguma forma de solucionar este problema de memória? Obrigado! Att, Greice.

3 4

glm.fit: fitted probabilities numerically 0 or 1 occurred
by Marcelo Luiz de Laia 13 Out '11

13 Out '11

Estou tentando analisar um grupo de dados de contagem da seguinte maneira: gm1 <- glmer(cbind(Germinacao, 20 - Germinacao) ~ Dia + Especie + Matriz + (Dia|Repeticao),family = binomial, data = ka.dat) Mas, esta mostrando uma advertencia: Warning message: glm.fit: fitted probabilities numerically 0 or 1 occurred Eu vi que nas paginas 197-9 do MASS4 existe uma explicação para esse fenômeno (que a principio somente ocorre quando se usa a binomial). Confesso que não entendi se a advertencia é grave ou se eu posso relevá-la. Eu tenho 3 repeticoes, somente. Seria isso? Se desejar ler as paginas acima, me peça que digitalizo e envio. Obrigado -- Marcelo Brazil Linux user number 487797

3 4

Como criar um vetor/matriz de valores NA?
by Thiago Touya 12 Out '11

12 Out '11

Senhores, Estou criando um banco de dados para fazer um estudo onde se levam em conta as safras. Na verdade, por não conseguir uma maneira mais elegante, usando o for, criei 12 bases de dados, cada uma representando um mês. Vou unir essas bases numa base maior que será o objeto de minha análise. Meu problema reside no fato de que no segundo mês terei que colocar um vetor-linha de valores NA para cada uma das variáveis; no terceiro mês uma matriz com duas linhas, no quarto mês uma matriz com 3 linha e assim até uma matriz com 11 linhas e quatro colunas com todas as entradas sendo NA's. Não estou conseguindo criar esses vetor/matrizes de valores faltantes. Tentei usar c(,,,,) ou c( , , , ) ou c(.,.,.,.) para o vetor e a partir dele estender o resultado para as matrizes, mas não deu certo. Alguém tem alguma sugestão? Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

2 3

[glmmADMB] Correta modelagem para dados de contagem
by Marcelo Luiz de Laia 12 Out '11

12 Out '11

Estou tentando analisar um conjunto de dados de contagem e fui cair na pagina do pacote glmmADMB. No entanto, achei pouca informação que pudesse me ajudar. Meus dados estão assim: Esp Mat Time Rep Germ Inf Sucupira M02 1 1 0 0 Sucupira M02 2 1 1 0 Sucupira M02 3 1 1 8 Sucupira M02 4 1 1 8 Sucupira M02 5 1 1 9 Sucupira M02 6 1 1 10 Sucupira M02 7 1 1 15 Sucupira M02 1 2 0 0 (...) Sucupira M02 7 3 1 10 Sucupira M04 1 1 0 0 Sucupira M04 2 1 0 1 (...) Vinhatico M15CV15 7 3 5 8 Esp = Espécie Mat = Matriz (três matrizes de locais diferentes) Time = Dias em que foi analisado Rep = Repetição (gerbox com 20 sementes cada) Germ = Quantidade de sementes que germinaram Inf = Quantidade de sementes com micro-organismos Trata-se de medidas repetidas no tempo, ou seja, no dia 1 contaram-se quantas sementes germinaram e quantas apresentavam micro-organismos. No dia 2 contaram-se quantas sementes germinaram e quantas apresentavam micro-organismos. Assim até o dia 7. Tentei o seguinte código: library(glmmADMB) fit = glmm.admb(Germ~Esp*Mat+Rep,random=~Rep,group="Time", zeroInflation=TRUE,data=ka.dat,family="poisson") Mas, além do erro: Maximum number of iterations exceeded in dvector eigen(_CONST dmatrix& m) Error in glmm.admb(Germ ~ Esp * Mat + Rep, random = ~Rep, : The function maximizer failed In addition: Warning message: running command './nbmm -maxfn 500 ' had status 1 não tenho certeza se a modelagem proposta por mim esta correta. Por exemplo, o 'group' deveria ser mesmo o tempo? E o efeito aleatorio seria mesmo so a repeticao? Eu tambem gostaria de ver se ha alguma relacao entre germinacao e infestacao. Se voce tiver experiencia com esse pacote, poderia dar algumas dicas? Obrigado -- Marcelo Brazil Linux user number 487797

2 1

[Dúvida] Revolution R Enterprise ou R?
by Pedro Rafael 12 Out '11

12 Out '11

Pessoal existe alguma vantagem de se utilizar o Revolution R Enterprise que é livre para uso acadêmico em vez de utilizar o R? Vejo que no site da empresa eles elogiam bastante o R mais que o Revolution é muito superior, principalmente no fato que no Revolution eu posso analisar base de dados gigantescas e por dar suporte à programação paralela. Agora uma inversão da pergunta feita anteriormente: Existe alguma vantagem de se utilizar o Revolution R Enterprise no lugar do R puramente? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 1