Julho 2012 - R-br - DINF UFPR

Combinar linhas entre conjuntos...
by andrebvs 17 Jul '12

17 Jul '12

3 6

Análise Fatorial com dummy
by Felipe Buchbinder 17 Jul '12

17 Jul '12

Professores, Eu tenho uma base de dados com variáveis dummies e gostaria de fazer uma redução de dimensionalidade. Inicialmente pensei em uma análise fatorial com correlação tetracórica. Porém, a matriz de correlação tinha correlações unitárias fora da diagonal principal sem razão aparente (olhando o gráfico, havia pontos nos quatro quadrantes, mas mesmo assim a correlação dava 1!) Como minhas variáveis não são discretizações de variáveis contínuas, não sei se o uso de correlação tetracórica é adequado. Pensei em usar o *phi* de Pearson como medida de correlação, mas ele dá idêntico à correlação de Pearson. Outra idéia foi fazer a extração usando análie por componentes principais para dados categóricos (CATPCA) mas, pelo que andei pesquisando, o resultado não dá muito diferente de uma análise fatorial convencional (com correlação de Pearson e extração por PCA) Eu queria poder contar com uma opinião sobre o que fazer: Usar correlação tetracórica (neste caso, como eu resolvo o fato da matriz de correlação não ser positiva semidefinida?) Usar o phi de Pearson como medida de correlaçao Usar a correlação de Pearson, pura se simplesmente Extrair os fatores por CATPCA ou outra alternativa, que é sempre bem-vinda. Desde já agradeço a atenção e a colaboração, Felipe Buchbinder

2 2

Re: [R-br] Ajuda programação!
by Fernando Mayer 17 Jul '12

17 Jul '12

(encaminhando para a lista para manter arquivado) O tab ali é o objeto criado pelo Ivan na sua mensagem inicial. A única modificação que fiz foi no nome das linhas: k1 <- c(1.8850,2.1575,2.4775,2.7975) k2 <- c(2.9375,2.2475,2.7300,2.9950) tab <- matrix(c(k1,k2),ncol=2) rownames(tab) <- c('p1n1','p1n2','p2n1','p2n2') []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/7/16 Mauro Sznelwar <sznelwar(a)uol.com.br>: > Estava olhando o vosso script, e não consegui reproduzir onde tem o tab, > alega que não reconhece este objeto. Precisa de alguma biblioteca pra isto? > > > Também não sei se entendi a pergunta direito, mas se vc quer apenas um > "atalho" para chamar determinada linha usando apenas os valores de "a" > e "b", pode fazer algo como > > f <- function(a, b){ > x <- paste("p", a, "n", b, sep = "") > return(x) > } > >> f(a = 1, b = 2) > [1] "p1n2" >> f(a = 1, b = 1) > [1] "p1n1" >> tab[f(1,1), ] > k1 k2 > 1.8850 2.9375 >> as.vector(tab[f(1,1),]) > [1] 1.8850 2.9375 > > > []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 > >

1 0

Duvidas em análise cluster
by edlayne moraes 17 Jul '12

17 Jul '12

Boa tarde, Inseri uma planilha com os dados em .csv, coloquei na linha 1º o nome dos meus dados. Ao importar a tabela eu mandei ler essa linha como rótulo das variáveis. Para fazer o tratamento eu instalei o pacote "vegan" e ao mandar plotar o gráfico o programa informa um erro de "rowsums (x, na.rm=TRUE): x deve ser numérico". Alguém sabe me dizer o que há de erro e se existe outra maneira de fazer o cluster no programa? Desde já agradeço a atenção Att, -- Edlayne Pimentel de Moraes**** Graduanda em Geologia-UFES**** (28) 81117-0585//(28) 3552-0194

4 3

Ajuda programação!
by Ivan Bezerra Allaman 17 Jul '12

17 Jul '12

