- R-br - DINF UFPR

Re: [R-br] Ajuda: duvida passando uma string como modelo para uma função
by Jose Claudio Faria 17 Abr '11

17 Abr '11

Walmes, Embora ache que essa função ainda possa ser melhorada, sua ajuda foi valiosa! Sua sugestão pode ser mais simples, e próximo do que realmente estou precisando, veja: mF <- function(x=seq(1, 5, by=1), y=seq(0, 10, by=5), model){ z <- outer(x, y, function(x, y){ eval(model) }) z } mF(model=quote(2 - 0.4*y - 0.5*x + 0.8*sqrt(x))) > mF(model=quote(2 - 0.4*y - 0.5*x + 0.8*sqrt(x))) [,1] [,2] [,3] [1,] 2.300000 0.3000000 -1.700000 [2,] 2.131371 0.1313708 -1.868629 [3,] 1.885641 -0.1143594 -2.114359 [4,] 1.600000 -0.4000000 -2.400000 [5,] 1.288854 -0.7111456 -2.711146 Talvez na segunda-feira outras pessoas possam dar alguma contribuição... Abs e muito obrigado! -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 3

Ajuda: duvida passando uma string como modelo para uma função
by Jose Claudio Faria 17 Abr '11

17 Abr '11

Pessoal, Gostaria de passar uma string como modelo para uma função, e não estou conseguindo. (Acho que já estou ficando velho...) De uma forma primitiva consigo: mF <- function(x=seq(1, 5, by=1), y=seq(0, 10, by=5), z=function(x, y) 2 - 0.4*y - 0.5*x + 0.8*sqrt(x)) { z <- outer(x, y, z) z } > mF() [,1] [,2] [,3] [1,] 2.300000 0.3000000 -1.700000 [2,] 2.131371 0.1313708 -1.868629 [3,] 1.885641 -0.1143594 -2.114359 [4,] 1.600000 -0.4000000 -2.400000 [5,] 1.288854 -0.7111456 -2.711146 Mas o que preciso (para automatizar um tanto de coisas é algo assim): mF <- function(x=seq(1, 5, by=1), y=seq(0, 10, by=5), z='2 - 0.4*y - 0.5*x + 0.8*sqrt(x)') { z <- outer(x, y, function(x, y) DUVIDA) z } com o mesmo resultado anterior. Alguma dica? Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

Spline Cúbico Interpolador
by Tiago de Oliveira Vieira 17 Abr '11

17 Abr '11

Bom dia a todos! Eu estou com uma problema para desenvolver um algoritmo, e implementá-lo no R, para criar passo-a-passo uma função de spline cúbico interpolador. Alguém, por favor, tem uma sugestão ou indicação bibliográfica. Desde de já, obrigado! Att, Tiago Vieira.

4 5

[Re: Spline Cúbico Interpolador]]
by Felipe Emanoel Barletta Mendes 17 Abr '11

17 Abr '11

Tiago, Segue link http://leg.ufpr.br/lib/exe/fetch.php/pessoais:referencial.pdf Com um pdf que escrevi sobre Modelos GAM e SPLINES. No livro do Simon há exemplos de implementações splines no R. Espero que ajude!! Att Paulo Justiniano > Tiago > > neste caso,onde vcquere programar e nao utilizar o que jah esta pronto > me parece que voce deveria comecar pelo > metodo, consultando alguma bibliografica como as indicadas pelo Felipe > ouaida Green e Silvermann. > > O spline pode ser escrito como um modelo linear segmentado. > Comece por ai: como escrever uma funcao segmentada (com nos fixos) > e dai por diante. > Comece com algo simples, em 2/3 partes, depois veja como generalizar. > > bom trabalho! > > > > > On Fri, 15 Apr 2011, Tiago de Oliveira Vieira wrote: > >> Eu preciso entender o algoritmo, já abri a função e não entendi. >> >> Em Sex, 2011-04-15 às 13:08 +0100, Benilton Carvalho escreveu: >>> essa funcao nao ja' esta' implementada? >>> >>> splines::spline >>> >>> b >>> >>> 2011/4/15 Tiago de Oliveira Vieira <tiagovieira7(a)yahoo.com.br>: >>>> Bom dia a todos! Eu estou com uma problema para desenvolver um >>>> algoritmo, e implementá-lo no R, para criar passo-a-passo uma função >>>> de >>>> spline cúbico interpolador. Alguém, por favor, tem uma sugestão ou >>>> indicação bibliográfica. >>>> >>>> Desde de já, obrigado! >>>> >>>> Att, >>>> >>>> Tiago Vieira. >>>> >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)listas.c3sl.ufpr.br >>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>> >>> >>> >>> >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>_______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Felipe E. Barletta Mendes (41)9189-5198 (41)3025-2150 (41)3328-7216 http://www.leg.ufpr.br/doku.php/pessoais:felipe Felipe E. Barletta Mendes (41)9189-5198 (41)3025-2150 (41)3328-7216 http://www.leg.ufpr.br/doku.php/pessoais:felipe Felipe E. Barletta Mendes (41)9189-5198 (41)3025-2150 (41)3328-7216 http://www.leg.ufpr.br/doku.php/pessoais:felipe

2 1

Res: Experimento Fatorial! Desdobramento da interação!
by Thiago De paula protásio 16 Abr '11

16 Abr '11

Walmes; Muito obrigado pela dica..... A função fat2.dic() funcionou perfeitamente...... E foi bem mais simples.... Bom final de semana! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

Re: [R-br] graficos de regressão linear
by Walmes Zeviani 16 Abr '11

16 Abr '11

