Janeiro 2013 - R-br

GLMM com dados categóricos
by Pavel Dodonov 24 Jan '13

24 Jan '13

Saudações, É a primeira vez que escrevo para cá, perdoem-me caso eu quebre alguma norma... Eu preciso ver os efeitos de algumas variáveis (contínuas e categóricas) sobre uma variável categórica com quatro níveis. Além disso, os dados estão aninhados em dois grupos, então acho que precisa ser um modelo misto. Explicando melhor: É um experimento com predação de ninho. A variável-resposta é o tipo de predador ("Ave", "Mamífero", "Desconhecido", "Nenhum"), as variáveis explanatórias são: Local (dois níveis - "Borda" e "Interior"), altura (variável contínua), e densidade (variável contínua, limitada entre 0 e 1 mas que na prática não ultrapassa 0.3). Tem também o fator aleatório "Campanha" (dois níveis, "Outubro" e "Novembro"), que até poderia entrar como variável explanatória, mas que seria melhor usada como fator aleatório porque não estamos interessados na diferença entre as campanhas. Eu pensei em fazer quatro análises usando cada nível da variável-resposta como uma variável binária e fazer um GLMM com distribuição binomial, MAS seria melhor se desse pra considerar os quatro níveis de uma vez. Talvez um GLMM com distribuição multinomial? Se alguém tiver alguma sugestão, ficarei imensamente agradecido. Att - Pavel ------- Pavel Dodonov Biologist, PhD student in Ecology and Natural Resources, São Carlos Federal University (UFSCar), SP, Brazil http://ufscar-br.academia.edu/pdodonov *"The highest function of ecology is understanding consequences." - Frank Herbert, Dune, 1965*

1 0

Mensagem de erro: cannot allocate memory block of size 3.2
by Cláudio Brisolara 23 Jan '13

23 Jan '13

Caros colegas Estou precisando de ajuda para identificar e solucionar um problema relativo a alocação de memória no R. Muitas mensagens já foram postadas sobre esse assunto, mas não consegui localizar uma solução. Por essa razão, solicito a ajuda de alguém. Em algumas operações que faço recebo a seguinte mensagem: “Error: cannot allocate memory block of size 3.2 GbIn addition: Warning message:restarting interrupted promise evaluation” A rotina que estou tentando executar é a seguinte: mu1<-c(8.035153, 3.931724) # fixa vetor de médias sigma1<-matrix(c(0.014179, 0.003042, 0.003042, 0.026259), 2, 2) # fixa matriz de variâncias set.seed(23456) n<-20000 dados.sim1<-rmvnorm(n, mu1, sigma1, method=c("chol")) # gera vetor normais correlacionados Ocorre que a questão da memória também é reportado em outras situações como por exemplo na simples chamada de um pacote: require(tseries) Loading required package: tseries Error in get(Info[i, 1], envir = env) : cannot allocate memory block of size 2.7 Gb Eu usei os comandos memory.limit() and memory.size() para checar a limitação de memória, mas não observei nenhum problema. Encaminho também as informações obtidas pelo sessionInfo(). > memory.limit() [1] 6004 > memory.size() [1] 1361.88 > sessionInfo() R version 2.15.2 (2012-10-26) Platform: x86_64-w64-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 loaded via a namespace (and not attached): [1] grid_2.15.2 quadprog_1.5-4 stabledist_0.6-5 tools_2.15.2 [5] xtable_1.7-0 O mais intrigante é que o problema realmente não é de memória, pois executei totalmente o mesmo script em um computador bem mais velho, com menos memória. É um problema de configuração da máquina, arquivo perdido algo do gênero? Já reinstalei o R e o RStudio, mas não obtive êxito. Por favor, algum colega poderia me ajudar a entender o que está acontecendo? Antecipadamente agradeço. Atenciosamente, Cláudio Brisolara

1 0

Re: [R-br] Legenda no eixo X entre dois gráficos
by Humberto Hazin 23 Jan '13

23 Jan '13

