- R-br - DINF UFPR

Ajuda em uma função
by Marcelo Zeuli 15 Jun '12

15 Jun '12

Bom dia, Algume conheçe como implementar função impulso, no R? Do tipo delta de Dirac: O sinal delta de Dirac δ(t) foi criado pelo físico inglês Paul A. Maurice Dirac (1902 - 1984) e pode ser visto como o equivalente contínuo do delta de Kronecker (matemático alemão,1823 - 1991) e como o limite de uma distribuição normal Grato, Marcelo.

1 0

Ajuda em uma função
by Sérgio Henrique almeida da silva ju 15 Jun '12

15 Jun '12

Olá Gostaria de fazer uma função e como não entendo nada, gostaria de pedir ajuda de vocês Eu uso o Alfa de Cronbach para analisar a consistência interna de uma escala, normalmente se calculo o alfa total e o alfa se uma variável for retirada. O alfa total é fácil de fazer, porém quando eu vou fazer o calculo de quando retiro uma variável tenho que fazer tudo manualmente. Ex.: a = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2) b = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1) c = c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5) d = c(2,2,3,3,4,3,2,3,4,5,5,5,4,3,2,3,4,5,4,3,2) dados = data.frame(a,b,c,d) library(psy) # Alfa Total cronbach(dados[,1:4]) alfa total = 0.5324005 Agora preciso retirar a variável "a" e rodar novamente o calculo, depois colocar ela denovo no banco e retirar a variável "b" e assim sucessivamente, dessa forma terei o valor do alfa caso eu retire uma variável. # Alfa retirando a variável a = 0.220339 b = 0.3239247 c = 0.7669716 d = 0.2198046 Qual a melhor forma de fazer isso através de uma função? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 5

Re: [R-br] Tornar mensagens de erro "invisíveis" em função
by Ivan Bezerra Allaman 15 Jun '12

15 Jun '12

Tente isto, d<-data.frame(x=1:10, y=rep(c(FALSE,TRUE),c(4,6))) d glm(y~x, data=d, family=binomial) keepWarnings <- function(expr) { localWarnings <- list() value <- withCallingHandlers(expr, warning = function(w) { localWarnings[[length(localWarnings)+1]] <<- w invokeRestart("muffleWarning") }) list(value=value, warnings=localWarnings) } z <- keepWarnings(glm(y~x, data=d, family=binomial)) z z$value (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}

3 2

R2 para modelos Poisson ou Binomial Negativo
by Sérgio Henrique almeida da silva ju 14 Jun '12

14 Jun '12

Olá Alguém sabe como estimar o R2 para modelos de Poisson ou Binomial negativo? Encontrei esse documento falando da possíbilidade dessa estimação: http://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0CFcQ… Porém não achei de como calcular no R. Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

7 7

Estacionariedade Regressão Espacial
by Isabela Zardo 14 Jun '12

14 Jun '12

Boa tarde, Alguém sabe como eu faço pra verificar a estacionariedade no R para a regressão espacial? Já que é preciso verificar para que o resulatdo do índice de Moran seja válido. Preciso de ajuda com urgência!!! Agradeço!! Isabela Zardo

1 0

Dúvida com quando o array tem NA
by Victor Eduardo 14 Jun '12

14 Jun '12

Galera, estava com um array de tamanho 150 onde tinham matrizes 1x1. Sendo que em alguns tinham NA. Dentro de minha função criei uma variável auxliar para identificar esses NA. Fiz do seguinte modo: aux<-NULL for(t in 1:n){ if(is.na(Y[,,t])==TRUE){ aux[t]<-0 }else{ aux[t]<-1 } } Ao identificar, o aux terá tamanho de 150, de 0 e 1. Feito essa verificação, estou pulando para a parte debaixo que é assim: for(t in 1:n){ if(aux[t]==1){ B<-matrix(Z[,,t],ncol=m) U[,,t]<-Y[,,t]-B%*%At[,,t]-D[,,t] F[,,t]<-B%*%Pt[,,t]%*%t(B)+H[,,t] K<-T[,,t]%*%Pt[,,t]%*%t(B)%*%ginv(F[,,t]) L[,,t]<-T[,,t]-K%*%B At[,,t+1]<-T[,,t]%*%At[,,t]+K%*%U[,,t]+C[,,t] Pt[,,t+1]<-T[,,t]%*%Pt[,,t]%*%t(L[,,t])+R[,,t]%*%Q[,,t]%*%t(R[,,t]) }else{ At[,,t+1]<-T[,,t]%*%At[,,t]+C[,,t] Pt[,,t+1]<-T[,,t]%*%Pt[,,t]%*%t(L[,,t])+R[,,t]%*%Q[,,t]%*%t(R[,,t]) } } return(list(At,Pt)) Só que o R está dando um erro que não estou conseguindo resolver, o erro é esse: Erro em svd(X) : infinite or missing values in 'x' Abraços!