GIlson, Segue algumas possibilidades. # opção 1 plot(1, xlab=expression(bold(PPDF)), sub=expression(bold((mu~mol~.~fótons~.~m^{-2}~.~s^{-1})))) # opção 2 plot(1, xlab="") mtext(text=expression(bold(PPDF)), side=1, line=2) mtext(text=expression(bold((mu~mol~.~fótons~.~m^{-2}~.~s^{-1}))), side=1, line=3) plot(1, xlab="") legend("topleft", bty="n", legend=expression(y=='49,81'~+~'52,24'^(-x/'1,31')~.~R^2=='0,71')) # não pode usar virgula, pois é separador de argumentos e não decimal # números com virgulas escrever como string # não usar ~ pre/pos sinal de == # ver demo(plotmath) # para mudar globalmente as saidas de texto/gráfica . para , use options(OutDec=",") Procure enviar e-mails para lista. Sua chance de ser respondido imediatamente aumentam consideravelmente, as contribuições dos demais inscritos podem ser bem mais interessantes, e toda troca de informações fica a disposição de todos para consulta futura. À 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 ========================================================================== Em 16 de abril de 2011 10:58, <gilsonsch(a)gmail.com> escreveu: > > Mais uma vez eu com minas perguntas, > > no gráfico que eu estou fazendo e que vc ajudou muito, como eu posso > quebrar o titulo do eixo > > por exemplo eu tenho o script (abixo) e quero quebrar aonde aparace \n > > xlab=expression(bold(PPDF(\n)~(mu~mol~.~fótons~.~m^{-2}~.~s^{-1}))) > > > > Outro, como eu posso colocar a equação ajustada da curva na legenda não > quero como texto (y=49,81+52,24^(-x/1,31), R2=0,71) > > legend("topright", pch=c(19,5,1,15), lty=c(1,2,3,4), > legend=c(expression (Antes~da~antese~.~-8*degree~C), > > Legend=c(expression(y=~49,81~+~52,24^(-x/1,31)~.~R^2~=~0,71) (eu quero > colocar esta equação na legenda) > > <quote author='Walmes Zeviani 2'> > Gilson, > > Fiz um CMR (código mínimo reproduzível) que está disponível ná página de > dicas curtas de R (ridiculas, > > http://www.leg.ufpr.br/doku.php/ridiculas?&#representacao_de_ajuste_de_curv… > ). > Confira. > > À disposição. > Walmes. > > Em 23 de março de 2011 16:02, gilsonsch <gilsonsch(a)gmail.com> escreveu: > > > Caro Wlames, grato pela ajuda > > mas ainda to com algumas duvidas, ainda pra crear gráficos, bom plotar 4 > > curvas ajustadas num mesmo gráfico > > > > tenho os seguintes dados > > > > antese semanas tem polen antes 1 -8 1 antes 2 -8 2 antes 3 -8 3 antes > 4 > > -8 4 antes 5 -8 5 antes 6 -8 6 antes 7 -8 7 antes 8 -8 8 antes 9 -8 9 > > antes 10 -8 10 antes 11 -8 11 antes 12 -8 12 antes 1 8 13 antes 2 8 14 > > antes 3 8 15 antes 4 8 16 antes 5 8 17 antes 6 8 18 antes 7 8 19 antes 8 > 8 > > 20 antes 9 8 21 antes 10 8 22 antes 11 8 23 antes 12 8 24 depois 1 -8 25 > > depois 2 -8 26 depois 3 -8 27 depois 4 -8 28 depois 5 -8 29 depois 6 -8 > 30 > > depois 7 -8 31 depois 8 -8 32 depois 9 -8 33 depois 10 -8 34 depois 11 -8 > > 35 depois 12 -8 36 depois 1 8 37 depois 2 8 38 depois 3 8 39 depois 4 8 > 40 > > depois 5 8 41 depois 6 8 42 depois 7 8 43 depois 8 8 44 depois 9 8 45 > > depois 10 8 46 depois 11 8 47 depois 12 8 48 > > > > Agradeçp ajuda > > > > > > 2011/3/23 Walmes Zeviani 2 [via R-br] <[hidden > > email]<http://user/SendEmail.jtp?type=node&node=3400321&i=0&by-user=t> > > > > > > >> Gilson, > >> > >> Pode-se fazer o gráfico dos dados observados (plot()) e usar lines() > para > >> colocar as curvas ajustadas. Também pode usar a função matplot() e > >> matlines(). Confira a documentação. > >> > >> À disposição. > >> > >> Em 23 de março de 2011 14:05, gilsonsch <[hidden > >> email]<http://user/SendEmail.jtp?type=node&node=3400100&i=0&by-user=t> > >> > escreveu: > >> > >>> Pessoal, > >>> to precissando fazer um grafico de analises de regressão de tres > >>> espécies > >>> dende, caiaue, e hi, eu já conseguir fazer os grafico, mas eu quero que > >>> as > >>> três linhas estejam num grafico só. o mesmo quando eu quero grafico de > >>> dispersão, quero colocar mas de uma variave no mesmo gráfico, como eu > >>> faço > >>> isso. > >>> > >>> exemplo: > >>> dende 2 > >>> dende 2 > >>> dende 3 > >>> caiuae 4 > >>> caiaue 3 > >>> caiaue 3 > >>> hi 2.5 > >>> hi 3.4 > >>> hi 5.4 > >>> > >>> Gilson > >>> > >>> -- > >>> View this message in context: > >>> > http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… > < > http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… > > > >>> Sent from the R-br mailing list archive at Nabble.com. > >>> _______________________________________________ > >>> R-br mailing list > >>> [hidden > >>> email]<http://user/SendEmail.jtp?type=node&node=3400100&i=1&by-user=t> > >>> > >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >>> > >> > >> > >> > >> -- > >> > >> > ========================================================================== > >> 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: [hidden email] > >> <http://user/SendEmail.jtp?type=node&node=3400100&i=2&by-user=t> > >> twitter: @walmeszeviani > >> > >> > >> > >> homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 > >> > >> > >> > ========================================================================== > >> > >> > >> > >> _______________________________________________ > >> R-br mailing list > >> [hidden > >> email]<http://user/SendEmail.jtp?type=node&node=3400100&i=3&by-user=t> > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >> > >> > >> ------------------------------ > >> If you reply to this email, your message will be added to the > discussion > >> below: > >> > >> > http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… > < > http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… > > > >> > > > > > > > > -- > > 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. > > > > > > ------------------------------ > > View this message in context: Re: [R-br] graficos de regressão > > linear< > http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… > > > > > > Sent from the R-br mailing list > > archive<http://r-br.2285057.n4.nabble.com/>at Nabble.com. > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > > > > > -- > > ========================================================================== > 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 > ========================================================================== > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > </quote> > Quoted from: > > http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… >

1 0