Tito muito obrigado consegui!!! André segue o código correto e com as tuas perguntas (que antes não havia entendido) respondidas Abr Humberto pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis <- c(length(tableMat1a$or):1) plot(tableMat1a$or,y.axis, type = "p", axes = T, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) #axis(1,at=seq(0,1.2,by=0.1),labels = T,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis <- c(length(tableMat1b$or):1) plot(tableMat1b$or,y.axis, type = "p", axes = T, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) #axis(1,at=seq(0,1.2,by=0.1),labels = T,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") mtext("Relative Catchability",side=1,line=3,adj =-0.6) dev.off() Enviado do Email do Windows De: Humberto Hazin Enviado: ‎23‎ de ‎janeiro‎ de ‎2013 ‎16‎:‎13 Para: R_br, R_br Assunto: Re: [R-br] Legenda no eixo X entre dois gráficos Olá André e Tito, Eu realmente não entendi tua pergunta André! Eu vou tentar aqui a solução pela dica do Tito e se conseguir posto aqui o resultado abr Humberto Enviado do Email do Windows De: andrebvs(a)bol.com.br Enviado: ‎23‎ de ‎janeiro‎ de ‎2013 ‎13‎:‎46 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Legenda no eixo X entre dois gráficos O seu gráfico ainda não está igual ao que deseja, falta colocar os valores no eixo X e fechar a caixa (box) onde ficam os intervalos. Como vc vai fazer isso? Em 23/01/2013 09:10, Tito Conte < tito.conte(a)gmail.com > escreveu: De um ?mtext Lá vai ter tudo o que você precisa! Tito Conte On Thu, Jan 17, 2013 at 12:37 PM, Humberto Hazin <hghazin(a)hotmail.com> wrote: Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo (https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Species lci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3 Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5 Sailfish 0.1257 0.25 0.498 6 White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci or uci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3 Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5 Sailfish 0.16000 0.20 0.4000 6 White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a<-read.csv("tableMetaGW.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill =TRUE) tableMat1b<-read.csv("tableMetaGB.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis plot(tableMat1a$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis plot(tableMat1b$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") dev.off() Enviado do Email do Windows _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Re: [R-br] Legenda no eixo X entre dois gráficos
by Humberto Hazin 23 Jan '13

23 Jan '13

