Junho 2012 - R-br - DINF UFPR

Densidade normal multivariada em área irregular
by Rodrigo Sant'Ana 11 Jun '12

11 Jun '12

Boa noite R's, estou buscando, homericamente, gerar uma distribuição de densidades multivariada em uma área irregular, porém estou dando murro em ponta de faca. Acredito que devo estar errando em algo simples, e fácil. O intuito é simular a distribuição de densidade de um peixe em uma área conhecida, para que posteriormente, estas densidades sejam alvos de simulações amostrais. Abaixo segue o CMR de minha tentativa! ###---------------------------------- Simulação Densidade Normal Multivariada # Carregando os dados em anexo! load("R-BR.RData") ## A área de interesse deveria ser esta plotada a seguir plot(bor, type = "l") # Separação dos vetores Lat e Long do objeto dados long <- dados$long lat <- dados$lat ## Utilizando um kernel de Densidade da Normal Multivariada # Função de amostragem estratificada! system.time( d <- kde2d(long, lat, n = 100) ) ## Função criada para gerar as figuras em perspectiva (persp) e contorno ## (contour) de mesma densidade (biomassa). kde2dplot <- function(d, # Densidades 2D computada pela # função kde2d ncol=50, # Número de cores zlim=c(0,max(z)), # Limites da coordenada z nlevels=20, # vide nlevels em contour theta=30, # vide theta em persp phi=60) # vide phi em persp { z <- d$z # Recebe as densidades (biomassa) nrz <- nrow(z) # Recebe o número de linhas de z ncz <- ncol(z) # Recebe o número de colunas de z ## O objeto couleurs determina o esquema de cores a ser utilizado, sendo ## as opções: rainbow, heat.colors, terrain.colors, topo.colors, ## cm.colors. couleurs <- tail(terrain.colors(trunc(1.4 * ncol)), ncol) fcol <- couleurs[trunc(z/zlim[2]*(ncol-1)) + 1] dim(fcol) <- c(nrz, ncz) fcol <- fcol[-nrz, -ncz] par(mar = c(5, 5, 2, 2)) image(d, col = couleurs, xlab = "Longitude", ylab = "Latitude") contour(d, add = T, nlevels = nlevels, xlab = "Longitude", ylab = "Latitude") box() } ## Salva a figura de saída da simulação gerada kde2dplot(d) polygon(bor) ##----------------------------------------------------------------------------------------------------- Se alguém tiver uma alternativa, realmente salvaria uma alma! Brincadeiras a parte, queria muito agradece-los pela atenção. ______________ Rodrigo Sant'Ana Oceanógrafo

1 0

Ajuda para rodar um de dois códigos R
by Maria Papa 11 Jun '12

11 Jun '12

Olá Pessoal, bom dia, Eu tenho dois códigos R que eu seleciono em função da minha aplicação, ou seja, do tipo de dados, eu gostaria de saber de alguma função que permita rodar um entre estes dois códigos R. Obrigada a todos! Célia

2 1

Teste de aleatoriedade
by Gilbert Queiroz 11 Jun '12

11 Jun '12

Olá! Alguém poderia indicar algum teste de aletoriedade no R? Abs. Gilbert.

2 1

Dunnet
by Peterson Gardin 11 Jun '12

11 Jun '12

Pessoal Preciso fazer um teste de Dunnet para verificar se existe diferença entre os tratamentos (níveis 0,002, 0,004 e 0,006) e a testemunha (dose 0,000). Esse experimento é um aditivo colocado em um alimento, mas é preciso que os consumidores (julgadores) não percebam a diferença. Ou seja, não seja diferente da testemunha. Realize a análise, mas o tete de Dunnet sai comparando cada julgador, gostaria que a saída, me desse, cada tratamento em relação a testemunha que é a dose zero do banco em anexo. Usei o pacote multicomp para fazer a comparação. Desde já agradeço dados3 <- read.table("prova3.txt", h=T, dec=",") attach(dados3) names(dados3) levels(trat) trat <- factor(trat) prova3.av <- aov(resp~julga+trat) anova(prova3.av) require(multcomp) summary(glht(prova3.av)) João Peterson Pereira Gardin Doutor em Fisiologia Vegetal, Eng. Agr., Pesquisador Fones: (49) 3566-0054 e 9933-4290 MSN: petersongardin(a)yahoo.com.br - SKYPE: petersongardin Endereço: Rua João Zardo, 1660, Bairro Campo Experimental, Videira - SC - CEP 89560000