Experimento Fatorial! Desdobramento da interação!
by Thiago De paula protásio 16 Abr '11

16 Abr '11

Prezado membros da R-Br; Recorro novamente a preciosa ajuda de vcs!! Estou analisando um experimento em esquema fatorial.... Tenho 3 clones em dois clones, com 4 repetições (fatorial completo). A interação clone*local é signficativa, mas não estou conseguindo desdobrar a interação... Estou usando os seguintes comandos....... anava1=aov(P~Local*Clone) summary(anava1) Df Sum Sq Mean Sq F value Pr(>F) Local 1 4.7082 4.7082 18.738 0.0004045 *** Clone 2 9.5002 4.7501 18.905 3.776e-05 *** Local:Clone 2 16.0320 8.0160 31.903 1.209e-06 *** Residuals 18 4.5227 0.2513 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 d2.anava1=aov(P~Clone/Local) effects(d2.anava1)[1:6] (Intercept) CloneI144 CloneI220 Clone3334:LocalPonte CloneI144:LocalPonte CloneI220:LocalPonte -69.489983 -1.850408 -2.465000 -1.502602 4.154252 1.106622 summary(d2.anava1, split = list("Clone:Local" = list(3334 = 1, I144 = 2, I220 = 3))) Erro: unexpected '=' in "summary(d2.anava1, split = list("Clone:Local" = list(3334 =" Mas não estou obtendo sucesso. Caso queiram posso enviar os dados! Agradeço-lhes a ajuda!! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

RES.: Gráficos de Superfície! Remover as linhas!
by Thiago De paula protásio 15 Abr '11

15 Abr '11

Muito obrigado Eder e demais membros da R-Br pela grande e valiosa ajuda!! Bom final de semana a todos!!! Att, 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

precisão numérica no R - uma falha?
by Cézar Freitas 15 Abr '11

15 Abr '11

Olá. Bem, escrevo apenas para ilustrar algo que me surpreendeu. Estava plotando gráficos e os colorindo segundo seu desvio padrão estava ou não acima de um ponto de corte. Reproduzo abaixo a parte do código que interessa: y=as.data.frame(matrix(c(4.96,5.02,4.97,5.03,4.96,5.03), ncol=3)) y$dp=apply(y,1,sd) > y V1 V2 V3 dp 1 4.96 4.97 4.96 0.005773503 2 5.02 5.03 5.03 0.005773503 refer=y$dp[1] y$dp<=refer [1] TRUE FALSE Sei que o R é extremamente preciso e até podemos alterar sua mantissa nos cálculos, mas isso não é rotineiro para cálculos tão simples. O pessoal do cálculo numérico saberia dizer o que evitar para não cair nessas armadilhas no futuro? Abraços, Cézar Freitas

6 10

RES.: Gráficos de Superfície! Remover as linhas!
by Thiago De paula protásio 15 Abr '11

15 Abr '11

Walmes e demais membros da R-Br; Utilizei os comandos baixos e obtive sucesso na elaboração de um gráfico de superfície. library(lattice) b0 <- 3696.046; b1 <- 62.183; b2 <- -33.414 g <- expand.grid(TCF = 68:80, RG = 32:37) g$PCS <- b0 + b1*g$TCF + b2*g$RG rgb.palette <- colorRampPalette(c("yellow", "green4", "blue"), space ="rgb") wireframe(PCS~TCF+RG, data=g, main="PCS=f(TCF+RG)", sub="Modelo Múltiplo", screen=list(z=30,x=-45), drape=T,border=F, col.regions=rgb.palette(90)) No entanto, eu gostaria de retirar as linhas do plano.Alguém conhece alguma função que me permite realizar tal tarefa??? Agradeço-lhes a ajuda!! 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

RES.:Gráficos de superfície!!
by Thiago De paula protásio 15 Abr '11

15 Abr '11

Walmes; Muito obrigado pela dica.. É isso mesmo que eu estava procurando.. 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

1 0

(sem assunto)
by Juraci Mendes Moreira 15 Abr '11

15 Abr '11

Estou tentando determinar um tamanho otimo de parcela com base em um experimento em branco, usando o modelo plato de resposta linear, mas deu esse erro e não sei como resolver. Usando a mesma rotina em analises anteriores deu certo. mas essa não. sem mais para o momento antecipo agradecimentos. jura Erro em qr.qty(QR, resid) : NA/NaN/Inf em chamada de função externa (argumento 5) Além disso: Warning message: In Ops.factor(lhs, rhs) : - not meaningful for factors > lrp.fit

2 1

Teste de t^2 multivariado de Hotellings
by Alexandre dos Santos 15 Abr '11

15 Abr '11

Boa tarde pessoal, Estava tentando comparar dois grupos de indivíduos com 5 variáveis, usando o teste de t^2 de Hotelling, ate ai tudo bem, mais quando simulei os indivíduos e suas variáveis com pelos menos um NA no vetor, a função não conseguia realizar o cálculo, dando a mensagem Erro em solve.default(v) : sistema é computacionalmente singular: condição recíproca número = 0 Tentei inserir na função na.omit(), sem sucesso, alguem poderia me dar um help, obrigado, segue a rotina: #####################Teste de t multivariado de Hotellings############################################ ##Depois de algum teste de cluster (ex. Mahalanobis) vou observar se os dois grupos possuem diferenças ##Grupo 1 ind.pop.1<-1:30 #######Os indivíduos med.pop.1<-data.frame(mma=Mod(rnorm(30)),mmc=Mod(rnorm(30)),mmd=Mod(rnorm(30 )),mca=Mod(rnorm(30)),mca2=Mod(rnorm(30)))#######As medidas dos indivíduos dados.pop.1<-cbind(ind.pop.1,med.pop.1)#######Planilha com os indivíduos e suas medidas head(dados.pop.1) ##Grupinho Oligo 3 ind.pop.3<-1:6 #######Os indivíduos med.pop.3<-data.frame(mma=Mod(rnorm(6)),mmc=Mod(rnorm(6)),mmd=Mod(rnorm(6)), mca=Mod(rnorm(6)))#######As medidas dos indivíduos mca2<-c(1.22222,1.666666,0.4500000,NA ,2.33333,0.789000) dados.pop.3<-cbind(ind.pop.3,med.pop.3,mca2)#######Planilha com os indivíduos e suas medidas head(dados.pop.3) hotelling = function(y1, y2) { ### Função que faz o cálculo de Hotellings by http://sas-and-r.blogspot.com/2010/05/example-737-calculation-of-hotellings. html # check the appropriate dimension k = ncol(y1) k2 = ncol(y2) if (k2!=k) stop("input matrices must have the same number of columns: y1 has ", k, " while y2 has ", k2) # calculate sample size and observed means n1 = nrow(y1) n2 = nrow(y2) ybar1= apply(y1, 2, mean); ybar2= apply(y2, 2, mean) diffbar = ybar1-ybar2 # calculate the variance of the difference in means v = ((n1-1)*var(y1)+ (n2-1)*var(y2)) /(n1+n2-2) # calculate the test statistic and associated quantities t2 = n1*n2*diffbar%*%solve(v)%*%diffbar/(n1+n2) f = (n1+n2-k-1)*t2/((n1+n2-2)*k) pvalue = 1-pf(f, k, n1+n2-k-1) # return the list of results return(list(pvalue=pvalue, f=f, t2=t2, diff=diffbar)) } ## Fazendo o teste entre o grupo 1 e 3 (result= hotelling(dados.pop.1,dados.pop.3)) Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29