1 1

Erro em plot()
by Ricardo Dobrovolski 13 Jun '12

13 Jun '12

Oi Pessoal, Eu quero plotar 18 gráficos de dispersão. No entanto, eu tive a seguinte mensagem de erro: "Erro em plot.new() : margens da figura muito grande" Os comandos que utilizei são os seguintes: par(mfrow=c(6,3)) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.na "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.na "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.na "]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.rn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.rn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.rn"]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.ra"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.ra"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.ra"]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) plot(agro.hdi[,"HDI_UN_2011"],agro.hdi[,"prop.agroprod.nn"]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"GDP_PerCapita_2010_WorldBank"],agro.gdp[,"prop.agroprod.gn "]*100,ylim=c(0,100),cex=1.4) plot(agro.gdp[,"Agro_GDP_2007_WorldBank"],agro.gdp[,"prop.agroprod.ga "]*100,ylim=c(0,100),cex=1.4) onde agro.gdp e agro.hdi são duas tabelas de dados e os termos entre aspas representam nomes das minhas colunas. Atenciosamente, RD.

3 2

Ordenar linhas em uma matriz
by Eric Ferreira 13 Jun '12

13 Jun '12

Prezados colegas, Estou com o seguinte problema: desejo ordenar uma matriz de acordo com os nomes de suas linhas. Seja o exemplo, M<-matrix(c(4,5,6,3,3),5,1) rownames(M)<-c('c','d','e','a','b') A matriz ordenada que eu desejo obter é: M1<-matrix(c(3,3,4,5,6),5,1) rownames(M1)<-c('a','b','c','d','e') Alguém conhece algum comando do tipo "sort by..."?? Obrigado e abraços, -- Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

3 3

Re: [R-br] R2 para modelos Poisson ou Binomial Negativo
by Ivan Bezerra Allaman 13 Jun '12

13 Jun '12

Sua resposta aqui, http://r.789695.n4.nabble.com/compute-coefficient-of-determination-R-square… (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

Tornar mensagens de erro "invisíveis" em função
by Emmanuel Arnhold 12 Jun '12

12 Jun '12

É possível tornar "invisíveis" as mensagens de erros de uma função? Incluindo mensagens de avisos perdidas. Emmanuel

1 0

Re: [R-br] Alteração em função de pacote (novo tópico)
by Junior Beleti 12 Jun '12

12 Jun '12

Okay Benilton, muito obrigado pela atenção. Sobre a questão do tempo ser linear, eu também acreditava nessa possibilidade, tanto é que em uma outra função do pacote geoComp, função "volta.quad" eu consegui tempo linear usando mclapply do multicore. A questão é a função "volta.cokri" mesmo, onde consegui tempos bons apenas para 2 núcleos. Aumentando a quantidade de núcleos não se tem um bom ganho. Eu acredito que seja devido a chamada que a "volta.cokri" faz dentro da "mvrnorm", onde faz o carregamento da "eigen.o" (eigen.f). Ou seja, ainda bato na tecla que o problema esta na chamada a função "eigen" fortran. -- Junior Beleti

2 1

Re: [R-br] Alteração em função de pacote (novo tópico)
by Junior Beleti 12 Jun '12

12 Jun '12