2 1

Usar o nome dos elementos da lista com sapply lapply
by Walmes Zeviani 09 Jun '12

09 Jun '12

Saudações, Estou fazendo análise descritiva de várias variáveis e gostaria de automatizar em parte o processo. Assim, uma lista contém cada variável como um elemento e eu gostaria de usar o nome do elemento da lista como título do gráfico. Veja o CMR mylist <- list(salario=sample(1:3,40,r=TRUE), dependentes=sample(1:3,40,r=TRUE)) par(mfrow=c(2,1)) sapply(mylist, function(x) barplot(table(x), main=deparse(substitute(x)))) for(i in names(mylist)) barplot(table(mylist[[i]]), main=i) Eu gostaria de pegar os nomes "salario" e "dependentes" e usá-los como titulo. Tentei várias combinações de parse(), deparse(), substitute(), etc. Quero evitar o for(). Alguma sugestão? Grato. 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 ==========================================================================

5 4

problemas no uso da função Wireframe
by Fernando Antonio de souza 08 Jun '12

08 Jun '12

Caros colegas, estou com dificuldades em utilizar a funçao wireframe para contruir um gráfico 3D. A função fica me retornando que não é permitido utilizar vetores de comprimento negativo. Entretanto não há (nem existe) vetores de comprimento negativo. Caso tenha experiência no uso desta função agradeço desde ja a ajuda GPD<-rep(100,66) FDN<-rep(15:80,each=1,1) dput(GRAF) CMSPM<- structure(c(41.4342238231054, 43.9158002951626, 46.3314684686358, 48.6812283435249, 50.9650799198301, 53.1830231975512, 55.3350581766883, 57.4211848572414, 59.4414032392104, 61.3957133225955, 63.2841151073965, 65.1066085936135, 66.8631937812465, 68.5538706702954, 70.1786392607604, 71.7374995526413, 73.2304515459382, 74.6574952406511, 76.01863063678, 77.3138577343248, 78.5431765332856, 79.7065870336625, 80.8040892354552, 81.835683138664, 82.8013687432888, 83.7011460493295, 84.5350150567862, 85.3029757656589, 86.0050281759476, 86.6411722876522, 87.2114081007729, 87.7157356153095, 88.1541548312621, 88.5266657486306, 88.8332683674152, 89.0739626876157, 89.2487487092323, 89.3576264322648, 89.4005958567133, 89.3776569825777, 89.2888098098582, 89.1340543385546, 88.913390568667, 88.6268185001954, 88.2743381331398, 87.8559494675001, 87.3716525032764, 86.8214472404687, 86.205333679077, 85.5233118191013, 84.7753816605415, 83.9615432033978, 83.08179644767, 82.1361413933581, 81.1245780404623, 80.0471063889825, 78.9037264389186, 77.6944381902707, 76.4192416430388, 75.0781367972229, 73.671123652823, 72.198202209839, 70.659372468271, 69.054634428119, 67.383988089383, 65.647433452063), label = "Predicted values") data<-data.frame(CMSPM=CMSPM,GPD=GPD,FDN=FDN) wireframe(CMSPM~FDN*GPD,data=data) structure(c(41 structure(c(41.4342238231054, 43.9158002951626, 46.3314684686358, 48.6812283435249, 50.9650799198301, 53.1830231975512, 55.3350581766883, 57.4211848572414, 59.4414032392104, 61.3957133225955, 63.2841151073965, 65.1066085936135, 66.8631937812465, 68.5538706702954, 70.1786392607604, 71.7374995526413, 73.2304515459382, 74.6574952406511, 76.01863063678, 77.3138577343248, 78.5431765332856, 79.7065870336625, 80.8040892354552, 81.835683138664, 82.8013687432888, 83.7011460493295, 84.5350150567862, 85.3029757656589, 86.0050281759476, 86.6411722876522, 87.2114081007729, 87.7157356153095, 88.1541548312621, 88.5266657486306, 88.8332683674152, 89.0739626876157, 89.2487487092323, 89.3576264322648, 89.4005958567133, 89.3776569825777, 89.2888098098582, 89.1340543385546, 88.913390568667, 88.6268185001954, 88.2743381331398, 87.8559494675001, 87.3716525032764, 86.8214472404687, 86.205333679077, 85.5233118191013, 84.7753816605415, 83.9615432033978, 83.08179644767, 82.1361413933581, 81.1245780404623, 80.0471063889825, 78.9037264389186, 77.6944381902707, 76.4192416430388, 75.0781367972229, 73.671123652823, 72.198202209839, 70.659372468271, 69.054634428119, 67.383988089383, 65.647433452063), label = "Predicted values") .4342238231054, 43.9158002951626, 46.3314684686358, 48.6812283435249, 50.9650799198301, 53.1830231975512, 55.3350581766883, 57.4211848572414, 59.4414032392104, 61.3957133225955, 63.2841151073965, 65.1066085936135, 66.8631937812465, 68.5538706702954, 70.1786392607604, 71.7374995526413, 73.2304515459382, 74.6574952406511, 76.01863063678, 77.3138577343248, 78.5431765332856, 79.7065870336625, 80.8040892354552, 81.835683138664, 82.8013687432888, 83.7011460493295, 84.5350150567862, 85.3029757656589, 86.0050281759476, 86.6411722876522, 87.2114081007729, 87.7157356153095, 88.1541548312621, 88.5266657486306, 88.8332683674152, 89.0739626876157, 89.2487487092323, 89.3576264322648, 89.4005958567133, 89.3776569825777, 89.2888098098582, 89.1340543385546, 88.913390568667, 88.6268185001954, 88.2743381331398, 87.8559494675001, 87.3716525032764, 86.8214472404687, 86.205333679077, 85.5233118191013, 84.7753816605415, 83.9615432033978, 83.08179644767, 82.1361413933581, 81.1245780404623, 80.0471063889825, 78.9037264389186, 77.6944381902707, 76.4192416430388, 75.0781367972229, 73.671123652823, 72.198202209839, 70.659372468271, 69.054634428119, 67.383988089383, 65.647433452063), label = "Predicted values")

