- 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