O "script" ao que eu me refiro nada mais é do que uma sequencia de comandos, que pode ser visto a seguir: require(geoComp) require(MASS) require(statmod) require(geoR) data(pivo) dados <- pivo[,c(6,7,8,1,2)] dados <- as.geoComp(dados) bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) estima <- mec(dados) gr <- pred_grid(bor, by=4) md.cov.ck <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados) preditos.gh <- volta.quad(md.cov.ck,n.pontos=7,Variancia=FALSE) preditos.gh <- data.frame(preditos.gh) write.table(preditos.gh,"pred_by4k7ns1000MBM.txt") system.time(preditos.simu <- volta.cokri(md.cov.ck ,num.simu=100,int.conf=0.95)) Informações sobre o computador e sessionInfo: Release: Ubuntu 11.10 (oneiric) Kernel: 3.0.0-19-generic (#33-Ubuntu SMP Thu Apr 19 19:05:14 UTC 2012) OS Type: Linux GCC version: 4.6.1 (x86_64-linux-gnu) Vendor: GenuineIntel CPUs: 8 Model name: Intel(R) Core(TM) i7 CPU Q 720 @ 1.60GHz Frequency: 933.000 MHz L2 cache: 6144 KB R version 2.15.0 (2012-03-30) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] geoR_1.7-2 sp_0.9-99 statmod_1.4.14 MASS_7.3-17 geoComp_1.0 loaded via a namespace (and not attached): [1] grid_2.15.0 lattice_0.20-6 RandomFields_2.0.54 [4] splancs_2.01-31 tools_2.15.0 Att, -- Junior Beleti

2 1

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 08 Jun '12

08 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 07 Jun '12

07 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

Duvida envelope simulado para ZINB
by Humberto Hazin 07 Jun '12

07 Jun '12

Olá Pessoal, Estou precisando fazer uma análiser de resíduo para um ZINB e gostaria de aplicar um envelope simulado! existe algum pacote ou codico que faça? Humberto

2 2

[OFF-TOPIC] Free Webinar: Spatial Analysis in R
by ASANTOS 07 Jun '12

07 Jun '12

Bom tarde Pessoal, Para quem trabalha com estatística espacial ou tem interesse sobre o assunto, encontrei um webinar grátis interessante na próxima semana (13/06 as 14:00h), intitulado: Finding Meaning in Points, Areas and Surfaces: Spatial Analysis in R com o professor David Unwin, segue o link para inscrição: http://www.revolutionanalytics.com/news-events/free-webinars/2012/spatial-a… Abraço, -- 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

1 0

Re: [R-br] Help me please! ^^
by Ivan Bezerra Allaman 07 Jun '12

07 Jun '12

Irei partir do princípios que as medidas repetidas são independentes. Se o interesse é comparar mais de duas médias podemos usar a técnica da análise de variância seguida de um teste de comparação múltipla. Segue o script. DIA1 <- c(7.306,7.296,7.326,7.334,7.334) DIA2 <- c(7.244,7.258,7.288,7.304,7.240) DIA3 <- c(6.998,7.044,7.044,7.076,7.086) DIA4 <- c(6.940,6.956,6.98,6.998,7.002) DIA5 <- c(7.116,7.130,7.142,7.148,7.150) pot <- c(DIA1,DIA2,DIA3,DIA4,DIA5) dias <- rep(c('DIA1','DIA2','DIA3','DIA4','DIA5'),rep(5,5)) dad <- data.frame(dias,pot) dad mod <- aov(pot ~ dias,data=dad) summary(mod) plot(mod)#Avaliando os pressupostos. library(agricolae) HSD.test(mod,'dias') (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

Help me please! ^^
by Patricia alcantara nguyen 07 Jun '12

07 Jun '12

Olá, Tudo bem? Preciso fazer um teste de repetitividade da potência de uma balança, mas estou tendo dificuldades com o R. Uma balança de radiação foi usada para medir a potência de um equipamento de ultrassom (3 MHz) durante 5 dias. Em cada dia foram feitas 5 medições de uma determinada potência. Como implementar isso no R??? Quero verificar se a média da potência dos 5 dias são iguais. Potência Nominal - 7 W MEDIDA 1 MEDIDA 2 MEDIDA 3 MEDIDA 4 MEDIDA 5 DIA 1 7,306 7,296 7,326 7,334 7,334 DIA 2 7,244 7,258 7,288 7,304 7,240 DIA 3 6,998 7,044 7,044 7,076 7,086 DIA 4 6,940 6,956 6,98 6,998 7,002 DIA 5 7,116 7,130 7,142 7,148 7,150 Muito Obrigada! Patrícia Nguyen

2 1

plotar valores específicos com função image( )
by Vitor Aguiar 06 Jun '12

06 Jun '12

Olá, desculpe por estar postando a pergunta novamente, mas não tive resposta da primeira vez. Se alguém puder ajudar, agradeço muito. Tenho uma matrix com que estou fazendo um heatmap usando a função image( ). Gostaria de plotar não apenas os valores observados na matrix, mas também valores que seriam possíveis mas não são observados. Isso significa o seguinte: Tenho uma matrix 16x16. Valores possíveis, mas não observados, correspondem a "0" (zeros) onde a soma (índice da linha + índice da coluna) são iguais ou menores que 15. Exemplo: Se eu tenho um valor "0" na posição [14, 3], eu NÃO quero plotar, pois 14 + 3 = 17. Se eu tenho um valor "0" na posição [11, 3], eu quero plotar, pois 11 + 3 = 14. Eu queria que esses valores fossem plotados em cor diferente do meu heatmap, por exemplo, em azul. É possível fazer isso? Matrix: x = c(3045, 893, 692, 830, 617, 155, 246, 657, 105, 60, 18, 7, 7, 4, 2, 11234, 2985, 2242, 2471, 1575, 366, 503, 1283, 170, 79, 32, 6, 4, 1, 3, 19475, 4756, 3233, 3251, 1810, 409, 575, 1210, 139, 41, 11, 4, 2, 0, 0, 20830, 4739, 2990, 2531, 1346, 298, 325, 612, 60, 17, 1, 0, 1, 0, 0, 15304, 3196, 1885, 1440, 610, 117, 115, 185, 14, 2, 0, 0, 0, 0, 0, 8026, 1535, 806, 539, 223, 33, 37, 39, 0, 0, 0, 0, 0, 0, 0, 3300, 562, 286, 141, 45, 14, 5, 12, 0, 0, 0, 0, 0, 0, 0, 1067, 160, 65, 40, 14, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 277, 47, 6, 2, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 72, 6, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) xmat = matrix(x, ncol = 12) xmat = cbind(xmat, rep(0, 15), rep(0, 15), rep(0, 15), rep(0, 15)) xmat = rbind(xmat, rep(0, 16)) dimnames(xmat) = list(0:15, 0:15) image(log(xmat)) Muito Obrigado, Vitor

2 1

Re: [R-br] Alteração em função de pacote (novo tópico)
by Junior Beleti 06 Jun '12

06 Jun '12

Olá Benilton.. Com as informações que você apresentou sobre as execuções, realizei a execução do script novamente. Passando como parâmetro primeiramente (num.simu=1000), a qual ficou executando por mais de 18 horas e não conseguiu chegar ao final. Na segunda execução, utilizei (num.simu=100) e obtive o seguinte resultado, quanto ao tempo: usuário sistema decorrido 60869.136 76.144 61005.444 A chamada a função foi a seguinte, modificando apenas o valor do parâmetro: system.time(preditos.simu <- volta.cokri(md.cov.ck ,num.simu=100,int.conf=0.95)) Gostaria de saber como você executou para ter esses tempos tão pequenos. Att, Junior Beleti.

2 1

Escrever diretorio formato windows em arquivo
by Leandro Marino 06 Jun '12

06 Jun '12

Caros, estou preparando uns arquivos para uma simulação é preciso que o R guarde um vetor da seguinte forma x <- 't:\pasta' Isto depois será usado por um outro programa. No entanto, recebo a mensagem de erro Erro: '\p' é uma seqüência de escape não reconhecida na cadeia de caracteres começando com "T:\p" Alguma sugestão? Abraços, Leandro

5 10

list para data.frame
by Leandro Marino 05 Jun '12

05 Jun '12

Caros, estou com um grande dúvida, tenho a seguinte lista x <- list() x$s1 <- list() x$s2 <- list() x$s1$tp1 <- data.frame(v1=1:10,v2=3:12,s=1,tp=1) x$s1$tp2 <- data.frame(v1=1:10,v2=3:12,s=1,tp=2) x$s1$tp3 <- data.frame(v1=1:10,v2=3:12,s=1,tp=3) x$s2$tp1 <- data.frame(v1=1:10,v2=3:12,s=2,tp=1) x$s2$tp2 <- data.frame(v1=1:10,v2=3:12,s=2,tp=2) x$s2$tp3 <- data.frame(v1=1:10,v2=3:12,s=2,tp=3) quero juntar de tal forma que : s1 <- rbind(x[[1]][[1]],x[[1]][[2]],x[[1]][[3]]) s2 <- rbind(x[[1]][[1]],x[[1]][[2]],x[[1]][[3]]) no entanto, s varia de 1 a 300. isto é possivel?

2 2

Ajuda em modelo não linear
by Leonardo Guimarães de Oliveira 05 Jun '12

05 Jun '12

Como faço para encontrar o ponto na curva obtida na CMR abaixo onde a curva muda de direção. x=structure(list(Vt = c(0, 0.34, 3.68, 10.38, 21.94, 32.08, 48.54, 79.69, 94.82, 107.99), t = c(0L, 3L, 6L, 9L, 12L, 16L, 24L, 36L, 48L, 72L)), .Names = c("Vt", "t"), class = "data.frame", row.names = c(NA, -10L)) m4=nls(Vt~(107.99)/(1+exp(2-4*k*(t-g))), data=x, start=c(k=0.05,g=0.1)) m4 f4=function(x){(107.99)/(1+exp(2-4*0.02853*(x-8.71643)))} curve(f4, 1,100) Leonardo G. Oliveira Médico Veterinário Doutorando em Ciência Animal - EVZ/UFG

2 1

pergunta simples - diagonal de matrix
by Vitor Aguiar 04 Jun '12

04 Jun '12

Boa noite, se eu tenho uma matrix, por exemplo, 7x7 eu posso facilmente extrair a diagonal usando diag(mat). Mas e se eu quiser algo como "uma segunda diagonal", uma diagonal composta de mat[1,2] mat[2, 3] mat[3, 4] mat[4, 5] mat[5, 6] mat[6, 7]? Tem alguma forma simples de fazer? obrigado, Vitor

2 1

Re: [R-br] Alteração em função de pacote (Junior Beleti)
by Junior Beleti 04 Jun '12

04 Jun '12

Olá Benilton, meus estudos se iniciaram após uma conversa com a professora Ana Beatriz, visto que ela tinha problemas de desempenho no trabalho desenvolvido por ela e pelo professor Paulo. Eu fui analisando os pontos onde ela me disse que teria problemas, e identifiquei na volta.cokri um maior tempo de processamento. Vou processar novamente os scripts, mas em minha máquina o tempo de processamento fui bem maior. Volto a postar os resultados. E o professor Paulo está certo pois esse trabalho não visava desempenho, e como ele mesmo citou, vi uma oportunidade para utilizar processamento paralelo e tentar ajudar no desempenho dessa aplicação. Muito obrigado pela ajuda. Junior Beleti.

1 0

borda externa de mapa
by Renato Braga 04 Jun '12

04 Jun '12

Lista, como eu faco para plotar somente a borda externa de um mapa no formato shapefile? Grato, Renato

2 2

Python para R
by Daniel Marcelino 03 Jun '12

03 Jun '12

Caros, alguem da lista conhece programação python? Estou há dias tentando desenvolver uma rotina em R para explorar dados eleitorais, mas não tenho conseguido ir muito longe. Eu tenho um modelo do que mais ou menos eu preciso, mas a rotina está em python. Alguem encara me ajudar nisso? Daniel ----------------------------------------------------- My contacts: Phone 1+514 343 6111 Ext 3799 Skype: d.marcelino 3200 Jean Brillant, Office C5075 Montreal, QC; H3T 1N8 Canada -----------------------------------------------------

3 2

Combinar vetores com tamanhos diferentes...
by andrebvs 03 Jun '12

03 Jun '12

7 13

borda externa
by André Figueiras Rutz 02 Jun '12

02 Jun '12

Caros, como é que faz para inserir num gráfico a moldura externa (só a linha preta em volta do gráfico?) [image: Inline image 3] grato!

3 8

Gráfico
by Paulo Ricardo Silva Moreira 01 Jun '12

01 Jun '12

Boa tarde prezados, Estou precisando fazer um gráfico, plotando duas variáveis com suas respectivas escalas, uma de cada lado do gráfico. Alguém sabe me dizer qual comando eu poderia usar!? Desde já agradeço, Paulo Ricardo.

4 5

Alteração em função de pacote (novo tópico)
by Junior Beleti 01 Jun '12

01 Jun '12

Olá Benilton, tentando explicar melhor meu problema: Faço uso de um pacote chamada geoComp. Tal pacote é utilizado de forma geral para modelar o padrão espacial de dados composicionais. Foi um trabalho da professora Ana Beatriz. O problema está quanto ao tempo de processamento para realizar o trabalho necessário. Analisando mais a fundo, verifiquei que o maior problema quanto a desempenho está na função "volta.cokri" desse pacote. Nessa função existe uma chamada à função "mvrnorm" e é justamente nessa chamada que ocorre o maior problema de desempenho. Verifiquei também que na função "mvrnorm" existe uma chamada para a tal função "eigen", que se trata de um código-objeto fortran. Inicialmente utilizei a função "mclapply" para tentar melhorar tal desempenho, mas não obtive muito sucesso. O que estava pensando quando insisti na ideia de criar a nova função no pacote base, foi o de que se chamando funções "diferentes", ou seja, criando funções idênticas cada uma com o espaço de endereçamento e nomes diferentes, poderia aumentar o desempenho. A seguir, está um script utilizado: require(geoComp) require(MASS) require(statmod) require(geoR) data(pivo) dados <- pivo[,c(6,7,8,1,2)] dados <- as.geoComp(dados) bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) ## bor é uma matriz 2x102 estima <- mec(dados) gr <- pred_grid(bor, by=4) md.cov.ck <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados) preditos.gh <- volta.quad(md.cov.ck,n.pontos=7,Variancia=FALSE) preditos.gh <- data.frame(preditos.gh) write.table(preditos.gh,"pred_by4k7ns1000MBM.txt") preditos.simu <- volta.cokri(md.cov.ck,num.simu=1000,int.conf=0.95) preditos.simu <- data.frame(preditos.simu[[1]]) preditos.simu.ic <- data.frame(preditos.simu[[2]]) write.table(preditos.simu,"predsimu_by4k7ns1000MBM.txt") write.table(preditos.simu.ic,"predsimuic_by4k7ns1000MBM.txt") O pacote geoComp pode ser adquirido no endereço: ( http://www.din.uem.br/~pg45178/geocomp/geoComp.tar.gz) Um resumo do trabalho pode ser encontrado em: ( http://www.leg.ufpr.br/lib/exe/fetch.php/pessoais:abtmartins:rbb09_2.pdf) Obrigado pela atenção. Carlos.