2 2

4 plots no mesmo frame em alta resolução no linux
by Marcelo Claro de Souza 08 Jun '12

08 Jun '12

Preciso plotar 4 dendrogramas no mesmo frame e salvar em alta resolução, mas somente a última figura está sendo salva. Onde está meu erro? require(graphics) X11() par(mfrow=c(2,2)) wid = 1500 hei = 1200 reso = 300 tiff(file="dendrograma.tiff",width = wid, height=hei, res=reso) plot(1...) plot(2...) plot(3...) plot(4...) dev.off() Muito obrigado pela ajuda Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Aumentar dpi de figura
by Marcelo Claro de Souza 08 Jun '12

08 Jun '12

A resolução padrão dos geráficos do R é de 96 dpi, entretanto preciso salvar minhas figuras com 300 dpi. Como posso fazer isso? Abraço Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

4 3

Cálculo de coordenadas finais à partir da distância e ângulo
by ASANTOS 08 Jun '12

08 Jun '12

Boa tarde pessoal, Gostaria de calcular as coordenadas e representá-las à partir de um ponto inicial (p0) com as distâncias (dist) e ângulos (ang) conhecidos, sendo: ##Cálculo das coordenadas à partir da distância e ângulo em relação à p0 p0<-c(100,100)##Coordenada conhecida dist<-c(5,15,7.5,3.5)##Distância à partir de p0 ang<-c(256,93,183,15) Fiz algumas pesquisas e encontrei a função grid.curve() do pacote grid, que faz a representação gráfica das distâncias(retas) dados o ponto inicial (p0), pontos finais (x,y) e o ângulo (ang), meu problema é que as coordenadas finais são as variáveis que desejo encontrar, para depois fazer a representação gráfica Alguém teria alguma sugestão de por onde eu poderia começar? Obrigado -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

3 3

Re: [R-br] Aumentar dpi de figura
by Ivan Bezerra Allaman 07 Jun '12

07 Jun '12

Há uma longa discussão aqui: http://r.789695.n4.nabble.com/300-dpi-and-eps-td3088001.html (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0