1 0

Gráficos de superfície!!
by Thiago De paula protásio 15 Abr '11

15 Abr '11

Pessoal; Estou tendo dificuldades em colocar no gráfico de superfície tonalidades de cor diferenciada em cada parte do gráfico, pois isso facilita a visualização e interpretação da relação de cada variável independente com a variável independente. Estou usando os seguintes comandos: TCF=seq(68,80) RG=seq(32,37) f1=function(x,y)3696.046+62.183*x-33.414*y PCS=outer(TCF,RG,f1) persp(TCF,RG,PCS,theta=-50,phi=-0,shade=2,ticktype = "simple",box=T,border=T,r = sqrt(20), d = 1,col="blue") Agradeço a ajuda! Atenciosamente, 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

Data - time Zones
by eder＠leg.ufpr.br 15 Abr '11

15 Abr '11

Bom dia pessoal, Eu tenho datas no seguinte formato: DATE <- "00Z29Sep2006" # hora Z dia mês ano # Eu consegui mudar o formato da data com o comando: Sys.setlocale("LC_TIME", "USA") DATE <- strptime(paste(substr(DATE,4,12),substr(DATE,1,2),sep=''),"%d%B%Y%H") Sys.setlocale("LC_TIME","Portuguese_Brazil.1252")## Volta ao original DATE Meu objetivo final e ter a data no formato dd-mm-aaaa hh:mm:ss ('13-04-2011 08:35:00') Alguma tem uma dica, de como mudar o formato desta data direto sem ter que usar o Sys.setlocate? Atenciosamente Éder

5 10

Dúvida sobre função
by Victor Eduardo 15 Abr '11

15 Abr '11

Pessoal, há a possibilidade de criar uma única função que tenha varias váriaveis sendo que depois dessa função eu queria que cada variavel analisada da função gerasse uma coluna para cada uma, por exemplo, os comandos seriam mais ou menos esses: # Diretório wrk.dir <- "/home/nathalia/Área de Trabalho/synop/" dat.dir <- "/home/nathalia/Área de Trabalho/synop/" setwd(wrk.dir) # Leitura dos dados coluna<-c("sid", "ano", "mes", "dia", "hora", "hh", "lon", "lat", "height", "hls", "ps", "pmsl", "dd", "ff", "t", "td") dados<-read.table('201103.txt',head=T,col.names=coluna) attach(dados) # Criar uma nova coluna para a planilha dados #dados<-transform(dados, temp=t-273) fx<-function(var.temp, var.pre){ ifelse(var.temp> -25 & var.temp <40 | var.pre>930 & var.pre<1050, 1, 0) } dados<-transform(dados, CQ_temp= fx(t-273)) attach(dados) No caso, está errado a minha função, pois estou usando o | e isso não vai ficar legal, pois acabam que tornam elas dependentes, e quero analisar de forma independente. Agora na hora do transform, eu queria ver se dava para pegar cada variavel da minha função e ir adionando os 0 e 1 em coluna das respectivas variaveis. É possível isso? Abraços!

2 1

Arquivo Binario
by eder＠leg.ufpr.br 14 Abr '11

14 Abr '11

Pessoal, Minha duvida agora e com arquivos binarios, contextualizando: Quero abrir arquivos de saida de modelos climáticos (Eta), estes modelos soltam arquivo no formato GRIB, em em palavras R-anas uma array onde cada variável meteorológica seria um z dento do array (dados[x,y,z]), no arquivo ctl (indica como o binario é) indica qual o código falha (NA), em arquivo que o código era -9999 a função readBin funcionou, em outro arquivo que o código era 1.0E20 ela não funcionou ou seja estou suspeitando que ela não esta entendo este código falha, exemplo: formato <- readLines('http://www.leg.ufpr.br/~eder/ll20gfsa2006120100+2006120100.ctl') formato cone <- file('http://www.leg.ufpr.br/~eder/ll20gfsa2006120100+2006120100.dat','rb') dados <- readBin(cone,numeric(),195*200*58,4) length(dados) head(dados) close(cone) pelas informações em formato era para ter no arquivo binário 195*200*58= 2262000 números, porem ela carrega apenas 1911000, faltando dados, desta forma. Alguma dica? Atenciosamente

3 6

Ajustar escala em plot
by Marcelo Claro de Souza 14 Abr '11

14 Abr '11

Estou plotando alguns gráficos de regressões lineares da seguinte forma: cor=read.csv2('correlações.csv',h=TRUE) cor=as.data.frame(cor) simbolos=c(19,19,19,19,19,19,19,2,2,2,2,2,2,2,3,3,3,3,3,3,3) N.P.cor=lm(N~P,data=cor) plot(N~P,data=cor,pch=simbolos,ylab=expression(N~(g~Kg^{-1})),xlab=expression(P~(g~Kg^{-1}))) abline(N.P.cor) legend("topleft",c("Ep","Es","Et"),pch=c(19,2,3)) Como eu faço para delimitar a escala dos eixos x e y deixando por exemplo x(0:2) e y(0:30) Muito obrigado Marcelo Claro de Souza Biólogo, Doutorando em Biologia Vegetal, Depto. de botânica, IBRC- UNESP.16 - 91914780 http://lattes.cnpq.br/5691388412841376