Bom dia senhores! Suponhamos que tem-se dois níveis de um fator "p" e dois níveis de um fator "n", no qual os níveis do fator "p" irei representar pela letra "a" e os níveis do fator "n" pela letra "b". Então tem-se a seguinte configuração: k1 <- c(1.8850,2.1575,2.4775,2.7975) k2 <- c(2.9375,2.2475,2.7300,2.9950) tab <- matrix(c(k1,k2),ncol=2) rownames(tab) <- c('panb','panb,'panb','panb') colnames(tab) <- c('k1','k2') tab k1 k2 panb 1,8850 2,9375 panb 2,1575 2,2475 panb 2,4775 2,7300 panb 2,7975 2,9950 Vamos supor que eu queira os valores do nível 1 do fator "p" e nível 1 do fator "n". Então, a=1 e b=1; p1n1 1,8850 2,9375 Agora nível 1 do fator "p" e nível 2 do fator "n". Então, a=1, b=2; p1n2 2,1575 2,2475 Tem como automatizar isto em função destes dois parâmetros? algo como: medias = as.vector(tab[ab, ])???? A coluna não interessa, podem ficar tranquilo! Apenas as linhas interessam de acordo com os parâmetros. Se alguém puder me dar uma ajuda ficarei grato! Abraços! S,f.P (Allaman) \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 3

Colocar elipse de confiança num gráfico
by Mauro Sznelwar 17 Jul '12

17 Jul '12

Como eu faço para inserir uma elipse de confiança neste gráfico acima de Análise discriminante, não sei como usar o comando: plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) library(ellipse) plot(ellipse(S, centre = , t = ), type = "l") Abaixo tem o script que gerou este gráfico, e estou enviando também o arquivo de dados. http://www.datafilehost.com/download-1eccadf2.html marajoara=read.csv('Marajoara Mario.csv',head=T,dec=',',sep=';') marajoara.norm=sapply(data.frame(marajoara[,2:15]),scale) #marajoara.norm=(data.frame(marajoara[,2:15])) #summary(marajoara.norm) #library(fBasics) #basicStats(marajoara.norm) marajoara.norm.h=hclust(dist(marajoara.norm),method='average') #marajoara.norm.h=hclust(dist(marajoara.norm)^2,method='ward') plot(marajoara.norm.h,main='Marajoara original-Ward Linkage\n Square Euclidian Distance\n Hclust Dendograma') #pdf('figuras/Marajoara-original-dendograma-hclust-average.pdf') library(cluster) marajoara.norm.d=diana(dist(marajoara.norm)) plot(marajoara.norm.d,main='Marajoara original-Divisive Dendograma') #pdf('figuras/Marajoara-original-dendograma-divisive.pdf') #Componentes Principais initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km1=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) library(MASS) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster), col=as.numeric(marajoara.norm.km1$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma hclust 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendograma hclust 3') initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.d,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km2=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km2$cluster), col=as.numeric(marajoara.norm.km2$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma diana 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendodiana-3.pdf') marajoara.norm.km=kmeans(marajoara.norm,3) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km$cluster), col=marajoara.norm.km$cluster) title('PCA Marajoara original - clusters kmeans 3 centros') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-clusters-kmeans-3-centros.pdf') #Análise Discriminante Linear marajoara.norm.ldakm1=lda(marajoara.norm,as.character(marajoara.norm.km1$cluster)) plot(marajoara.norm.ldakm1,dimen=2,cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-Grupo do K-Means - corte dendograma hclust em 3') #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendohclust-3.pdf') #fazendo um zoom no gráfico corta um ponto plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-parte do gráfico para visualização\n Grupos do K-Means-corte dendograma hclust em 3') #identify(marajoara.norm.km1$cluster,n=5) #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendodiana-3.pdf') library(ellipse) plot(ellipse(S, centre = , t = ), type = "l")

1 0

Re: [R-br] Digest R-br, volume 17, assunto 14
by Gilenio Borges Fernandes 17 Jul '12

17 Jul '12

