- R-br - DINF UFPR

[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 16 Abr '11

16 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