3 2

Res: configuracoes do RStudio
by Samuel Carvalho 14 Abr '11

14 Abr '11

Obrigado Lucas. O problema era mesmo a versão do RStudio que estava utilizando, tanto que foi discutido isso no blog do software Benilton o problema não era como sugerido na lista Valeu ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 1

Sugestões
by Alexandre Galvão Patriota 14 Abr '11

14 Abr '11

Caros, Gostaria de parabenizar aos organizadores pela nova lista. Só tenho alguma sugestão: Existem muitas informações, o que é importante para os leitores. Porém, algumas informações importantes poderiam ser mais diretas e de fácil acesso sem precisar ler todo o conteúdo. Por exemplo, a senha de acesso. Está no meio do texto e uma batida de olho rápida eu não consegui acha-la. Só depois de ler o conteúdo quase inteiro do email eu percebi qual era a minha senha. Poderia ser uma coisa direta: Login: XXX Senha: YYY Quem não quiser ler o texto todo não precisará. Isso se aplica a outras informações importantes também. Fica a dica, Abraço a todos, Alexandre.

1 0

Intercept
by Edson Lira 14 Abr '11

14 Abr '11

Boa tarde a todos. Como faço para retirar o intercep do modelo de regressão, não consegui ver isso na ajuda online. Edson Lira Estatístico Manaus-Amazonas

3 2

configuracoes do RStudio
by Samuel Carvalho 13 Abr '11

13 Abr '11

Caros Alguem saberia me dizer como modificar as alterações no RStudio para que o mesmo aceite por exemplo acentos (^, ~) ou ainda cedilha (ç)? na leitura/importação dos dados... Valeu ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

3 2

Re: [R-br] Desdobramento da interação via glht:multcomp
by Ivan Bezerra Allaman 13 Abr '11

13 Abr '11

Só por curiosidade, por qual motivo utilizar a função glht? Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

teste de hipoteses
by Daniela Recchia 13 Abr '11

13 Abr '11

Boa tarde pessoal, desejo fazer um teste de hipoteses bem simples, porem nao consigo visualizar como deveria fazer isso no R, uma vez que os testes que ja fiz sempre foram utilizando um modelo (lm) e a funcao linearHypothesis. Nesse caso tenho os dados da seguinte forma: > cov V1 V2 V3 V4 V5 V6 45 65 81 55 31 0 41 80 73 51 36 0 42 79 74 40 45 0 43 77 76 62 * 0 Onde as colunas V1...V6 sao niveis de agua adicionados a uma planta, e cada linha eh uma caixa onde o experimento foi realizado, nesse caso 4 caixas para cada nivel de agua. Os valores sao o numero de semestes que germinaram. Quero testar a hipotese de que o nivel de agua nao influencia na germinacao das sementes. A principio pensei em usar anova, mas nao tenho claro como devo escrever meu modelo. Alguem poderia me dar uma dica? Abracos e bom dia! Daniela -- Daniela Rodrigues Recchia Master Student of Statistics - Technische Universität Dortmund. "Not everything that can be counted counts; and not everything that counts can be counted.” George Gallup

2 1

Re: [R-br] Res: Blocos incompletos no R
by Julia Maria 13 Abr '11

13 Abr '11

Walmes obrigado pela ajuda. A proposito, nao sei se deu tempo, mas naquele topico antigo vc pediu pra eu lhe mandar um .xls com os dados, ta lembrado? Voce chegou a dar uma olhada? Estou perguntando porque nao sei se recebi alguma mensagem sua, o modo com que essas mensagens chegam aqui estao meio bagunçadas. POr exemplo essa duvida que vc me respondeu hoje, eu usei o topico que eu criei la atras pra responder, mas ela apareceu sozinha na lista.

1 0

Blocos incompletos no R
by julia_est 13 Abr '11

13 Abr '11

Pessoal, tomei conhecimento dessa lista através do site do leg. Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III... Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai. Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov. Alguém poderia me ajudar? Grata Julia -- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407… Sent from the R-br mailing list archive at Nabble.com.

7 12

loop em nchar
by Daniel Marcelino 13 Abr '11

13 Abr '11

Caros, etou tentando obter o número de caracteres para cada célula para descobrir se trata-se de CNPJ ou CPF, mas acho que meu looping está com problemas pois não para de computar e não me retorna nada. O que está errado? for (i in CPF_CNPJ) { receitas2002$len = nchar(CPF_CNPJ) } Daniel

3 3

Importação correta de palavras acentuadas com a gdata::read.xls()
by Walmes Zeviani 13 Abr '11

13 Abr '11

Saudações, Tem pouco tempo que importo arquivos de extensão *xls. Estou com um problema que ao importar um arquivo possuo nomes de cidades brasileiras com acentos no nomes. Ao importar, as letras acentuadas são representadas de outra forma. Quero saber como importar corretamente ou se é possível devolver os nomes certos após importação. Veja o CMR que preparei require(gdata) # contém a função read.xls() para importar planilhas teste <- read.xls("http://www.leg.ufpr.br/~walmes/foo/teste.xls", sheet=1, method="tab") head(teste) # abaixo o print desse comando Sigla Codigo mun a2000 1 AC 1200013 Acrel\xe2ndia 1042.38 # â trocado por \xe2 2 AC 1200054 Assis Brasil 429.22 3 AC 1200104 Brasil\xe9ia 2243.25 4 AC 1200138 Bujari 786.68 5 AC 1200179 Capixaba 603.84 6 AC 1200203 Cruzeiro do Sul 9363.32 Os demais acentos possuem o seu "\xe?" correspondente. Alguém tem uma solução? Outro pacote/função? Meu sistema é Ubuntu 10.04. Os dados são um subconjunto dos dados do ipeadata, baixados no formato.xls. Outra coisa que percebi é que se o nome das colunas tiver acentos, a importação falha, ex "municipio" (sem acento) importa, mas com "município" (com acento) no nome da coluna não importa. Como evitar isso sem ter que sair renomeando as colunas manualmente? Planilhas de sites coleta de dados do Brasil sempre tem acentos por causa do nomo idioma. Minhas buscas no nabble não retornaram resultado. # sheet 2 está escrito "Município" como nome da coluna teste <- read.xls("http://www.leg.ufpr.br/~walmes/foo/teste.xls", sheet=2, method="tab") Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 3 Certo sua atenção, agradeç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 ==========================================================================

