Junho 2011 - R-br - DINF UFPR

Plotar dois tipos de gráficos em uma mesma janela gráfica com a lattice.
by Ivan Bezerra Allaman 15 Jun '11

15 Jun '11

Bom dia senhores! Estou tentando plotar dois tipos de gráficos em uma mesma janela gráfica, mais tá difícil. Não sei até onde vale a pena o 'tempo/benefício' em se usar a lattice. Segue abaixo o CMR. require(plotrix) jet.colors <- colorRampPalette( c("blue", "green") ) x <- seq(-1.95, 1.95, length=30) y <- seq(-1.95, 1.95, length=35) da <- expand.grid(x=x, y=y) da$z <- with(da, x*y^2) require(lattice) panel.3d.contour <- function(x, y, z, rot.mat, distance, nlevels = 20, zlim.scaled, ...) { add.line <- trellis.par.get("add.line") panel.3dwire(x, y, z, rot.mat, distance, zlim.scaled = zlim.scaled, ...) clines <- contourLines(x, y, matrix(z, nrow = length(x), byrow = TRUE), nlevels = nlevels) for (ll in clines) { m <- ltransform3dto3d(rbind(ll$x, ll$y, zlim.scaled[1]), rot.mat, distance) panel.lines(m[1,], m[2,], col = add.line$col, lty = add.line$lty, lwd = add.line$lwd) } } g1 <- wireframe(z~x+y, da, drape=TRUE, scales=list(arrows=FALSE),col.regions=jet.colors(100),panel.3d.wireframe="panel.3d.contour") g2 <- contourplot(z ~ x+y,da) library(grid) pushViewport(viewport(layout=grid.layout(2,1))) print(g1,vp=viewport(layout.pos.row=1,layout.pos.col=1)) print(g2,vp=viewport(layout.pos.row=2,layout.pos.col=1)) Este último comando não plota os dois gráficos em uma mesma janela gráfica. Eu gostaria de um comando equivalente ao par(mfrow). Desde já agradeço novamente a atenção de todos. 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

Funcionamento (implementação) da função apply (lapply) em loops
by Junior Beleti 14 Jun '11

14 Jun '11

Boa noite, gostaria de entender o funcionamento da função apply, mais precisamente lapply, para loops. Já conheço a descrição de seu funcionamento: ‘lapply’ returns a list of the same length as ‘X’, each element of which is the result of applying ‘FUN’ to the corresponding element of ‘X’. Mas o que eu gostaria de saber é como ele funciona em sua essência. Por exemplo, utilizar um lapply ao invés de um loop "for" é mais rápido. Eu gostaria de saber o porque, ou seja, como o lapply é implementado. Não sei se fui claro, mas gostaria de encontrar a implementação da função lapply, o que não encontrei no diretório do pacote "base". Se alguém puder ajudar. Também peço desculpas por não ter dado grandes contribuições para a lista, mas creio ter mais experiência em programação (lógica) de fato, do que no ambiente R. Obrigado. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

5 5

krige.conv() do geoR
by Alexandre dos Santos 14 Jun '11

14 Jun '11

Boa noite a todos, Transformei minha variável resposta (data) do geoR em log(x+1), mas na hora de fazer a krigagem desejo que o $predict da krige.conv() fosse com os valores da minha variável sem tranformação, então gostaria de saber se tem alguma forma de fazer exp(x)+1 no $predict. Uma outra pergunta também, tem como fazer interpolação pelo inverso do quadrado da distância com o krige.control(), Obrigado, 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

3 4

número mínimo de casos em glm e betafit
by Daniel Marcelino 14 Jun '11

14 Jun '11

Caros, alguém saberia me informar se existe algum limite mínimo de casos para uma análise de glm, gls e betafit? Não achei nenhuma informação que me levasse diretamente a este ponto.

5 9

erro padrão no grafico
by Gilson Sanchez 13 Jun '11

13 Jun '11

Pessol mais uma vez consultando vc's Estou querendo fazer um gráfico de linhas, e em cada ponto quero colocar o erro padrão (sd), alguem pode ajudar O CMR que estou trabalhando require(XLConnect) especie1<-loadWorkbook("foto.xls", create= FALSE) especie1<-readWorksheet(folha.numero1, sheet="especie1") especie1 attach(especie1) plot(par,foto, type= "p", ylim=c(-2,10), xlim=c(0,2000), pch=16) abline(h=0, lty=4) lines(par,hib) dados: http://www.datafilehost.com/download-94c764d9.html 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.

2 1

Converter datas em nº de dias do ano
by Cristiano Santos 13 Jun '11

13 Jun '11

Caros, Sou novo no R e não estou conseguindo manipular datas. Tenho uma variável com data e hora, e quero transformar no número de dias do ano, por exemplo: "4/02/2011 03:24:13" deverá transformar-se em "35" Desde já agradeço pela ajuda. Cristiano Santos Economia-UFC

5 8

geoestatística - ajuste de modelo
by nani.design＠bol.com.br 13 Jun '11

13 Jun '11

2 2

Re: [R-br] Atualização do R no Ubuntu 11.04
by Junior Beleti 13 Jun '11

13 Jun '11

Instalei o aptitude e subversion, e após consegui instalar o R. Muito obrigado... abraços -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

1 0

Re: [R-br] Eliminar NA e merge com colunas de dif tamanho
by Alexandre Heinemann 12 Jun '11

12 Jun '11

Olá Eder, Muito obrigado. Alexandre

1 0

Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo)
by Marcelo Cardoso mello 12 Jun '11

12 Jun '11

Estou com dificuldades em desdobrar a interação em uma regressão em fatorial. Segue um exemplo abaixo. Como obter um modelo de regressão de r em função de b para a1 e a2?????? a<-as.factor(c(rep(1,9),rep(2,9))) b<-c(rep(0,3),rep(5,3),rep(10,3), rep(0,3),rep(5,3),rep(10,3)) r<-c(1,1.2,1.3,2,2.3,2.25,4,3.9,3.8,5,4.8,4.7,3.5,3,3.4,1.7,1.5,1.2) b<-as.factor(b) m<-lm(r~a*b) anova(m) b<-as.numeric(b) regressão<-lm(r~b) regressão anova(regressão) Como desdobrar b dentro de a1 e a2????? Obrigado. Marcelo.

4 9