Muito obrigado ao Elias T. Krainski. A sugestão foi fundamental para a simulação que pretendo fazer Gilenio Fernandes Em 16 de julho de 2012 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 programação! (Ivan Bezerra Allaman) > 2. duvida controle de qualidade (Natalia Martins) > 3. Re: duvida controle de qualidade (Eder David Borges da Silva) > 4. Re: Ajuda programação! (Lucas Barbosa) > 5. Re: duvida controle de qualidade (Natalia Martins) > 6. Re: duvida controle de qualidade > (henrique jose de paula alves paula alves) > 7. Re: duvida controle de qualidade (Elton Gean Araújo Araújo) > 8. Re: Ajuda programação! (Walmes Zeviani) > 9. Re: Ajuda para Simulação de variáveis binárias > correlacionadas no R (Elias T. Krainski) > 10. Duvidas em análise cluster (edlayne moraes) > 11. Re: Duvidas em análise cluster (Natalia Martins) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 15 Jul 2012 09:07:18 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Ajuda programação! > Message-ID: > <1342368438.19845.YahooMailNeo(a)web161801.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia senhores! > > Suponhamos que tem-se dois níveis de um fator "p" e dois níveis de um > fator "n", no qual os níveis do fator "p" irei representar pela letra "a" e > os níveis do fator "n" pela letra "b". Então tem-se a seguinte configuração: > > k1 <- c(1.8850,2.1575,2.4775,2.7975) > k2 <- c(2.9375,2.2475,2.7300,2.9950) > tab <- matrix(c(k1,k2),ncol=2) > rownames(tab) <- c('panb','panb,'panb','panb') > colnames(tab) <- c('k1','k2') > tab > k1 k2 > panb 1,8850 2,9375 > panb 2,1575 2,2475 > panb 2,4775 2,7300 > panb 2,7975 2,9950 > > > Vamos supor que eu queira os valores do nível 1 do fator "p" e nível 1 do > fator "n". Então, a=1 e b=1; > > p1n1 1,8850 2,9375 > > Agora nível 1 do fator "p" e nível 2 do fator "n". Então, a=1, b=2; > > p1n2 2,1575 2,2475 > > Tem como automatizar isto em função destes dois parâmetros? > > algo como: > > medias = as.vector(tab[ab, ])???? > > A coluna não interessa, podem ficar tranquilo! Apenas as linhas interessam > de acordo com os parâmetros. > > Se alguém puder me dar uma ajuda ficarei grato! > > Abraços! > > S,f.P > (Allaman) > > > > > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120715/e65d4cd4/atta… > > > > ------------------------------ > > Message: 2 > Date: Sun, 15 Jul 2012 17:36:31 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: Lista R <R-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] duvida controle de qualidade > Message-ID: > < > CAGPGL2FkzWBmCnkL96JZpG9-VQ06EyU-x-3OyanyM34xR+1DBA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Por favor estou com um dúvida em controle de qualidade e ficaria muito > grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > Bem é o seguinte, tenho uma máquina que é programada para colocar uma > quantidade exata de adubo em um ponto especificado no campo via GPS, ou > seja, a maquina vai passando pelo campo e derramando o adubo na quantidade > exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a > cada ponto. > Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as > quantidades que a maquina *realmente derruba.* > Quando a maquina termina o seu serviço um equipamento instalado nela me > informa o que ela realmente derrubou de adubo em cada coordenada, assim eu > tenho que verificar se esta maquina está sob controle estatístico de > qualidade. > - tenho o que a maquina deve derramar por coordenada no campo todo > - tenho o que a maquina realmente derramou por coordenada no campo todo > > Minhas duvidas em relação a esse processo são: > > Para verificar se meu processo esta sob controle eu fiz a diferença entre o > que a maquina realmente aplicou e o que deveria ser aplicado, eu posso > aplicar a teoria de controle de qualidade a esta diferença? > Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também pode > ser aplicada a todos os dados? > Quando faço o teste de normalidade de shapiro wilks este mostra me que meus > dados não apresentam normalidade, devo transformá-los? > > Muito obrigada, > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120715/2a5a2916/atta… > > > > ------------------------------ > > Message: 3 > Date: Sun, 15 Jul 2012 19:50:19 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > <CALKkXXqsbB7ACDufTkRL= > mHngB06gtx-196ptJHKHVxpBhsTYw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Natalia, > Estes dados são legais....eu seguiria uma outro caminho de analise, como a > distribuição é desenvolvida no campo por faixas (A maquina vai passando) e > muito possivelmente em regiões distintas da lavoura por questão de > topografia o erro deve ser diferente, uma analise espacial de seus dados > seria muito legal, pois supostamente o erro não > é distribuído aleatoriamente dentro do talhão, de uma olhada neles por este > ponto de vista podem aparecer informações legais. > Att > > Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com > >escreveu: > > > Por favor estou com um dúvida em controle de qualidade e ficaria muito > > grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > > Bem é o seguinte, tenho uma máquina que é programada para colocar uma > > quantidade exata de adubo em um ponto especificado no campo via GPS, ou > > seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > > exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a > > cada ponto. > > Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as > > quantidades que a maquina *realmente derruba.* > > Quando a maquina termina o seu serviço um equipamento instalado nela me > > informa o que ela realmente derrubou de adubo em cada coordenada, assim > eu > > tenho que verificar se esta maquina está sob controle estatístico de > > qualidade. > > - tenho o que a maquina deve derramar por coordenada no campo todo > > - tenho o que a maquina realmente derramou por coordenada no campo todo > > > > Minhas duvidas em relação a esse processo são: > > > > Para verificar se meu processo esta sob controle eu fiz a diferença entre > > o que a maquina realmente aplicou e o que deveria ser aplicado, eu posso > > aplicar a teoria de controle de qualidade a esta diferença? > > Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > > pode ser aplicada a todos os dados? > > Quando faço o teste de normalidade de shapiro wilks este mostra me que > > meus dados não apresentam normalidade, devo transformá-los? > > > > Muito obrigada, > > > > -- > > > > Natália da Silva Martins > > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > > Contato: (19) 8306-4743 > > > > > > > > _______________________________________________ > > 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/20120715/7975b619/atta… > > > > ------------------------------ > > Message: 4 > Date: Sun, 15 Jul 2012 21:36:24 -0300 > From: Lucas Barbosa <lucasbcr(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Ajuda programação! > Message-ID: > <CAPhSkDFxTn-e8=ZtoKmctr6MqOqYK-ff0cak1= > MVdyENYF9tjg(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Vê aí se é isso mesmo o que vc quer: > > > expand.grid(1:2, 1:2) > Var1 Var2 > 1 1 1 > 2 2 1 > 3 1 2 > 4 2 2 > > apply(expand.grid(1:2, 1:2), 1, function(x){paste('p', x[1], 'n', x[2], > sep='')} ) > [1] "p1n1" "p2n1" "p1n2" "p2n2" > > > Em 15 de julho de 2012 13:07, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> escreveu: > > Bom dia senhores! > > > > Suponhamos que tem-se dois níveis de um fator "p" e dois níveis de um > fator > > "n", no qual os níveis do fator "p" irei representar pela letra "a" e os > > níveis do fator "n" pela letra "b". Então tem-se a seguinte configuração: > > > > k1 <- c(1.8850,2.1575,2.4775,2.7975) > > k2 <- c(2.9375,2.2475,2.7300,2.9950) > > tab <- matrix(c(k1,k2),ncol=2) > > rownames(tab) <- c('panb','panb,'panb','panb') > > colnames(tab) <- c('k1','k2') > > tab > > k1 k2 > > panb 1,8850 2,9375 > > panb 2,1575 2,2475 > > panb 2,4775 2,7300 > > panb 2,7975 2,9950 > > > > Vamos supor que eu queira os valores do nível 1 do fator "p" e nível 1 do > > fator "n". Então, a=1 e b=1; > > > > p1n1 1,8850 2,9375 > > > > Agora nível 1 do fator "p" e nível 2 do fator "n". Então, a=1, b=2; > > > > p1n2 2,1575 2,2475 > > > > Tem como automatizar isto em função destes dois parâmetros? > > > > algo como: > > > > medias = as.vector(tab[ab, ])???? > > > > A coluna não interessa, podem ficar tranquilo! Apenas as linhas > interessam > > de acordo com os parâmetros. > > > > Se alguém puder me dar uma ajuda ficarei grato! > > > > Abraços! > > > > S,f.P > > (Allaman) > > > > > > > > > > \begin{signature} > > <<>>= > > Prof. Dr. Ivan Bezerra Allaman > > Universidade Estadual de Santa Cruz > > Departamento de Ciências Exatas e Tecnológicas > > Ilhéus/BA - Brasil > > Fone: +55 73 3680-5596 > > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > > @ > > \end{signature} > > > > _______________________________________________ > > 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: Mon, 16 Jul 2012 08:33:58 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > <CAGPGL2FJvm5SxhU8p4RUJWEAu5v3MphJ-crxoBH8= > y-QXiTFmw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Oi Eder, bom dia! > muito obrigada pela sugestão, mas uma análise espacial já foi realizada > sim, agora o objetivo era verificar realmente a eficiência da maquina, e > como análises preliminares já fiz um teste de comparação de médias o qual > resultou que a média de aplicação real é diferente da média da aplicação > que deveria ser realizada. > > Att, > > > > > Em 15 de julho de 2012 19:50, Eder David Borges da Silva > <eder(a)leg.ufpr.br>escreveu: > > > Natalia, > > Estes dados são legais....eu seguiria uma outro caminho de analise, como > a > > distribuição é desenvolvida no campo por faixas (A maquina vai passando) > e > > muito possivelmente em regiões distintas da lavoura por questão de > > topografia o erro deve ser diferente, uma analise espacial de seus dados > > seria muito legal, pois supostamente o erro não > > é distribuído aleatoriamente dentro do talhão, de uma olhada neles por > este > > ponto de vista podem aparecer informações legais. > > Att > > > > Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com > >escreveu: > > > >> Por favor estou com um dúvida em controle de qualidade e ficaria muito > >> grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > >> Bem é o seguinte, tenho uma máquina que é programada para colocar uma > >> quantidade exata de adubo em um ponto especificado no campo via GPS, ou > >> seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > >> exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a > >> cada ponto. > >> Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as > >> quantidades que a maquina *realmente derruba.* > >> Quando a maquina termina o seu serviço um equipamento instalado nela me > >> informa o que ela realmente derrubou de adubo em cada coordenada, assim > eu > >> tenho que verificar se esta maquina está sob controle estatístico de > >> qualidade. > >> - tenho o que a maquina deve derramar por coordenada no campo todo > >> - tenho o que a maquina realmente derramou por coordenada no campo todo > >> > >> Minhas duvidas em relação a esse processo são: > >> > >> Para verificar se meu processo esta sob controle eu fiz a diferença > entre > >> o que a maquina realmente aplicou e o que deveria ser aplicado, eu posso > >> aplicar a teoria de controle de qualidade a esta diferença? > >> Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > >> pode ser aplicada a todos os dados? > >> Quando faço o teste de normalidade de shapiro wilks este mostra me que > >> meus dados não apresentam normalidade, devo transformá-los? > >> > >> Muito obrigada, > >> > >> -- > >> > >> Natália da Silva Martins > >> Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > >> Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > >> Contato: (19) 8306-4743 > >> > >> > >> > >> _______________________________________________ > >> 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. > > > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/d4006b68/atta… > > > > ------------------------------ > > Message: 6 > Date: Mon, 16 Jul 2012 09:12:52 -0300 > From: henrique jose de paula alves paula alves <jpahenrique(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > < > CAEfjKJTXP8R3BEH_C_WkxuLDieL+qjDZX9_6WxwgBNPTTzuxOA(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > escreva pra este professor aqui da ufop, ele pode lhe ajudar melhor. o > cara é fera nesta área. > > > fernandoest(a)ig.com.br > > > Em 16 de julho de 2012 08:33, Natalia Martins <nsmbarreto(a)gmail.com> > escreveu: > > Oi Eder, bom dia! > > muito obrigada pela sugestão, mas uma análise espacial já foi realizada > sim, > > agora o objetivo era verificar realmente a eficiência da maquina, e como > > análises preliminares já fiz um teste de comparação de médias o qual > > resultou que a média de aplicação real é diferente da média da aplicação > que > > deveria ser realizada. > > > > Att, > > > > > > > > > > Em 15 de julho de 2012 19:50, Eder David Borges da Silva < > eder(a)leg.ufpr.br> > > escreveu: > > > >> Natalia, > >> Estes dados são legais....eu seguiria uma outro caminho de analise, > como a > >> distribuição é desenvolvida no campo por faixas (A maquina vai > passando) e > >> muito possivelmente em regiões distintas da lavoura por questão de > >> topografia o erro deve ser diferente, uma analise espacial de seus dados > >> seria muito legal, pois supostamente o erro não é distribuído > aleatoriamente > >> dentro do talhão, de uma olhada neles por este ponto de vista podem > aparecer > >> informações legais. > >> Att > >> > >> Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com> > >> escreveu: > >>> > >>> Por favor estou com um dúvida em controle de qualidade e ficaria muito > >>> grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > >>> Bem é o seguinte, tenho uma máquina que é programada para colocar uma > >>> quantidade exata de adubo em um ponto especificado no campo via GPS, ou > >>> seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > >>> exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas > a > >>> cada ponto. > >>> Eu sei a quantidade que a maquina deve derrubar em todo o campo e as > >>> quantidades que a maquina realmente derruba. > >>> Quando a maquina termina o seu serviço um equipamento instalado nela me > >>> informa o que ela realmente derrubou de adubo em cada coordenada, > assim eu > >>> tenho que verificar se esta maquina está sob controle estatístico de > >>> qualidade. > >>> - tenho o que a maquina deve derramar por coordenada no campo todo > >>> - tenho o que a maquina realmente derramou por coordenada no campo todo > >>> > >>> Minhas duvidas em relação a esse processo são: > >>> > >>> Para verificar se meu processo esta sob controle eu fiz a diferença > entre > >>> o que a maquina realmente aplicou e o que deveria ser aplicado, eu > posso > >>> aplicar a teoria de controle de qualidade a esta diferença? > >>> Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > >>> pode ser aplicada a todos os dados? > >>> Quando faço o teste de normalidade de shapiro wilks este mostra me que > >>> meus dados não apresentam normalidade, devo transformá-los? > >>> > >>> Muito obrigada, > >>> > >>> -- > >>> > >>> Natália da Silva Martins > >>> Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > >>> Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > >>> Contato: (19) 8306-4743 > >>> > >>> > >>> > >>> _______________________________________________ > >>> 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. > > > > > > > > > > -- > > > > Natália da Silva Martins > > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > > Contato: (19) 8306-4743 > > > > > > > > _______________________________________________ > > 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. > > > > -- > Henrique José de Paula Alves > Graduando em Estatística > UFOP - MG > > > ------------------------------ > > Message: 7 > Date: Mon, 16 Jul 2012 05:20:05 -0700 (PDT) > From: Elton Gean Araújo Araújo <egarauj(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > <1342441205.79483.YahooMailNeo(a)web126004.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia Natália! > Dá uma olhada no pacote qcc do R, esse pacote realiza analises de controle > de qualidade e pode gerar alguns gráficos para verificar se a máquina > adubadeira está bem regulada ou não. > Elton. > > > ________________________________ > De: henrique jose de paula alves paula alves <jpahenrique(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Segunda-feira, 16 de Julho de 2012 9:12 > Assunto: Re: [R-br] duvida controle de qualidade > > escreva pra este professor aqui da ufop, ele pode lhe ajudar melhor. o > cara é fera nesta área. > > > fernandoest(a)ig.com.br > > > Em 16 de julho de 2012 08:33, Natalia Martins <nsmbarreto(a)gmail.com> > escreveu: > > Oi Eder, bom dia! > > muito obrigada pela sugestão, mas uma análise espacial já foi realizada > sim, > > agora o objetivo era verificar realmente a eficiência da maquina, e como > > análises preliminares já fiz um teste de comparação de médias o qual > > resultou que a média de aplicação real é diferente da média da aplicação > que > > deveria ser realizada. > > > > Att, > > > > > > > > > > Em 15 de julho de 2012 19:50, Eder David Borges da Silva < > eder(a)leg.ufpr.br> > > escreveu: > > > >> Natalia, > >> Estes dados são legais....eu seguiria uma outro caminho de analise, > como a > >> distribuição é desenvolvida no campo por faixas (A maquina vai > passando) e > >> muito possivelmente em regiões distintas da lavoura por questão de > >> topografia o erro deve ser diferente, uma analise espacial de seus dados > >> seria muito legal, pois supostamente o erro não é distribuído > aleatoriamente > >> dentro do talhão, de uma olhada neles por este ponto de vista podem > aparecer > >> informações legais. > >> Att > >> > >> Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com> > >> escreveu: > >>> > >>> Por favor estou com um dúvida em controle de qualidade e ficaria muito > >>> grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > >>> Bem é o seguinte, tenho uma máquina que é programada para colocar uma > >>> quantidade exata de adubo em um ponto especificado no campo via GPS, ou > >>> seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > >>> exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas > a > >>> cada ponto. > >>> Eu sei a quantidade que a maquina deve derrubar em todo o campo e as > >>> quantidades que a maquina realmente derruba. > >>> Quando a maquina termina o seu serviço um equipamento instalado nela me > >>> informa o que ela realmente derrubou de adubo em cada coordenada, > assim eu > >>> tenho que verificar se esta maquina está sob controle estatístico de > >>> qualidade. > >>> - tenho o que a maquina deve derramar por coordenada no campo todo > >>> - tenho o que a maquina realmente derramou por coordenada no campo todo > >>> > >>> Minhas duvidas em relação a esse processo são: > >>> > >>> Para verificar se meu processo esta sob controle eu fiz a diferença > entre > >>> o que a maquina realmente aplicou e o que deveria ser aplicado, eu > posso > >>> aplicar a teoria de controle de qualidade a esta diferença? > >>> Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > >>> pode ser aplicada a todos os dados? > >>> Quando faço o teste de normalidade de shapiro wilks este mostra me que > >>> meus dados não apresentam normalidade, devo transformá-los? > >>> > >>> Muito obrigada, > >>> > >>> -- > >>> > >>> Natália da Silva Martins > >>> Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > >>> Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > >>> Contato: (19) 8306-4743 > >>> > >>> > >>> > >>> _______________________________________________ > >>> 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. > > > > > > > > > > -- > > > > Natália da Silva Martins > > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > > Contato: (19) 8306-4743 > > > > > > > > _______________________________________________ > > 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. > > > > -- > Henrique José de Paula Alves > Graduando em Estatística > UFOP - MG > _______________________________________________ > 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/20120716/2a471d52/atta… > > > > ------------------------------ > > Message: 8 > Date: Mon, 16 Jul 2012 10:48:16 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Ajuda programação! > Message-ID: > <CAFU=EkY3OHAz8KJCLv3QZ+0ZSRd2= > Ak3KGJNSHKFVDoi7beWHQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Coloca as colunas que identificam os fatores junto com as respostas num > data.frame() e usa tapply(), aggregate(), by(), etc para obter resumos por > grupos. Será que entendi a pergunta? > > À 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/20120716/b10f32cd/atta… > > > > ------------------------------ > > Message: 9 > Date: Mon, 16 Jul 2012 07:26:43 -0700 (PDT) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Ajuda para Simulação de variáveis binárias > correlacionadas no R > Message-ID: > <1342448803.39401.YahooMailNeo(a)web163403.mail.gq1.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Olá Gilenio, seria algo assim: > > ### simula dados > > n <- 1000; beta <- c(-3,1) > y <- integer(n) > y[1] <- rbinom(1,1,.5) > for (i in 2:n) > y[i] <- rbinom(1, 1, 1/(1+exp(-beta[1]-beta[2]*y[i-1]))) > > ### estima modelo > > m <- glm(y ~ y2, family=binomial, data=data.frame(y=y[2:n], y2=y[2:n-1])) > summary(m) > > > Elias T. Krainski > > > >________________________________ > > De: Gilenio Borges Fernandes <gilenio(a)ufba.br> > >Para: r-br(a)listas.c3sl.ufpr.br > >Enviadas: Sábado, 14 de Julho de 2012 20:52 > >Assunto: [R-br] Ajuda para Simulação de variáveis binárias > correlacionadas no R > > > > > >Simulação de variáveis binárias > correlacionadas no R é bem popular, usando os pacotes bindata, bild, > binarySimCLF, > BinNor e binGroup. Mas eu estou precisando simular variáveis binárias > correlacionadas num modelo auto regressivo de primeira ordem de regressão > logística, > a partir dos parâmetros do preditor linear, incluindo o referente à > variável dependente > defasada (Yi-1). Alguém poderia me indicar algum pacote ou função no > R para este fim? > >Desde já, agradeço-- > >Gilenio Borges Fernandes > >Universidade Federal da Bahia > >Instituto de Matemática > >Departamento de Estatística > >Av. Adhemar de Barros, s/n ? Ondina. > >40.170-110 - Salvador - BA, Brasil > >Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 > >URL: http://lattes.cnpq.br/6764860618464860 > > > > > >_______________________________________________ > >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/20120716/0532d17a/atta… > > > > ------------------------------ > > Message: 10 > Date: Mon, 16 Jul 2012 11:45:41 -0300 > From: edlayne moraes <moraes.epgeo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Duvidas em análise cluster > Message-ID: > <CAPbDHAgFkwcrt= > BzyHsgeQRhwYgBWOaWBG-aUzjGgqYbey1-Tg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde, > Inseri uma planilha com os dados em .csv, coloquei na linha 1º o nome dos > meus dados. Ao importar a tabela eu mandei ler essa linha como rótulo das > variáveis. > Para fazer o tratamento eu instalei o pacote "vegan" > e ao mandar plotar o gráfico o programa informa um erro de "rowsums (x, > na.rm=TRUE): x deve ser numérico". > Alguém sabe me dizer o que há de erro e se existe outra maneira de fazer o > cluster no programa? > Desde já agradeço a atenção > Att, > > -- > > Edlayne Pimentel de Moraes**** > > Graduanda em Geologia-UFES**** > (28) 81117-0585//(28) 3552-0194 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/1476e8f7/atta… > > > > ------------------------------ > > Message: 11 > Date: Mon, 16 Jul 2012 11:56:22 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Duvidas em análise cluster > Message-ID: > <CAGPGL2EzM27XQ-EpcfKRvT+YzuxX+PD=Ena+Cob= > 8jHSJMJApg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Oie, tenta na hora de ler os dados colocar o comando head=T, assim ele lerá > a 1 linha como o nome das suas variaveis, ai na 1 coluna só haverá valores > numericos. > > Att, > Em 16 de julho de 2012 11:45, edlayne moraes <moraes.epgeo(a)gmail.com > >escreveu: > > > Boa tarde, > > Inseri uma planilha com os dados em .csv, coloquei na linha 1º o nome dos > > meus dados. Ao importar a tabela eu mandei ler essa linha como rótulo das > > variáveis. > > Para fazer o tratamento eu instalei o pacote "vegan" > > e ao mandar plotar o gráfico o programa informa um erro de "rowsums (x, > > na.rm=TRUE): x deve ser numérico". > > Alguém sabe me dizer o que há de erro e se existe outra maneira de fazer > o > > cluster no programa? > > Desde já agradeço a atenção > > Att, > > > > -- > > > > Edlayne Pimentel de Moraes**** > > > > Graduanda em Geologia-UFES**** > > (28) 81117-0585//(28) 3552-0194 > > > > > > _______________________________________________ > > 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. > > > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/429ebf6b/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 17, assunto 14 > ***************************************** > -- Gilenio Borges Fernandes Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

1 0

Ajuda para Simulação de variáveis binárias correlacionadas no R
by Gilenio Borges Fernandes 16 Jul '12

16 Jul '12

Simulação de variáveis binárias correlacionadas no R é bem popular, usando os pacotes bindata, bild, binarySimCLF, BinNor e binGroup. Mas eu estou precisando simular variáveis binárias correlacionadas num modelo auto regressivo de primeira ordem de regressão logística, a partir dos parâmetros do preditor linear, incluindo o referente à variável dependente defasada (Yi-1). Alguém poderia me indicar algum pacote ou função no R para este fim? Desde já, agradeço -- Gilenio Borges Fernandes Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

3 2

Colocar elipse de confiança num gráfico
by Mauro Sznelwar 14 Jul '12

14 Jul '12

Como eu faço para inserir uma elipse de confiança neste gráfico acima de Análise discriminante, não sei como usar o comando: plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) library(ellipse) plot(ellipse(S, centre = , t = ), type = "l") Abaixo tem o script que gerou este gráfico, e estou enviando também o arquivo de dados. http://www.datafilehost.com/download-1eccadf2.html marajoara=read.csv('Marajoara Mario.csv',head=T,dec=',',sep=';') marajoara.norm=sapply(data.frame(marajoara[,2:15]),scale) #marajoara.norm=(data.frame(marajoara[,2:15])) #summary(marajoara.norm) #library(fBasics) #basicStats(marajoara.norm) marajoara.norm.h=hclust(dist(marajoara.norm),method='average') #marajoara.norm.h=hclust(dist(marajoara.norm)^2,method='ward') plot(marajoara.norm.h,main='Marajoara original-Ward Linkage\n Square Euclidian Distance\n Hclust Dendograma') #pdf('figuras/Marajoara-original-dendograma-hclust-average.pdf') library(cluster) marajoara.norm.d=diana(dist(marajoara.norm)) plot(marajoara.norm.d,main='Marajoara original-Divisive Dendograma') #pdf('figuras/Marajoara-original-dendograma-divisive.pdf') #Componentes Principais initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km1=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) library(MASS) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster), col=as.numeric(marajoara.norm.km1$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma hclust 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendograma hclust 3') initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.d,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km2=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km2$cluster), col=as.numeric(marajoara.norm.km2$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma diana 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendodiana-3.pdf') marajoara.norm.km=kmeans(marajoara.norm,3) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km$cluster), col=marajoara.norm.km$cluster) title('PCA Marajoara original - clusters kmeans 3 centros') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-clusters-kmeans-3-centros.pdf') #Análise Discriminante Linear marajoara.norm.ldakm1=lda(marajoara.norm,as.character(marajoara.norm.km1$cluster)) plot(marajoara.norm.ldakm1,dimen=2,cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-Grupo do K-Means - corte dendograma hclust em 3') #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendohclust-3.pdf') #fazendo um zoom no gráfico corta um ponto plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-parte do gráfico para visualização\n Grupos do K-Means-corte dendograma hclust em 3') #identify(marajoara.norm.km1$cluster,n=5) #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendodiana-3.pdf') library(ellipse) plot(ellipse(S, centre = , t = ), type = "l")

1 0

Dúvida sobre o ANOSIM
by Débora Rodrigues de Souza 13 Jul '12

13 Jul '12

Boa tarde a todos, Acabei de encontrar essa página de discussões, e achei muito interessante a ideia, parabéns aos organizadores. Bom, gostaria de auxílio para executar uma análise de similaridade do tipo ANOSIM. Sei que o teste é feito a posteriori o NMDS, ou alguma análise de ordenação, mas tenho dúvida quando ao script do mesmo. Já tentei seguindo os passos indicados no manul do programa R, mas fiquei em dúvida em qual comando devo alterar, para então inserir os meus dados. Estou mexendo nesse programa a pouco tempo, por isso a dificuldade. Ao que me parece, devo entrar com dois tipos de planilhas. Mas mesmo colando meus dados, o resultado que aparece é o da memória do programa. Como alterar isso??? Esse é o script que vejo no programa, data(dune) data(dune.env) dune.dist <- vegdist(dune) attach(dune.env) dune.ano <- anosim(dune.dist, Management) summary(dune.ano) plot(dune.ano) Alguém poderia me ajudar??? Desde já agradeço a atenção, e aguardo resposta. Att, Débora

2 5