8 12

Anova multivariada.
by Edson Lira 12 Abr '11

12 Abr '11

Boa noite a todos. Estou replicando um exemplo de um livro que verifica a dependência entre o tempo de viagem(y) ~ número de entregas(x1) +número de milhas percorridas. Quando peço a anova no R. anova(mod.cam), o R me dá a saída abaixo. Analysis of Variance Table Response: cam$tem_via Df Sum Sq Mean Sq F value Pr(>F) cam$n_entrega 1 9.0449 9.0449 27.535 0.001190 ** cam$milhas 1 12.5556 12.5556 38.222 0.000453 *** Residuals 7 2.2994 0.3285 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Eu gostaria da ANOVA com o Df da Regressão e do Resíduo, ou seja a anova reduzida. Alguém poderia me indicar qual a função ou rotina para isso? Em anexo o banco de dados em csv e a rotina. []' s Edson Lira Estatístico Manaus-Amazonas

2 2

estimadores
by Killa p'unchay Morales 12 Abr '11

12 Abr '11

--- Em ter, 12/4/11, Killa p'unchay Morales <killa.punchay(a)yahoo.com.br> escreveu: De: Killa p'unchay Morales <killa.punchay(a)yahoo.com.br> Assunto: help Para: "RRRR" <r-br-request(a)listas.c3sl.ufpr.br> Data: Terça-feira, 12 de Abril de 2011, 3:04 olá amigos estou com dificultade de encontrar os 4 parâmetros desta função de verossimilhança....primeiro to tentando derivar uma somatoria da acumulada da Gama que chamei de XX com respeito de alpha e derivar a somatoria do log da acumulada da Gama com respeito de alpha = YY, fez analiticamente as derivadas com respeito dos 4 parâmetros......não sei utilizar o optim para maximizar essa função. estou anexando o que fez. Obrigada alguem me ajude.....

2 1

Problemas write.table
by Alan Rodrigo Panosso 12 Abr '11

12 Abr '11