3 3

Re: [R-br] Alteração em função de pacote
by Junior Beleti 01 Jun '12

01 Jun '12

Bom dia Fernando Mayer, segui suas orientações, mas obtive o seguinte erro ao tentar instalar o pacote (R CMD INSTALL base_2.15.0.tar.gz). * installing to library ‘/usr/local/lib/R/library’ * installing *source* package ‘base’ ... ** R ** demo ** inst ** preparing package for lazy loading Erro em asNamespace(ns, base.OK = FALSE) : operation not allowed on base namespace * removing ‘/usr/local/lib/R/library/base’ * restoring previous ‘/usr/local/lib/R/library/base’ Obrigado pela atenção e se tiver alguma outra sugestão. Att, -- Carlos

4 3

Re: [R-br] Alteração em função de pacote
by Junior Beleti 31 Mai '12

31 Mai '12

Olá Leonard. Estou criando uma função no pacote base, função esta semelhante à função eigen.f já existente. A nova função teria o nome eigem.f Entendi Benilton. Vou tentar colocar um código reproduzivel aqui. Obrigado, Carlos.

1 0

Criar um loop
by Humberto Hazin 31 Mai '12

31 Mai '12

Caros, Venho mais uma vez solicitar ajuda a vocês tenho um banco de dados com cerca de 80 mil linhas conforme o exemplo abaixo. gostaria de atribuir um id para cada barco, exemplo toda a coluna que tivesse o Barco Vieirasa a coluna ID = 1 e assim por diante. Alguem pode me ajudar? Humberto boat Id VIEIRASA COSTA.DE.SAO.MIGUEL VIEIRASA GUARISTE.PRIMEIRO COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS GAVILAN.DEL.MAR VIEIRASA GUARISTE GUARISTE RAYMI MAR.DE.LOS.SARGAZOS RAYMI GUARISTE RAYMI COSMOS COSMOS COSMOS

