- R-br - DINF UFPR

ajuste com nls
by Gilson Sanchez 11 Jul '11

11 Jul '11

Pessoal mais uma vez procuro ajuda de vc's, consegui ajustar um modelo com nls, mas agora o que quero saber como eu faço para estimar o coeficiente (R2) e ver se a equação ajustou corretamente. obrigado -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

3 2

[Notícias-Dúvidas] Novidades de interfaces gráficas no R.
by Pedro Rafael 11 Jul '11

11 Jul '11

Pesquisando mais sobre algumas novidades do R vi que existem vários projetos de construção de interfaces gráficas para o R fora o RStudio, R Comander, Tinn-R e RKward que sãos os mais conhecidos. Existe o Red-R que tenta facilitar a vida do usuário através de uma programação visual, uma espécie de SAS GUIDE que é uma funcionalidade do SAS bastante elogiada pelos usuários SAS. Observe o vídeo do Red-R em: http://www.youtube.com/watch?v=A4WwhOnuC10&feature=player_embedded. Existe um outro projeto também livre chamado R Analytic Flow, veja: http://www.ef-prime.com/products/ranalyticflow_en/. Amigos vocês conhecem mais alguns projetos como esse? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico (Secretaria de Saúde de Estado da Paraíba).

2 1

Re: [R-br] Identidade de modelos!!!
by Thiago De paula protásio 10 Jul '11

10 Jul '11

Walmes; No meu caso aqui ajustei modelos não lineares para 11 tratamentos e gostaria de saber se posso utilizar um único modelo para todos eles.... Muito obrigado pela informação.... Abraço! 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

Identidade de modelos!!!
by Thiago De paula protásio 10 Jul '11

10 Jul '11

Prezados membros da R-br; Venho novamente pedir a ajuda de vcs...Gostaria de saber como posso proceder ao teste de identidade de modelos no R?Ajustei algumas equações aqui para 6 tratamentos e visualmente percebi certa similaridade entre elas....dessa forma, gostaria de proceder a algum teste para testar a hipótese de que apenas uma equação seria necessária para todos os tratamentos.... Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Extrair células de uma matrix
by Diego Pujoni 10 Jul '11

10 Jul '11

Olá Pessoal, Tenho uma matriz onde cada coluna contém valores de probabilidade variando de 0 a 1. Como posso excluir os valores 0's e 1's de cada coluna, ficando apenas com as probabilidades intermediárias? Toda coluna possui apenas um valor 0 e apenas um valor 1, mas eles estão em linhas diferentes em cada coluna. Fiz um exemplo abaixo > A = matrix (c(1,.5,.8,.9,0.,.8,.5,0,1,.4,.8,0,.7,.3,1),5,3) > A [,1] [,2] [,3] [1,] 1.0 0.8 0.8 [2,] 0.5 0.5 0.0 [3,] 0.8 0.0 0.7 [4,] 0.9 1.0 0.3 [5,] 0.0 0.4 1.0 Quero obter essa matriz > A [,1] [,2] [,3] [1,] 0.5 0.8 0.8 [2,] 0.8 0.5 0.7 [3,] 0.9 0.4 0.3 Obrigado Diego PJ

2 1

regressão linear pelo origem
by Gilson Sanchez 09 Jul '11

09 Jul '11