Caros amigos, Como salvar a lista gerada pelo código abaixo, utilizando a função write.table? Estou encontrando o erro abaixo! set.seed(200) dados<-rnorm(8836,36,6) NL=94 NC=94 Xmin=0 Xmax=64 Ymin=0 Ymax=64 Zmin=min(dados) Zmax=max(dados) ma<-matrix(dados,ncol=NC) mc<-matrix(c(NL, Xmin,Ymin,Zmin,NC,Xmax,Ymax,Zmax),ncol=2) lac<-list(mc,ma) write.table(lac,"C://r//Saida.txt",sep="\t",col.names=F,row.names=F) Erro em data.frame(c(94, 0, 0, 12.7834627990016, 94, 64, 64, 61.968478487737 : arguments imply differing number of rows: 4, 94 Abraços a todos. Alan Rodrigo Panosso

3 6

Correlação de Spearman - "Não é possível calcular valores p exatos quando há empates"
by Edson Ferraz Evaristo de Paula 12 Abr '11

12 Abr '11

Boa noite! O que seria essa "warning message" que aparece quando rodo cor.test (pelo método de spearman)? "Não é possível calcular valores p exatos quando há empates" Obrigado -- Edson F. Evaristo de Paula (41) 9669-3646 (42) 8813-7268

3 2

Como responder a um tópico da lista!
by Ivan Bezerra Allaman 11 Abr '11

11 Abr '11

Boa noite senhores! Me parece que muitos (inclusive eu) estavam tendo problemas ao responder um determinado tópico da lista, gerando inúmeros tópicos repetidos. Eu particulamente, não estou a contento da maneira que está, em só responder os tópicos por email, quando poderíamos responder tais tópicos diretamente pela lista, como acontece no R internacional. Mais já foi justificado o porque disto. Logo, para evitarmos estes inúmeros tópicos repetidos, o que polui a lista, ao responder tal tópico, basta fazer o seguinte: Suponhamos que eu queira responder o tópico: [R-br] Anova duas vias com medidas repetitivasBasta enviarmos para a lista desta maneira: Para: r-br(a)listas.c3sl.ufpr.br Cc: Assunto: Re: [R-br] Anova duas vias com medidas repetitivas Pronto, basta colocar "Re:" na frente do título do tópico que sua resposta irá para o mesmo tópico da pessoa que a enviou, evitando deste modo outro tópico com o mesmo nome. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Anova duas vias com medidas repetitivas
by Ivan Bezerra Allaman 11 Abr '11

11 Abr '11

Chute doisssss! M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Anova duas vias com medidas repetitivas
by Ivan Bezerra Allaman 11 Abr '11

11 Abr '11

Chute um! M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Converter Shapefile em arquivo do pacote maps
by Edson Lira 11 Abr '11

11 Abr '11

Antonio, veja se o material em anexo te ajuda de alguma forma. Também sou de Manaus, talvez possamos trocar algumas idéias sobre o R. Meu email pessoal é edinhoestat(a)yahoo.com.br Edson Lira Estatístico Manaus-Amazonas --- Em seg, 11/4/11, Antônio Silva <toinho_sb(a)yahoo.com.br> escreveu: De: Antônio Silva <toinho_sb(a)yahoo.com.br> Assunto: [R-br] Converter Shapefile em arquivo do pacote maps Para: r-br(a)listas.c3sl.ufpr.br Data: Segunda-feira, 11 de Abril de 2011, 12:26 Bom dia! Estou buscando uma forma de plotar em 3D o mapa do Amazonas que está no formato shape.file juntamente com outra informação tipo densidade demografica por municipio, onde seja possivel girar o mapa a tal ponto de ser possivel visualizar o grafico de barras verticais em 3D, com a densidade, em cima do mapa. Encontrei um exemplo onde o autor utiliza o mapa dos EUA, que é oriundo do pacote maps, a diferença é que o mapa dos EUA já estava no pacote maps. Como faço para converter o arquivo shape.file no formato que o pacote maps leia? Att, Antonio Alcirley -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

1 0

Erro em Redes Neurais
by Alan Rodrigo Panosso 11 Abr '11

11 Abr '11

Caros amigos, estamos utilizando o pacote AMORE do R, e gostaríamos de saber como incluir um novo neurônio numa camada. Inserimos a função Gamma como ativação num novo neurônio e veio o seguinte erro na função train: Integer() can only be applied to a integer, not a double. Muito obrigado a todos Alan Rodrigo Panosso

2 1

Comando para cada alteração na variável calcular média
by Alexandre dos Santos 11 Abr '11

11 Abr '11

Bom Dia Pessoal, Fiz vários testes de Kolmogorov-Smirnov e o resultado é mostrado abaixo, para valor na coluna raio tenho 99 simulações (Nsim) e o resultado do teste p.value, gostaria de tirar a média de p.value a cada mudança na coluna raio, sem ter que especificar raio por raio como por exemplo: raio.0.025=res.ks[res.ks[,2]==0.025,], uma vez que tenho 40 diferentes raios, alguem poderia me dar uma luz, Obrigado, Alexandre #Criando uma tabela de resultados > write.table(res.ks, file ="res.ks", row.names=F, col.names=F) > res.ks=read.table("res.ks",header=F) > colnames(res.ks)=c("p.value","raio","Nsim") > head(c) p.value raio Nsim 1 3.671638e-02 0.025 1 2 2.396494e-06 0.025 2 3 3.053739e-07 0.025 3 4 7.559981e-05 0.025 4 5 1.147039e-04 0.025 5 6 1.071125e-01 0.025 6 > #

3 4

RES.: Função runs.test!! Modelo de médias!
by Thiago De paula protásio 11 Abr '11

11 Abr '11

Prezados colegas da R-Br; Escrevo-lhes pois estou utilizando o teste de falta de ajuste indica pelo Walmes e no qual estou comparando o meu modelo não linear com um modelo de médias. Mas infelizmente não estou conseguindo um bom ajuste quando comparo o meu modelo com o de médias. Alguém poderia me indicar uma maneira de contornar esse problema? O Walmes já esclareceu muita coisa.......mas gostaria de saber se existe um outro teste o que posso fazer? Ou algo semelhante? Os dados referentes ao ajuste e os comandos indicados pelo Walmes estão sendo de grande valia e encontram-se nas postagens anteriores. Caso necessitem novamente dos comandos e dados ...me avisem! Grato pela preciosa ajuda!!! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Anova duas vias com medidas repetitivas
by Daniela Faria Florencio 11 Abr '11

11 Abr '11

Boa tade, estou precisando realizar uma análise com observações repetidas ao longo de 3 tempos e em quatro diferentes tratamentos. Além disso, os dados são de proporção. Como posso realizar esta análise no R??? Obrigada -- Daniela Faria Florencio Universidade Federal de Viçosa Pós-doutorado em Entomologia

3 2

Anova duas vias com medidas repetitivas
by Ivan Bezerra Allaman 11 Abr '11

11 Abr '11

Veja se uma abordagem via Modelos Lineares Generalizados não seria mais interessante. Há vários materiais na net sobre o assunto. Caso contrário, é possível fazer a análise via LM, se assegurando é claro, que os resíduos atendem as pressuposições. ANOVA no R é simples. Busque por "Paulo Justiniano Ribeiro Júnior" no google, que na página dele há vários materias sobre o assunto. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

trabalhando na area de trabalho do R - part2
by Samuel Carvalho 11 Abr '11

11 Abr '11

Caros Venho novamente recorrer sobre um questionamento de como se trabalhar com data.frames diferentes na area de trabalho do R e poder manipula-los como o exemplo a seguir ### Código R dados <- data.frame(parcela=rep(1:5, each=3), rep=rep(1:3, 5)) dados$resposta <- rnorm(nrow(dados)) dados tmp <- dados[dados$rep == 3,] tmp dados$resposta[which(dados$parcela %in% tmp$parcela)] <- tmp$resposta dados ### A idéia é repetir o valor da resposta de rep=3 para todos os outros valores de rep daquela parcela ### Ex: Se na parcela 1 o valor da resposta na rep=3 seja -0.8543 então os valores na rep 1 e 2 ficariam ### também como -0.8543. Porém neste código acima ele repete a ultima medição sequencialmente nao levando ### em consideração se a parcela muda ou nao. ### Para os adeptos de planilha seria algo semelhante ao "procv" []'s Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

1 0

RStudio Beta 2
by Benilton Carvalho 11 Abr '11

11 Abr '11

Aos entusiastas: http://blog.rstudio.org/2011/04/11/rstudio-beta2/ b

1 0

Converter Shapefile em arquivo do pacote maps
by Antônio Silva 11 Abr '11

11 Abr '11

Bom dia! Estou buscando uma forma de plotar em 3D o mapa do Amazonas que está no formato shape.file juntamente com outra informação tipo densidade demografica por municipio, onde seja possivel girar o mapa a tal ponto de ser possivel visualizar o grafico de barras verticais em 3D, com a densidade, em cima do mapa. Encontrei um exemplo onde o autor utiliza o mapa dos EUA, que é oriundo do pacote maps, a diferença é que o mapa dos EUA já estava no pacote maps. Como faço para converter o arquivo shape.file no formato que o pacote maps leia? Att, Antonio Alcirley

1 0

Re: [R-br] Simbolos pré-definidos em função das linhas de uma tabela
by Marcelo Claro de Souza 11 Abr '11

11 Abr '11

Olá, infelizmente não deu certo. Continua aparecendo o nome da linha na figura ao invés de aparecer os símbolos. Eu tentei os comandos sugeridos tanto p biplot quanto p plot e mesmo assim não deu certo. Tem alguma outra forma que eu possa tentar? Em anexo segue o conjunto de dados que estou utilizando para vcs fazerem alguns testes. Muito obrigado. Marcelo Claro de Souza Biólogo, Doutorando em Biologia Vegetal, Depto. de botânica, IBRC- UNESP.16 - 91914780 http://lattes.cnpq.br/5691388412841376 --- Em sex, 8/4/11, Walmes Zeviani <walmeszeviani(a)gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani(a)gmail.com> Assunto: Re: [R-br] Simbolos pré-definidos em função das linhas de uma tabela Para: r-br(a)listas.c3sl.ufpr.br Data: Sexta-feira, 8 de Abril de 2011, 0:42 Marcelo, Seria isso biplot(dados.pca, xlab="pc1", ylab="pc2", pch=rep(c(19,1,2), each=6)) Espero que dê certo. À 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 ========================================================================== Em 7 de abril de 2011 23:12, Augusto Ribas <ribas.aca(a)gmail.com> escreveu: Tipo, um numero corresponde um simbulo. Por exemplo 19 é bolinha preta fechada 1 é bolinha com contorno preto e interior branco. digita ?pch ou ?points que vai ter uma parte que fala qual numero faz cada simbulo. então segundo seu exemplo vc deveria fazer algo como: biplot(dados.pca,xlab="pc1",ylab="pc2", pch(rep(19,6),rep(1,6),rep(2,6)) ai no caso repete o valor 19 que é bolinha preta 6 vez que é sua primeira especies e assim vai. a segunda especie vai ser bolinha aberta, que é o numero 1. O vertor tem que ter o mesmo tamanho do numero de pontos que vai ter no grafico, numero de linhas no caso. Acho que é mais ou menos isso :) Em 7 de abril de 2011 22:01, Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> escreveu: Walmes, Eu sou bem novato em R e não entendi direito. No caso eu devo usar da seguinte forma por exemplo: biplot(dados.pca,xlab="pc1",ylab="pc2", pch(1:7,8:14,15:21)) é mais ou menos isso que tenho de fazer? Muito obrigado. Marcelo Claro de Souza Biólogo, Doutorando em Biologia Vegetal, Depto. de botânica, IBRC- UNESP. 16 - 91914780 http://lattes.cnpq.br/5691388412841376 --- Em qui, 7/4/11, Walmes Zeviani <walmeszeviani(a)gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani(a)gmail.com> Assunto: Re: [R-br] Simbolos pré-definidos em função das linhas de uma tabela Para: r-br(a)listas.c3sl.ufpr.br Data: Quinta-feira, 7 de Abril de 2011, 14:51 Marcelo, O controle dos símbolos é com o argumento pch= {p}oint {ch}aracter dentro das funções gráficas. Em gráficos de dispersão, basta passar um vetor com os índices que identificam os símbolos. Para conhecer os símbolos faça plot(1:30, pch=1:30) Espero que no biplot() seja assim também. Caso não dê certo, forneça um exemplo reproduzível e destaque o erro para que possamos ajudar. À 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 ========================================================================== Em 7 de abril de 2011 13:44, Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> escreveu: Possuo um conjunto de dados composto por 21 linhas e 3 colunas. Ao importar os dados eu utilizo o comando row.names e atribuo nome as linhas para facilitar a visualização das espécies na pca (no caso utilizo o pacote vegan e uso o comando rda(dados,scale=TRUE) nesse conjunto de dados da linha 1 a 7 são plantas da mesma espécie, de 8 a 14 outra e de 16 a 21 a outra. Como eu faço p atribuir um símbolo único para cada espécie (ex: quadrado, circulo e triângulo) ao plotar um biplot? Muito obrigado. Marcelo Claro de Souza Biólogo, Doutorando em Biologia Vegetal, Depto. de botânica, IBRC- UNESP. 16 - 91914780 http://lattes.cnpq.br/5691388412841376 _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056 _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

3 3

Escrever valores sub e sobrescritos em legendas
by Marcelo Claro de Souza 11 Abr '11

11 Abr '11

Como eu faço p indicar em uma legenda que determinado número deve se sub ou sobrescrito para ficar na seguinte forma xlab="N (mg.Kg-1)" ao invés de xlab="N (mg.Kg-1)" e CO2 ao invés de CO2? Muito obrigado Marcelo Claro de Souza Biólogo, Doutorando em Biologia Vegetal, Depto. de botânica, IBRC- UNESP.16 - 91914780 http://lattes.cnpq.br/5691388412841376

3 3

pacote DTK, teste Tukey-Kramer...
by Edson Ferraz Evaristo de Paula 11 Abr '11

11 Abr '11

Bom dia pessoal! Alguem pode me ajudar com o pacote DTK? Tenho um experimento com tratamentos que possuem nº diferentes de repetições (T1=8, T2=8 e T3=4), queria fazer o teste de Tukey-Kramer mas não to me acertando com o comando DTK.test(x, f, a) e o gl.unequal(n, k)... Gostaria de saber como entro com os dados nesse caso. Obrigado!! -- Edson F. Evaristo de Paula (41) 9669-3646 (42) 8813-7268

2 1

RES.: Função sigmóide! Função nls.lm
by Thiago De paula protásio 10 Abr '11

10 Abr '11

Allaman e demais membros da lista; Obrigado pela dica da função nls.lm(). Vc pode nos dar mais detalhes de como podemos proceder para ajustar um modelo não linear a partir dessa função? Se vc puder dar uma olhada nos dados.......eu já mandei para a lista os dados e o Walmes ele fez a gentileza de anexá-la. Muito obrigado!!!! Bom final de semana a todos!

3 2

Aplicademic - Nova postagem
by Vinicius Brito Rocha 10 Abr '11

10 Abr '11

Em fim uma nova postagem em www.aplicademic.blogspot.com Att. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

1 0

RES.: Função runs.test! Falta de Ajuste!!!!
by Thiago De paula protásio 10 Abr '11

10 Abr '11

Walmes e demais membros da r-br; Muito obrigado pela ajuda!! Suas dicas e ensinamentos estão sendo muito válidos..... Vou estudar com calma o porquê dessa falta de ajuste e tentar ajustar outros modelos para comparar.... Calculei aqui o R² do modelo que vc sugeriu e foi bem superior a o modelo inicial....Acho que estamos no caminho certo. Já encontrei a sua dissertação e estou lendo a parte do teste da falta de ajuste. Vc tem outros trabalhos para me indicar? Gostaria de entender melhor esse teste e o modelo de médias..... Novamente, só tenho que agradecer a ajuda. Att., 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