2 2

função image( ) - plotar valores selecionados
by Vitor Aguiar 31 Mai '12

31 Mai '12

Boa tarde, alguém poderia me ajudar com a função image( )? Me desculpe por abusar do tempo de vocês, o texto abaixo é grande porém simples. Tenho uma matrix de valores com que estou fazendo um heatmap usando a função image( ). Gostaria de plotar não apenas os valores observados na matrix, mas também valores que seriam possíveis mas não são observados. Isso significa o seguinte: Tenho uma matrix 16x16. Valores possíveis, mas não observados, correspondem a "0" (zeros) onde a soma (índice da linha + índice da coluna) são iguais ou menores que 15. Exemplo: Se eu tenho um valor "0" na posição [14, 3], eu NÃO quero plotar, pois 14 + 3 = 17. Se eu tenho um valor "0" na posição [11, 3], eu quero plotar, pois 11 + 3 = 14. Eu queria que esses valores fossem plotados em cor diferente do meu heatmap, por exemplo, em azul. Matrix: x = c(3045, 893, 692, 830, 617, 155, 246, 657, 105, 60, 18, 7, 7, 4, 2, 11234, 2985, 2242, 2471, 1575, 366, 503, 1283, 170, 79, 32, 6, 4, 1, 3, 19475, 4756, 3233, 3251, 1810, 409, 575, 1210, 139, 41, 11, 4, 2, 0, 0, 20830, 4739, 2990, 2531, 1346, 298, 325, 612, 60, 17, 1, 0, 1, 0, 0, 15304, 3196, 1885, 1440, 610, 117, 115, 185, 14, 2, 0, 0, 0, 0, 0, 8026, 1535, 806, 539, 223, 33, 37, 39, 0, 0, 0, 0, 0, 0, 0, 3300, 562, 286, 141, 45, 14, 5, 12, 0, 0, 0, 0, 0, 0, 0, 1067, 160, 65, 40, 14, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 277, 47, 6, 2, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 72, 6, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) xmat = matrix(x, ncol = 12) row15 = rep(0,16) col12 = rep(0,15) col13 = rep(0,15) col14 = rep(0,15) col15 = rep(0,15) xmat = cbind(xmat, col12, col13, col14, col15) xmat = rbind(xmat, row15) dimnames(xmat) = list(0:15, 0:15) xmat # minha matrix original é a primeira xmat (15x12), mas estou criando essa outra com linhas e colunas de zeros pois, por exemplo, 0 + 15 = 15, e preciso plotar esse valor também. Mas se alguém tiver solução mais simples... Desde já agradeço, Vitor Vitor Rezende da Costa Aguiar --------------------------------------------------- PhD candidate in biotechnology Rede Nordeste de Biotecnologia UFES, Brazil --------------------------------------------------- Current Adress: Department of Integrative Biology University of California, Berkeley 2033 Valley Life Sciences Building #4134 Berkeley, CA 94720 Phone: 1 (510) 473-2070