Pessoal mais uma vez procuroa juda de vc's, bom a verdade é pra um colega ele é novo com o R, se alguem puder ajudar. Estou precisando da ajuda de você, sou novo no grupo e estou tentando fazer um grafico de regressao linear simples ajustando o intercepto da reta deve passar pela origem dos eixos (0,0) mas não estou conseguindo. O CMR que estou trabalhando é o seguinte. # REGRESSÃO EI30 VS DS require(XLConnect) reg<-loadWorkbook("Estadistica Final.xls", create= FALSE) reg1<-readWorksheet(reg, sheet="Plan1") reg1 attach(reg1) cor.test(EI30, DS, alternative = c("two.sided"), method = c("pearson"), conf.level = 0.95) reglin1<-lm(DS~EI30, data=reg1) anova(reglin1) reglin1 summary(reglin1) # GRAFICOS par(mfrow = c(3,1), mar=c(4.5,5,2,2)) plot(EI30,DS, type="p", pch=15, lwd=1.5, ylim=c(0,100), xlim=c(0,50), cex.main=1.2, cex.lab=0.9, cex.axis=0.9, xlab="", ylab="",cex.lab=2, main="Descoberto", font.main=2) abline(reglin1, lwd=2.0, col="black", lty=1) ## pode-se usar lty=9 legend("bottomright", bty="n", lwd=c(1.5), pch=c(15), col=c("black"), legend=c(expression(y=='1,516x'~+~'9,241'~','~R^2=='0,58'^'*')), cex=1) plot(EI30, LS, pch=16, col="black", ylim=c(0,3), xlim=c(0,50), cex.main=1.2, cex.lab=0.9, cex.axis=0.9, xlab="", ylab=expression(bold(Perda~do~solo~(g~L^{-1}))), font.lab=2, main="Liteira", font.main=2) abline(reglin2, lwd=2.0, col="black", lty=2) legend("bottomright", bty="n", lwd=c(1.5), pch=c(16), col=c("black"), legend=c(expression(y=='0,018x'~+~'0,323'~','~R^2=='0,16'^'*')), cex=1) Aguardo suas consideraçoes. Att. Omar Cubas Encinas Instituto Nacional de Pesquisas da Amazônia / INPA Aluno de Mestrado em Agricultura do Tropico Úmido / ATU Manaus - AM - Brasil (092) 9102-2542 -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

3 3

Multivariada - Clustering
by Emerson Cotta Bodevan 08 Jul '11

08 Jul '11

Prezados, qual pacote do R é mais adequado para análise multivariada? Especificamente, clustering (gráficos de dendograma,etc). Abs, -- *Emerson*

4 4

Kappa
by Edson Lira 08 Jul '11

08 Jul '11

Boa tarde R-Br lista!!! Para as tabelas abaixo, como faço para calcular o índice de concordância de Kappa e sua signficância. TABELA 1 OURO negative positive OB1 negative 449 3 positive 6 7 TABAELA2 OURO negative positive OB2 positive 1 8 TABELA 3 OURO negative positive OB3 negative 0 1 positive 1 7 Estou usando o pacote Psych, mas não vi nenhuma medida de signficância do Kappa. Alguem tem alguma sugestão? Edson Lira Estatístico Manaus-Amazonas

4 3

Re: [R-br] funcao tableStack do epicalc
by Pedro Emmanuel Alvarenga Americano do Brasil 07 Jul '11

07 Jul '11

Paula, Algumas dicas pra poder rodar essa funcao. Em primeiro lugar ela foi feita pra rodar tudo dE uma vez,nao faz sentido chamar pra grupos de variageis diferentes. Nessa funcao as coisas dao pau por conta do formato de uma ou mais variaveis. Nao use character, somente factor. Confira se as factor possui algum nivel que nao tem dado. Faca a funcao progressivamente inserindo as variaveis uma a uma. De uma olhada nos exemplos pra ver te da uma - paula - luz. :-) Pedro Brasil via android (:)= Em 07/07/2011 17:49, "Paula Mendes Luz" <luzpaulamendes(a)gmail.com>escreveu: oi benilton com certeza. por acaso não estou com o código em mãos agora e por isso tentei reproduzir em palavras... mas entendo q é mais complicado assim. então depois envio outro email. p. 2011/7/7 Benilton Carvalho <beniltoncarvalho(a)gmail.com>: > Paula, > > se vc puder seguir as recomendacoes do guia de postagem (vide rodape') > e fornecer um ...

1 0

funcao tableStack do epicalc
by Paula Mendes Luz 07 Jul '11

07 Jul '11

oi pessoal estou tentando usar a funcao tableStack da library epicalc mas tá complicado... os erros que recebo nao sao mto informativos... alguem tem alguma dica? tem algum tipo de variavel que a funcao nao aceita? acho que o problema tá nisso (apesar de q os erros sao diversos)... lendo o help acho que entendi que a funcao só aceita um tipo de var de cada vez... entao se tenho vars numericas e character, devo separar as vars... por ex, tenho idade (numeric), gender (character: 1/2), education (character: 1/2/3/4/9), altura (numeric), centro (character: 1/2) faco o tableStack de idade e altura usando o by = centro depois o tableStack de gender e education, by = centro, com vars.to.factor de gender e education (preciso transformar o centro tb?) ex de erros recebidos incluem: object 'table0' not found e object 'test.method' not found qq ajuda é bem vinda! valeu! paula.

2 2