Olá André e Tito, Eu realmente não entendi tua pergunta André! Eu vou tentar aqui a solução pela dica do Tito e se conseguir posto aqui o resultado abr Humberto Enviado do Email do Windows De: andrebvs(a)bol.com.br Enviado: ‎23‎ de ‎janeiro‎ de ‎2013 ‎13‎:‎46 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Legenda no eixo X entre dois gráficos O seu gráfico ainda não está igual ao que deseja, falta colocar os valores no eixo X e fechar a caixa (box) onde ficam os intervalos. Como vc vai fazer isso? Em 23/01/2013 09:10, Tito Conte < tito.conte(a)gmail.com > escreveu: De um ?mtext Lá vai ter tudo o que você precisa! Tito Conte On Thu, Jan 17, 2013 at 12:37 PM, Humberto Hazin <hghazin(a)hotmail.com> wrote: Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo (https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Species lci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3 Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5 Sailfish 0.1257 0.25 0.498 6 White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci or uci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3 Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5 Sailfish 0.16000 0.20 0.4000 6 White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a<-read.csv("tableMetaGW.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill =TRUE) tableMat1b<-read.csv("tableMetaGB.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis plot(tableMat1a$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis plot(tableMat1b$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") dev.off() Enviado do Email do Windows _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Legenda no eixo X entre dois gráficos
by Humberto Hazin 23 Jan '13

23 Jan '13

Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo (https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Species lci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3 Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5 Sailfish 0.1257 0.25 0.498 6 White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci or uci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3 Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5 Sailfish 0.16000 0.20 0.4000 6 White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a<-read.csv("tableMetaGW.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill =TRUE) tableMat1b<-read.csv("tableMetaGB.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis <- c(length(tableMat1a$or):1) plot(tableMat1a$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis <- c(length(tableMat1b$or):1) plot(tableMat1b$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") dev.off() Enviado do Email do Windows

3 2

vaga de emprego
by Manoel Galdino 23 Jan '13

23 Jan '13

Olá, a empresa onde trabalho tem uma vaga em aberto. Como nós usamos bastante o R aqui na empresa, pensei que essa lista poderia ser um lugar para achar um bom candidato. Segue a descrição da vaga. O trabalho é em São Paulo. Descrição do trabalho: A área de Analytics da Hyperativa <http://hyperativa.com.br/>trabalha com base de dados de clientes para desenvolver análises de perfil da base, modelos preditivos, plataformas de CRM e ajudar na tomada de decisão estratégica. As bases de dados são gerenciadas em SQL pela equipe de Banco de dados e o analista trabalha em conjunto com a equipe de banco de dados para extrair os dados e em seguida analisá-los. Análises típicas incluem análises descritivas, regressões linerares e modelos Bayesianos. Boa parte do tempo do candidato será gasto em operações de limpeza, agregações e transformações dos dados, bem como transferência de dados entre bases e/ou softwares (por exemplo, do SQL para R e do R para Excel). É comum que nós tenhamos que resolver problemas que não têm soluções prontas na empresa e é necessário aprender rapidamente novas ferramentas e/ou técnicas, e a verdade é que nós adoramos aprender novas técnicas e gostamos de desafios. Para nós é importante que, além de ser comprometido com o trabalho e com as entregas, você tenha uma vida com interesses além das típicas habilidades vendáveis no mercado. Gostamos de pessoas com bom humor, que curtam cachorro, joguem xadrez, andem de bicicleta, vão ao cinema, torçam pela desgraça do time rival e adorem tecnologia. Requisitos: - Superior completo – Preferencialmente com formação forte em métodos quantitativos - Inglês intermediário (Desejável leitura avançada) - Conhecimentos de Estatística (aplicada e teoria) - Office - Desejável Conhecimento de R - Conhecimentos de SQL – diferencial - Conhecimentos de programação em geral – Diferencial -- Manoel Galdino https://sites.google.com/site/galdinomcz/

4 5

Extrair diagonais de Matriz ao selecionar um elemento.
by Augusto Ribas 23 Jan '13

23 Jan '13

Bom dia a todos, estou com uma duvida aqui. Eu tenho varias matrizes dessa forma. > exemplo A B C D 1 0 1 0 0 2 0 0 2 0 3 1 0 0 3 4 0 2 0 0 Aqui ela esta 4x4, mas ela pode ser de qualquer forma, não necessariamente quadrada assim. O maior numero que eu tenho ali dentro é 3. Existe uma forma simples de obter a posição do 3? Eu estava tentando o seguinte: Obtinha o maior elemento max(exemplo) Depois olhava onde ele estava. > which(exemplo==3) [1] 15 Mas ele fala um número apenas, que ai tem que dividir o número de colunas, tem alguma função que me retornaria a posição, qual linha e qual coluna ele esta? Tipo [3,4], e que seja independente se a matrix for diferente quanto ao numero de linhas ou colunas? A segunda duvida é, tem alguma forma simples de extrair a diagonal dos números a partir do 3? Principalmente as posições dos elementos. No exemplo eu gostaria do seguinte: A diagonal do 3 assim: [1,2],[2,3],[3,4]. Que ai eu poderia obter quais as colunas que estão nessa diagonal, tipo B,C,D, e as linhas 1,2,3. Esse é o objetivo final, mas eu estou apanhando para obter um jeito simples de extrair a posição do ultimo elemento da diagonal ali na matriz, ou a posição dos elementos da diagonal em si. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 3

Comando que mantém a ordem natural dos dados
by Gisele Correia 23 Jan '13

23 Jan '13

Boa tarde, Gostaria de saber como plotar dados mantendo a ordem das variáveis. O R está ordenando automaticamente. exemplo: x <- c(2,1,4,2) plot(x,y) (aqui o x é mostrado na ordem 1,2,3 e 4 no gráfico). Desde já agradeço, obrigada!! -- Gisele Correia Engenheira de Pesca PPG-CIPET/UFAM "Agradeço por todas as dificuldades que enfrentei (e ainda enfrento)... se não fosse por elas, não teria saído do lugar. A facilidade nos impede de caminhar. Mesmo as críticas nos auxiliam muito..." Chico Xavier

4 4

Problema com TL 50 em análise de sobrevivência
by ASANTOS 22 Jan '13

22 Jan '13

Boa tarde pessoal, Estou com um problema para calcular a TL50 (tempo de morte de 50% do indivíduos de uma população), usando o modelo de sobrevivência paramétrico de Weibull, pois o valor extraído do summary do modelo com exp(Intercept) não esta batendo com o curve do modelo, para explicar melhor segue CRM: require(survival) ##Dados obtidos na web A <- read.table("http://www.mas.ncl.ac.uk/~nmf16/teaching/mas3311/A.dat") attach(A) ## Ajuste do modelo mod = survreg(Surv(time, status)~1,dist="weibull") summary(mod) ########################### Calculando os parâmetros alpha e tempo médio ########## Valor de Alpha (1/Scale) 1/0.891 # =1.122334 ##TL 50 exp do (Intercept) do modelo exp(3.887) #= 48.76437 ##Plotando par(cex.lab=1.8,bty="n",yaxs="r",xaxs="r",lwd=2,mar=c(4.5,4.5,1,1),mex=1.3,family="") curve(exp((-(48.76437)^(-1.122334))*(x^1.122334)),xaxt="n",yaxt="n",from=0,to=170,ylim=c(0,1),xlab="",ylab="") x=seq(0,170,10) y=seq(0,1,0.1) axis(1,x,lwd=2,cex.axis=1.5) axis(2,y,lwd=2,cex.axis=1.8) # abline(h=0.5, col="red") abline(v=48.76437, col="red") ## O que não bate o valor encontrado com a curva Alguém poderia de dar uma ajuda, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

1 0

Dúvida sobre intervalos....
by andrebvs＠bol.com.br 22 Jan '13

22 Jan '13

4 4