1 1

Re: [R-br] Alteração em função de pacote
by Junior Beleti 31 Mai '12

31 Mai '12

Entendido Benilton. Mudando o foco, você pode me sugerir alguma biblioteca paralela para trabalhar com a função eigen? Fiz experimentos com a multicore, mas não obtive grandes resultados. Att, Carlos

6 7

rcorr do Hmisc
by Edson Lira 31 Mai '12

31 Mai '12

Gente, quero fazer uma correlacção entre n variáveis, no Hmisc tem a função rcorr,vejam o exemplo abaixo: require(HmisC) x <- c(-2, -1, 0, 1, 2) y <- c(4, 1, 0, 1, 4) z <- c(1, 2, 3, 4, NA) v <- c(1, 2, 3, 4, 5) rcorr(cbind(x,y,z,v)) No meu banco são 90 variáveis como fzer o teste acima sem ter que declarar uma por uma? Algo do tipo rcorr(cbind(x:v)) Bom dia! Edson Lira Estatístico Manaus-Amazonas

2 1

Retirar dados de uma função no R
by Ana Lúcia Souza Silva Mateus 29 Mai '12

29 Mai '12

Saudações, Tenho a seguinte dúvida, se algúem poder me responder ficaria grata. Estou trabalhando com a seguinte função no R (stdiagn), apresentada abaixo: *> stdiagn function (pts, stkh, stse, stmc = 0, Dzero = FALSE) { oldpar <- par(mfrow = c(2, 2)) pointmap(as.points(pts)) title("Data map") oprod <- outer(stkh$ks, stkh$kt) st.D <- stkh$kst - oprod st.R <- st.D/(stse) if (!Dzero) { persp(stkh$s, stkh$t, st.D, xlab = "Distance", ylab = "Time", zlab = "D", expand = 0.5, ticktype = "detailed", theta = -30, shade = 0.4, cex = 0.7) title("D plot") } else { persp(stkh$s, stkh$t, st.D/oprod, xlab = "Distance", ylab = "Time", zlab = "Dzero", expand = 0.5, ticktype = "detailed", theta = -30, shade = 0.4, cex = 0.7) title("Dzero plot") } plot(outer(stkh$ks, stkh$kt), st.R, xlab = "K(s)K(t)", ylab = "R") abline(h = 0) title("Residual Plot") if (length(stmc) > 1) { hist(c(stmc$t0, stmc$t), nclass = 50, include.lowest = TRUE, xlab = "Test statistic", main = "") y.75 <- par()$usr[4] * 0.75 y.8 <- par()$usr[4] * 0.8 lines(rep(stmc$t0, 2), c(0, y.75), lwd = 2) text(stmc$t0, y.8, "Data Statistic") title("MC results") } else { plot(c(0, 1), c(0, 1), type = "n", axes = FALSE, xlab = "", ylab = "") ad <- par()$adj par(adj = 0.5) text(0.5, 0.5, "No MC available") par(adj = ad) } par(oldpar) } <environment: namespace:splancs> > * Gostaria de obter os valores de * st.D e st.R . Desde já agradeço pelo auxilio. *

5 8

R com php
by Vinicius Brito Rocha 29 Mai '12

29 Mai '12

Senhores, Estou precisando usar o R como uma engine por trás de um portal em desenvolvido em php. Estou precisando de orientação/material. Alguém já fez algum tipo de integração R + web? abs -- 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.

2 1

frequência de vetores
by Márcio Augusto Ferreira Rodrigues 29 Mai '12

29 Mai '12

Prezados, tenho um arquivo da seguinte forma z10z11z12z13z20z21z22z23z30z31z32z33 000100100010 010000010010 010000010001 001000100001 010000100100 010000010010 010001000100 cada linha eu considero como um vetor e gostaria de saber a frequência de cada um, pois tenho 1000 linhas e sei que so possui 64 tipos diferentes de vetores. Assim gostaria de saber a frenquência de cada um dos 64 vetores diferentes que tenho . Grato, Márcio

5 6

Survey
by Emerson Cotta Bodevan 28 Mai '12

28 Mai '12

Prezados, preciso criar uma máscara de um questionário para aplicá-lo numa comunidade. Aquela máscara que o Epi-info nos permite criar. O R tem recurso para isso? Senão, alguém teria indicação de um outro software que não seja o Epi-info? É que gostaria de algum que fosse multiplataforma ou que rodasse pelo menos no LINUX Ubuntu. Agradeço qualquer ajuda, -- *Emerson*

6 6

Fwd: Ajuda com vetores
by Fernando Colugnati 28 Mai '12

28 Mai '12

Prezados, tenho a seguinte matriz de estrutura Id a1 p1 p2 p3 p4 p5 p6 p7 1 2 3 4 0 0 0 0 0 1 6 3 0 0 0 0 0 0 2 7 5 3 2 0 0 0 0 2 3 6 0 0 0 0 0 0 Preciso de uma matriz a partir desta da seguinte forma: Id a1 p 1 2 3 1 2 4 1 6 3 2 7 5 2 7 3 2 7 2 2 3 6 Ou seja, preciso gerar quantas linhas forem diferentes de 0 os p1 - p7. É possível fazendo via um apply, sem loop? Abs e obrigado -- Fernando A.B. Colugnati -- Fernando A.B. Colugnati

2 2