- R-br - DINF UFPR

Mandar email pelo R
by Pedro Emmanuel Alvarenga Americano do Brasil 15 Mar '13

15 Mar '13

Amigos de R, Sei que essa questão volta algumas vezes na lista mas estou empacado. Estou fazendo um data management de um projeto em que consiste em repetidas vezes, tipo toda semana, verificar quantos voluntários entraram no proejto, dados inconsistentes, dados ausentes e coisas assim. Então, eu consegui fazer um script para esse fim em que gera um arquivo com algumas tabelas e graficos. Agora eu gostaria de incluir nesse script o envio de um email pra mim mesmo com esse arquivo, para que eu possa verificar e depois encaminhar para os outros participantes. Sei que existe o sendmailR, mas não estou conseguindo fazer funcionar. Encontrei um catatau de gente referindo o mesmo erro em diversos diferentes foruns: > sendmail(from, to, subject, body,control=list(smtpServer=" ASPMX.L.GOOGLE.COM")) Erro em socketConnection(host = server, port = port, blocking = TRUE) : não é possível abrir a conexão Além disso: Mensagens de aviso perdidas: In socketConnection(host = server, port = port, blocking = TRUE) : ASPMX.L.GOOGLE.COM:25 cannot be opened Depois, consegui encontrar essa função em um dos foruns. O autor fala que funciona, mas comigo tem dado um erro de conectividade e eu não sabendo como dar a volta nesse erro. Será que tem alguma relação com o proxy? > send.email function(to, from, subject, message, attachment=NULL, username, password, server="smtp.gmail.com:587", confirmBeforeSend=TRUE){ # to: a list object of length 1. Using list("Recipient" = " recip(a)somewhere.net") will send the message to the address but # the name will appear instead of the address. # from: a list object of length 1. Same behavior as 'to' # subject: Character(1) giving the subject line. # message: Character(1) giving the body of the message # attachment: Character(1) giving the location of the attachment # username: character(1) giving the username. If missing and you are using Windows, R will prompt you for the username. # password: character(1) giving the password. If missing and you are using Windows, R will prompt you for the password. # server: character(1) giving the smtp server. # confirmBeforeSend: Logical. If True, a dialog box appears seeking confirmation before sending the e-mail. This is to # prevent me to send multiple updates to a collaborator while I am working interactively. if (!is.list(to) | !is.list(from)) stop("'to' and 'from' must be lists") if (length(from) > 1) stop("'from' must have length 1") if (length(to) > 1) stop("'send.email' currently only supports one recipient e-mail address") if (length(attachment) > 1) stop("'send.email' can currently send only one attachment") if (length(message) > 1){ stop("'message' must be of length 1") message <- paste(message, collapse="\\n\\n") } if (is.null(names(to))) names(to) <- to if (is.null(names(from))) names(from) <- from if (!is.null(attachment)) if (!file.exists(attachment)) stop(paste("'", attachment, "' does not exist!", sep="")) if (missing(username)) username <- winDialogString("Please enter your e-mail username", "") if (missing(password)) password <- winDialogString("Please enter your e-mail password", "") require(rJython) rJython <- rJython() rJython$exec("import smtplib") rJython$exec("import os") rJython$exec("from email.MIMEMultipart import MIMEMultipart") rJython$exec("from email.MIMEBase import MIMEBase") rJython$exec("from email.MIMEText import MIMEText") rJython$exec("from email.Utils import COMMASPACE, formatdate") rJython$exec("from email import Encoders") rJython$exec("import email.utils") mail<-c( #Email settings paste("fromaddr = '", from, "'", sep=""), paste("toaddrs = '", to, "'", sep=""), "msg = MIMEMultipart()", paste("msg.attach(MIMEText('", message, "'))", sep=""), paste("msg['From'] = email.utils.formataddr(('", names(from), "', fromaddr))", sep=""), paste("msg['To'] = email.utils.formataddr(('", names(to), "', toaddrs))", sep=""), paste("msg['Subject'] = '", subject, "'", sep="")) if (!is.null(attachment)){ mail <- c(mail, paste("f = '", attachment, "'", sep=""), "part=MIMEBase('application', 'octet-stream')", "part.set_payload(open(f, 'rb').read())", "Encoders.encode_base64(part)", "part.add_header('Content-Disposition', 'attachment; filename=\"%s\"' % os.path.basename(f))", "msg.attach(part)") } #SMTP server credentials mail <- c(mail, paste("username = '", username, "'", sep=""), paste("password = '", password, "'", sep=""), #Set SMTP server and send email, e.g., google mail SMTP server paste("server = smtplib.SMTP('", server, "')", sep=""), "server.ehlo()", "server.starttls()", "server.ehlo()", "server.login(username,password)", "server.sendmail(fromaddr, toaddrs, msg.as_string())", "server.quit()") message.details <- paste("To: ", names(to), " (", unlist(to), ")", "\n", "From: ", names(from), " (", unlist(from), ")", "\n", "Using server: ", server, "\n", "Subject: ", subject, "\n", "With Attachments: ", attachment, "\n", "And the message:\n", message, "\n", sep="") if (confirmBeforeSend) SEND <- winDialog("yesnocancel", paste("Are you sure you want to send this e-mail to ", unlist(to), "?", sep="")) else SEND <- "YES" if (SEND %in% "YES"){ jython.exec(rJython,mail) cat(message.details) } else cat("E-mail Delivery was Canceled by the User") } from <- list("Recipient" = "emmanuel.brasil(a)gmail.com") to <- list("Recipient" = "emmanuel.brasil(a)gmail.com") subject <- paste0("Pendencias Nomograma de ",Sys.time()) message <- "Mensagem automática com relatório de pendências do projeto Nomograma." attachment <- 'c:/banco/nomograma/Controle_recrutamento.doc' > send.email(to,from,subject,message,attachment,'emmanuel.brasil(a)gmail.com ','...',confirmBeforeSend=FALSE) Erro em jython.exec(rJython, mail) : (62, 'Connection refused') Qualquer luz é bem vinda. Abraço forte, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

3 7

Series temporais
by edmar Caldas 14 Mar '13

14 Mar '13

Estou rodando na versão 2.11.1 e quando peço a curva normal no histograma dá uma mensagem de erro (plot.new ainda não foi chamado). Mas na versão 2.8 funciona. Estranho. Pode me ajudar por favor. Segue o código em TXT.

1 0

"Ordenação Mergesort"
by Jackeline Bonetti Campos 14 Mar '13

14 Mar '13

Olá!Procurei no fórum e não encontrei nenhum post a respeito desse tipo de ordenação.Esse algorítimo quebra o vetor de tamanho n em subvetores e se chama recursivamente até que os subvetores tenham tamanho 1.Depois, compara dois subvetores colando na posição 1 o subvetor menor e na posição 2 o subvetor maior.E depois volta a comparar os subvetores à direita e à esquerda até ter o vetor original ordenado. Estou implementando esse algoritmo porém, não estou conseguindo dividir o vetor até que length seja = 1. mergesort=function(x){ # Função que irá dividir o subvetor ao meio subsort=function(y,z){ #Função que vai juntando ordenadamente as partes da entrada resultado=NULL tam=length(x) # Se a entrada tiver mais de um elemento, o algoritmo irá "quebrar" a entrada ao "meio" para ordenar individualmente cada entrada e depois juntar as metades já ordenadas if(tamanho>1) { meio= tam/2 esq=c[1:flor(meio)] dit=c[(meio+1):tam] esq=mergesort(esq) #vou ter que ficar aplicando o mergesort até o vetor da esquerda ter um elemento. dit=mergesort(dit) #vou ter que ficar aplicando o mergesort até o vetor da direita ter um elemento if(esq[length(esq)]<=dit[1]) #se o vetor da esquerda do tamanho da esquerda (que deve ser um) for menor ou igual a 1 { return(c(esq,dit)) #na hora de juntar o subvetor da esquerda fica primeiro. } else { subsort(esq,dit) #subsort é o vetor temporário. } } else # Se a entrada for um único elemento não existe ordenação. Se TAM<1. devolve o elemento. { return(x) } Obrigada! Att,Jackeline

4 14

Re: [R-br] Limite de confiança em uma distribuição Beta
by Juliano Farias 14 Mar '13

14 Mar '13

Prezados, Obrigado a todos pela ajuda. A função qbeta(c(0.025,0.975),7,241) não funcionou para o que preciso. Porém, consegui solucionar utilizando o pacote "binom" com a função " binom.bayes". Atenciosamente, Juliano

1 0

[Dúvida] Inversa da função gamma incompleta.
by Pedro Rafael 14 Mar '13

14 Mar '13

Pessoal existe algum pacote que tenha implementada a função inversa da função gama incompleta? Obrigado, Pedro Rafael Diniz Marinho.

1 0

Limite de confiança em uma distribuição Beta
by Juliano Farias 14 Mar '13

14 Mar '13

Prezados, Vou ser mais específico. *Preciso estimar o *limite* de confiança de **uma distribuição **beta.* * * *Estou estimando a frequência de alelos resistentes a um determinado agente de controle.* * * *Para isso uso o método de inferência **bayesiana sem prévia informação.* * * *O valor de frequência é 0,008 e os parâmetros Beta (7, 241). Preciso determinar o limite inferior e superior 95% (0,025 - 0,975).* * * *Grato pela atenção.* * * *Atenciosamente,* *Juliano*

2 1

aRT
by Gilbert Queiroz 14 Mar '13

14 Mar '13

Pessoal, estou utilizando o pacote aRT para trabalhar com dados de um db criado no TerraView. Consegui chegar até os dados das tabelas do layer. Porém, surgiu uma dúvida: é possivel modelar esses dados no R e criar uma nova coluna, ainda no R, com os resultados dessa modelagem, para ser introduzida no db original, via R ? E como seria isso??? Abs. Gilbert. PS: tem um bd chamado "peixe-boi"na página do aRT e o script correspondente que tem um comando errado...

1 0

Compilacao do fonte, dependencias e "vanilla systems"
by Benilton Carvalho 14 Mar '13

14 Mar '13

Ola', algum de voces tem algum script, colecao de comandos, guia ou qq coisa que o valha que possa auxiliar na automacao de resolucao de dependencias e compilacao do R a partir dos fontes em ambientes "vanilla" (que nao tem nada instalado) e baseados no Fedora 18? Obrigado, b

1 0

Limite de confiança em uma distribuição Beta
by Juliano Farias 14 Mar '13

14 Mar '13

Boa noite prezados, Preciso estimar o limite de confiança (95%) de um distribuição Beta. Alguém sabe algo sobre isso? Grato pela atenção. Juliano

4 3

Pacote drc
by Tiago Souza Marçal 13 Mar '13

13 Mar '13

Bom dia pessoal, gostaria de saber se alguem da lista utiliza o pacote drc para estimar a CL-50. Caso tenha gostaria que me esclarecesse as seguintes dúvidas: Ele calcula a mortalidade corrigida pelo critério de Abbott? Porque as estimativas geradas por este pacote são diferentes das geradas pelo SAS e o Polo? Att. Tiago.

2 5

sub escrito
by Ze Henrique 13 Mar '13

13 Mar '13

Pessoal como que eu deixo um texto em Italico e/ou negrito e/ou subscrito. Obrigado a todos, desde já. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

4 6

erros no codigo do R
by alanarocha＠sapo.pt 13 Mar '13

13 Mar '13

Boa tarde, alguem me pode ajudar a corrigir estes erros? DADOS: Clínicas X2008 X2009 X2010 X2011 X2012 1 PBarca 82 87 87 87 87 2 Braga 84 85 86 89 87 3 Fafe 85 87 88 86 94 4 Maia 81 86 84 90 95 5 VNGaia 88 93 92 93 95 6 SMFeira 77 83 81 84 92 # Estatística por ano porano<-apply(dados[,-1],2,function(x)c(summary(x),0)[1:7]) desvio padrão<-sd(dados[,-1], na.rm= TRUE) Error: unexpected symbol in "desvio padrão" variância<-diag(var(dados[,-1], na.rm= TRUE)) quantile(dados[,-1], probs = c(0.25, 0.75, NA)/100) ampinterquartil<-Iqr(matrix(as.numeric(dados[,-1], na.rm= TRUE))) Error in quantile(as.numeric(x), c(0.25, 0.75), na.rm = na.rm, names = FALSE, : (list) object cannot be coerced to type 'double' amplitude<-apply(dados[,-1],2,range, na.rm= TRUE) amplitudediff<-c(amplitude[2,] - amplitude[1,]) rownames(porano) <- c('Minimo','Primeiro quartil','Mediana','Média','Terceiro Quartil','Máximo','NAs') poranol<-rbind(porano,desvio padrão,variância,amplitude,amplitudediff,ampinterquartil) > colnames(porano) = 2008:2012 > poranol X2008 X2009 X2010 X2011 X2012 Minimo 66.000000 65.000000 62.000000 59.000000 60.000000 Primeiro quartil 76.000000 84.000000 86.000000 88.250000 91.000000 Mediana 80.000000 87.000000 89.500000 91.500000 94.000000 Média 80.270000 86.290000 88.500000 90.500000 91.840000 Terceiro Quartil 85.000000 89.500000 92.000000 95.000000 96.000000 Máximo 94.000000 98.000000 98.000000 98.000000 99.000000 NAs 4.000000 3.000000 3.000000 3.000000 0.000000 desviopadrão 7.362374 6.254838 6.542402 7.439656 7.879345 variância 54.204545 25.746212 21.530303 25.130682 12.780303 66.000000 65.000000 62.000000 59.000000 60.000000 94.000000 98.000000 98.000000 98.000000 99.000000 amplitudediff 28.000000 33.000000 36.000000 39.000000 39.000000 > porano 2008 2009 2010 2011 2012 Minimo 66.00 65.00 62.0 59.00 60.00 Primeiro quartil 76.00 84.00 86.0 88.25 91.00 Mediana 80.00 87.00 89.5 91.50 94.00 Média 80.27 86.29 88.5 90.50 91.84 Terceiro Quartil 85.00 89.50 92.0 95.00 96.00 Máximo 94.00 98.00 98.0 98.00 99.00 NAs 4.00 3.00 3.0 3.00 0.00 > eu preciso de ter o terceiro quartil menos o primeiro quartil e de ter mínimo menos 1,5 vezes terceiro quartil menos o primeiro quartil e o máximo mais 1,5 vezes o terceiro quartil menos o primeiro quartil ou seja obter a amplitude interquartil e os outliers. alguem tem algum exemplo? e que no grafico me apareça mínimo 1ºquartil mediana media 3ºquartil máximo e as dois tipos de amplitude. boxplot(dados[,-1], axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(main= 'indicador eKtv sup 1_2 anual') # Acrescenta o título do gráfico title(xlab= 'Ano') # Acrescenta o titulo do eixo X title(ylab= 'Percentual') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gráfico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:5, 2008:2012) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 5, as etiquetas 2008 a 2012. obrigada Ana Rocha

2 1

código do R
by alanarocha＠sapo.pt 13 Mar '13

13 Mar '13

Eu encontrei a forma de calcular os quantis que precisava mas não conheço nem percebo os erros que me deram e os resultados de "IRQ" não me aparecem. código: # Estatística por ano porano<-apply(dados[,-1],2,function(x)c(summary(x),0)[1:7]) desvio padrão<-sd(dados[,-1], na.rm= TRUE) Error: unexpected symbol in "desvio padrão" variância<-diag(var(dados[,-1], na.rm= TRUE)) IQR) <- quantile(dados[,-1], 3/4) - Error: unexpected ')' in "IQR)" quantile(dados[,-1], 1/4) Error in quantile.default(dados[, -1], 1/4) : missing values and NaN's not allowed if 'na.rm' is FALSE amplitude<-apply(dados[,-1],2,range, na.rm= TRUE) amplitudediff<-c(amplitude[2,] - amplitude[1,]) rownames(porano) <- c('Minimo','Primeiro quartil','Mediana','Média','Terceiro Quartil','Máximo','NAs') poranol<-rbind(porano,desvio padrão,variância,amplitude,amplitudediff,IQR) Error: unexpected symbol in "poranol<-rbind(porano,desvio padrão" > colnames(porano) = 2008:2012 > poranol X2008 X2009 X2010 X2011 X2012 Minimo 66.000000 65.000000 62.000000 59.000000 60.000000 Primeiro quartil 76.000000 84.000000 86.000000 88.250000 91.000000 Mediana 80.000000 87.000000 89.500000 91.500000 94.000000 Média 80.270000 86.290000 88.500000 90.500000 91.840000 Terceiro Quartil 85.000000 89.500000 92.000000 95.000000 96.000000 Máximo 94.000000 98.000000 98.000000 98.000000 99.000000 NAs 4.000000 3.000000 3.000000 3.000000 0.000000 desviopadrão 7.362374 6.254838 6.542402 7.439656 7.879345 variância 54.204545 25.746212 21.530303 25.130682 12.780303 66.000000 65.000000 62.000000 59.000000 60.000000 94.000000 98.000000 98.000000 98.000000 99.000000 amplitudediff 28.000000 33.000000 36.000000 39.000000 39.000000 > porano 2008 2009 2010 2011 2012 Minimo 66.00 65.00 62.0 59.00 60.00 Primeiro quartil 76.00 84.00 86.0 88.25 91.00 Mediana 80.00 87.00 89.5 91.50 94.00 Média 80.27 86.29 88.5 90.50 91.84 Terceiro Quartil 85.00 89.50 92.0 95.00 96.00 Máximo 94.00 98.00 98.0 98.00 99.00 NAs 4.00 3.00 3.0 3.00 0.00

2 1

cor no titulo
by Ze Henrique 13 Mar '13

13 Mar '13

Pessoal gostaria de saber se alguem sabe como mudar a cor do título do gráfico e dos eixos x e y. Obrigado. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

3 3

gráficos no R
by alanarocha＠sapo.pt 13 Mar '13

13 Mar '13

Olá, eu queria fazer um gráfico boxplot do sumário das estatísticas descritivas, tenho um exemplo mas não aparece a média no gráfico, o que terei de acrescentar ao seguinte código: boxplot(dados[,-1], axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(main= 'indicador mortalidade anual') # Acrescenta o título do gráfico title(xlab= 'Ano') # Acrescenta o titulo do eixo X title(ylab= 'Percentual') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gráfico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:5, 2008:2012) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 5, as etiquetas 2008 a 2012. para que a média apareça no gráfico juntamente com as outras estatísticas. beijinhos Ana

1 0

Formatação de números
by Edmar Caldas 12 Mar '13

12 Mar '13

Boa Tarde, Pessoal preciso de uma ajuda, tenho um arquivo em txt dessa forma ,100 preciso de algum comando que faz com que ele fique assim 0,100 com o ZERO na frente. Edmar

3 2

Criar objeto com operações distintas
by ASANTOS 12 Mar '13

12 Mar '13

Bom dia pessoal, Estou tentando obter um novo objeto com alguns dados que sera ó resultado de várias operações, sendo: # dados<-expand.grid(especie=c("sp1","sp2","sp3","sp4","sp5", "sp6","sp7","sp8"),##8 espécies de inseto dossel=rnorm(5,75,15), ## Cobertura de dossel elevacao=rnorm(5,750,50),##Elevação do terreno argila=rnorm(5,45,15),## Teor de argila presenca=rbinom(5,1,0.6)) ## Presença ou ausência do inseto # Gostaria de explicitar que a cada variação em dados$especie, fosse realizado a soma de dados$presenca e a média de dados$dossel, dados$elevacao e dados$argila, na criação do novo objeto, mas com o comando tapply(), não estou tendo sucesso em uma única linha de comando, alguém teria alguma sugestão para fazer esta operação de maneira direta, Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n 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 ======================================================================

3 8

Re: [R-br] Atualização R 2.15.3i
by Sérgio Henrique almeida da silva ju 12 Mar '13

12 Mar '13

Oi Fernando Obrigado pela ajuda. Não entendi muito bem o que você quis dizer... Abraços Em 11 de março de 2013 21:05, FHRB Toledo <fernandohtoledo(a)gmail.com>escreveu: > Use outro repositório, nesse caso o repositório da versão mais atual e não > o padrão da sua distro! > > 2013/3/11 Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> > >> Prezados >> >> Estou tentando atualizar meu R para a versão 2.15.3, porém não estou >> conseguindo. Já utilizei os tutorias anteriores aqui da lista mais não deu >> certo. >> >> Utilizo o Ubuntu 12.10 (quantal) 64 bit e meu R está na versão 2.15.1. >> >> Alguém tem alguma ideia de como posso atualizar? >> >> Att >> >> -- >> 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) 68463637 / 94429486 >> >> _______________________________________________ >> 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. >> > > > _______________________________________________ > 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. > -- 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) 68463637 / 94429486

2 1

Outra dica de compartilhamento de arquivo
by Benilton Carvalho 12 Mar '13

12 Mar '13

Pessoal, aqui vai uma outra dica (além das que já discutimos antes) de como compartilhar arquivos com outros membros da lista. A) File Dropper: http://www.filedropper.com vcs não precisam assinar... basta escolher o arquivo e fazer o upload... um link será gerado e vocês compartilham apenas o link. Nao consegui gerar um link direto para o arquivo. B) Dropbox: http://www.dropbox.com precisa de conta e é possível obter um link direto para leitura a partir do R (desde q o arquivo esteja em Public) C) Ge.tt: http://www.ge.tt nao precisa de conta e vc tbm pode obter um link direto para o arquivo (clique o botao direito sobre "Download" e use o link dado)... Por exemplo: cat(paste(readLines('http://ge.tt/api/1/files/9EuKLua/0/blob?download'), collapse='\n')) Poderia usar mais tempo descobrindo outros serviços, mas acho que por aí basta... Sintam-se à vontade para adicionar outras sugestões. b

1 0

aRT, show de bola
by Gilbert Queiroz 12 Mar '13

12 Mar '13

Pessoal, estou utilizando o pacote aRT para trabalhar com dados de um db criado no TerraView. Consegui chegar até os dados das tabelas do layer. Porém, surgiu uma dúvida: é possivel modelar esses dados no R e criar uma nova coluna, ainda no R, com os resultados dessa modelagem, para ser introduzida no db original, via R ? Abs. Gilbert.

1 0

Atualização R 2.15.3
by Sérgio Henrique almeida da silva ju 12 Mar '13

12 Mar '13

Prezados Estou tentando atualizar meu R para a versão 2.15.3, porém não estou conseguindo. Já utilizei os tutorias anteriores aqui da lista mais não deu certo. Utilizo o Ubuntu 12.10 (quantal) 64 bit e meu R está na versão 2.15.1. Alguém tem alguma ideia de como posso atualizar? Att -- 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) 68463637 / 94429486

3 3

[OFF-TOPIC] Livros de divulgação científica em estatística
by ASANTOS 11 Mar '13

11 Mar '13

Prezados membros da lista, Saberiam indicar bons livros de divulgação científica em estatística ao exemplo de Uma senhora toma chá de Salsburg e O andar do bêbado de Mlodinow, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n 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 ======================================================================

3 3

Substituir informações de base de dados
by Marcos Coque Jr 11 Mar '13

11 Mar '13

Colegas, Gostaria de uma sugestão: qual função é mais adequada para fazer a substituição de informações em uma base de dados carregada por informações de outra base de dados? Há uma variável ID que relaciona as duas bases. Seria um procedimento semelhante ao PROCV do Excel... Obrigado. Marcos

3 2

chooseCRANmirror() não funciona :(
by Pedro Emmanuel Alvarenga Americano do Brasil 11 Mar '13

11 Mar '13

Amigos de R, Há alguns dias eu venho tentando fazer algumas atualizações e não consigo. Só ontem no final da tarde eu percebi que o R não consegue ir até o espelho, mas se eu baixar o zip e instalar dá certo. Como a nova versão saiu eu achei que pudesse ser alguma coisa de pacotes de uma versão em outra. Mas tambem não é isso. Ainda, a saída do erro não ajuda muito. > chooseCRANmirror() Mensagens de aviso perdidas: In open.connection(con, "r") : não foi possível abrir: status HTTP foi '0 (nil)' Achei em mensagens antigas de outros foruns que isso tem relação de como a função chama o arquivo que contem a lisa de repositorios ou com a configuração do proxy. Eu testei getCRANmirrors e parece funcionar bem. Desativei o proxy das configurações da internet mas continuo com o mesmo erro. Há alguma maneira de fazer o R ignorar o proxy? Alguma outra dica para esse erro? Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 1

(50% OffTopic) Ajuda com análise
by Rodrigo Coster 10 Mar '13

10 Mar '13

Caros, venho pedir ajuda para vocês. O local onde trabalho realizou uma pesquisa com os servidores, sendo que a amostra foi estratificada em 1 variável (foi estratificada em 2 variáveis na verdade, mas devidos a problemas perdemos a informação de uma delas), sendo todas as perguntas objetivas. As minhas dúvidas são: 1) Que tipo de análise da para fazer em dados estratificados (a ideia é fazer cruzamentos, então preciso de algo similar ao qui-quadrado) 2) Como fazer isso no R (quais pacotes trabalham com planos amostrais complexos)? Qualquer ajuda/sugestão vai ser extremamente útil e bem vinda! []'s

4 4

para plotar biomas brasileiros
by D'Apolito 10 Mar '13

10 Mar '13

Olá pessoal, Alguém ai tem ou sabe como conseguir dados de coord. geográficas para plotar os biomas brasileiros? Por exemplo, Cerrado, Mata Atlantica, Amazonia, etc agradeço desde ja Carlos

2 1

dados DTK.test
by Hélio Gallo Rocha 10 Mar '13

10 Mar '13

Não sei se a forma de inserir os dados está correta. Anexei o fsd.txt que tem objetivo de estudo de geoestatistica, mas estou tentando rodar o teste de tukey. script: library(geoR) library(DTK) dados=read.geodata("fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "dados$data", f = "f$data", a = "0.05") grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 9

cbind sem número de linhas
by Mauricio Cardeal 09 Mar '13

09 Mar '13

Caros useRs, por favor, como faço para que não apareça no resultado os números das linhas dos vetores ([1,], [2,], etc), ao usar o cbind ? Valeu! #CMR idade <- c(2,3,6,3,1,8,6,4,3,5) peso <- c(56,65,78,65,49,120,54,76,78,88) cbind(idade,peso) idade peso [1,] 2 56 [2,] 3 65 [3,] 6 78 [4,] 3 65 [5,] 1 49 [6,] 8 120 [7,] 6 54 [8,] 4 76 [9,] 3 78 [10,] 5 88 Se uso a função data.frame() os números continuam saindo (apesar de não se tratar de vetor) mas só quero as duas colunas idade e peso. Deve ser algo bem simples, mas não consegui achar a solução. data.frame(idade,peso) idade peso 1 2 56 2 3 65 3 6 78 4 3 65 5 1 49 6 8 120 7 6 54 8 4 76 9 3 78 10 5 88 -- Mauricio Cardeal UFBA

6 9

R sendo executado de arquivos .bat ou gerenciador de tarefas do sistema.
by Pedro Emmanuel Alvarenga Americano do Brasil 09 Mar '13

09 Mar '13

Amigos de R, Estou fazendo um data management de um projeto e fiz um script para importar os dados, fazer algumas conferencias, depois algumas tabelas de controle como quantos pacientes entraram no projeto e se ha informações ausentes e coisas assim e depois colar essas coisas num arquivo de saida. Agora eu gostaria de ter um arquivo .bat ou de alguma forma agendar no sistema (no meu caso um windows XP velho amigo), para com apenas um clique, ou sem clique algum, o R abra esse script, rode todo o script, e eu tenha o arquivo com as saidas pronto para conferencia posterior. O problema é que eu consigo fazer o R abrir com um arquivo .bat ou com o gerenciador de tarefas, mas não consigo fazer o R abrir o script com o mesmo recurso. Alguem tem algum material que eu possa ler ou saberia apontar um link inspirador para essa tarefa? Abraço, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 2

usando DTK.test - Erro em input[, 1] : número incorreto de dimensões
by Hélio Gallo Rocha 09 Mar '13

09 Mar '13

Bom dia a todos. O que poderia ser esta msg de erro? Erro em input[, 1] : número incorreto de dimensões grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 2

Curso de Modelagem Estatística com R
by Gelso Pedrosi Filho 08 Mar '13

08 Mar '13

Olá! Alguma instituição teria interesse em sediar um Curso de Modelagem Estatística com R? O curso é ministrado pelos professores Graeme D. Hutcheson da Universidade de Manchester e Luiz Moutinho da Universidade de Glasgow. O curso foi realizado com sucesso na Universidade Federal de Roraima em 2011. Interessados podem solicitar informações detalhadas pelo meu e-mail. Atenciosamente, Prof. Gelso Pedrosi Filho Universidade Federal de Roraima

3 2

teste de esfericidade
by Ze Henrique 08 Mar '13

08 Mar '13

Boa a tarde a todos. Gostaria de saber como faço um teste de esfericidade pelo programa R. Pois vi em alguns lugares mas não consegui executar. Obrigado desde já a todos. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

2 2

'subsetting' de um dataframe, como fazer uma função ou loop?
by Carlos pombo sonderblohm 08 Mar '13

08 Mar '13

Bom dia pessoal, estou com uma dúvida no 'subsetting' de um dataframe de 44 variaveis ambientais, onde as primeras dos são 'MONTH' e 'YEAR' (categóricas), e o resto são numericas (algumas contem NAs) Estou a precisar usar uma função para criar um novo data frame onde as variaveis numericas [ , 4:44] sejam colocadas em promedios por 'MONTH', bom aqui vai um exemplo do CMR que estou a usar: #o data frame str(OCTOPUS_S_2003) 'data.frame': 115 obs. of 44 variables: $ ID : int 157 158 159 160 161 162 163 164 165 166 ... $ Year : Factor w/ 10 levels "2003","2004",..: 1 1 1 1 1 1 1 1 1 1 ... $ Month : int 1 2 3 4 5 6 7 8 9 10 ... $ RainFall : num 88.4 123.2 98.2 NA 8.3 ... $ PLOBOI : int 128301 245454 280876 220158 25878 2546 $ PLOBOII : num 47.9 101.46 104.87 84.94 9.66 ... $ ALQUEVA : int 395168 91474 24386 28427 53978 28278 #eu consigo fazer para cada variable, mais uma de cada vez com o seguinte código (p.e. variavel nº 5 = PLOBOI: PLOBOI <- tapply(OCTOPUS_S_2003[ ,5], OCTOPUS_S_2003$Year, mean, na.rm =TRUE) #mas são 40 variaveis tirando ID, Year e Month ##Será que posso/devo fazer um loop? o uma função? obrigado Carlos P.D.: estou a usar o R R-215~1.2\\bin\\x64\ -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

6 9

ANOVA PARA MODELOS NÃO LINEARES
by Tiago Souza Marçal 07 Mar '13

07 Mar '13

Boa tarde pessoal, eu dei uma olhada no post que o Walmes fez no dicas ridículas usando o modelo y = A*t/(V+t) + D*t e ele justificou que se V e D fossem zero este se tornaria um modelo de intercepto, ou seja, y = A*t . Entretanto, estou trabalhando com um modelo logístico y = A/(1+B*exp(-C*x)) que acredito não apresentar intercepto a não ser que fosse linearizado. Tem algum outro procedimento para obter a anova deste modelo? CMR x<-c(5,10,15,20,25) y<-c(4,7,18,31,33) Agradeço desde já por qualquer ajuda. Att. Tiago.

2 4

Como gerar um gráfico de densidade
by Gabriel Escodino 07 Mar '13

07 Mar '13

Olá, iniciei recentemente os estudos com R e surgiu uma dúvida. Tenho um data.frame com 3 colunas e um número significativo de dados da seguinte maneira: -113.568 2.91482 ARG-380 107.41 -7.87132 GLY-381 -136.334 154.558 LYS-382 -105.745 1.29007 THR-383 -156.755 164.53 GLU-384 71.2662 31.4028 GLU-385 -115.43 159.42 TRP-386 -58.5111 -51.3003 ASP-387 -104.479 111.676 VAL-388 -94.2292 109.706 PRO-389 No momento, estou gerando um gráfico onde a coordenadas X do mesmo são compostas pelos dados da primeira coluna e as coordenadas Y pelos dados da segunda. Gostaria de obter como resultado um gráfico no modelo da imagem(arquivo anexo), porém, não utilizo mapas e gostaria de saber quais são as maiores áreas de concentração dos pontos no gráfico. Alguém pode me informar se existe uma função, ou qualquer material de estudos que possa ajudar a gerar tal tipo de gráfico de densidade? Obrigado pela atenção. -- Gabriel Almeida Escodino da Silva Graduando em Ciência da Computação - UENF

2 1

Tradução para português do R Reference Card
by Walmes Zeviani 06 Mar '13

06 Mar '13

Prezados Listeiros, Venho divulgar a tradução do R Reference Card que um aluno meu, Conrado Oliveira (estatística-ufpr), fez e me passou. Fiquei entusiasmado e agradecido com a iniciativa e envolvimento dele que não poderia deixar de compartilhar o resultado dessa dedicação. Confesso que não tive tempo de revisar o trabalho com detalhe mas quero que se sintam à vontade para usar, distribuir e fazer sugestões. tradução Conrado Oliveira: http://www.leg.ufpr.br/~walmes/cursoR/guia_rapido_R.pdf original Tom Short: http://cran.r-project.org/doc/contrib/Short-refcard.pdf À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

5 4

Re: [R-br] Tradução para português do R Reference Card
by guiferraudo 06 Mar '13

06 Mar '13

Parabens aos envolvidos. Walmes parabens pela motivacao ao aluno....isso ja vale muito. Abs Enviado do tablet Samsung -------- Mensagem original -------- De : Walmes Zeviani <walmeszeviani(a)gmail.com> Data: 05/03/2013 18h55 (GMT-03:00) Para: r-br(a)listas.c3sl.ufpr.br Assunto: [R-br] Tradução para português do R Reference Card Prezados Listeiros, Venho divulgar a tradução do R Reference Card que um aluno meu, Conrado Oliveira (estatística-ufpr), fez e me passou. Fiquei entusiasmado e agradecido com a iniciativa e envolvimento dele que não poderia deixar de compartilhar o resultado dessa dedicação. Confesso que não tive tempo de revisar o trabalho com detalhe mas quero que se sintam à vontade para usar, distribuir e fazer sugestões. tradução Conrado Oliveira: http://www.leg.ufpr.br/~walmes/cursoR/guia_rapido_R.pdf original Tom Short: http://cran.r-project.org/doc/contrib/Short-refcard.pdf À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

3 2

Comando write
by Julia Maria 05 Mar '13

05 Mar '13

Pessoal, uma dúvida. Há um tempo atrás, utilizei o comando write pra ir salvando em uma planilha eletronica os resultados de uma simulação. Rodei a simulação, e fiz um vetor com os parâmetros que eu necessitava e a cada rodada o R gravava pra mim esse vetor em uma linha na planilha. Na época isso me facilitou muito. Minha pergunta é: existe algum outro comando mais interessante para se fazer isso? Eu só preciso de gravar um vetor com 5 elementos em uma planilha a cada vez que eu rodar o programa.

2 1

Erro em separadores: Ubuntu 12.04 LTS (Português)
by Vinícius Lionel Mateus 05 Mar '13

05 Mar '13

Prezados useRs, Recentemente estou usando o R no Ubuntu 12.04 LTS (64 bits). Uso o RStudio e estou com problemas para atribuir objetos. > x<- (2, 9, 10, 13) Erro: ',' inesperado em "x<- (2," Tenho o mesmo problema ao usar o Terminal. Alguém já passou por este problema? Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

2 2

teste de Bartlett
by Ze Henrique 04 Mar '13

04 Mar '13

Boa tarde a todos, gostaria de saber se alguem sabe se para teste de Bartlett quando faz para fatorial usa os dois fatores juntos ou faz para cada um dos fatores separados? Obrigado a todos deste já. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

2 2

Re: [R-br] teste de Bartlett
by Ivan Bezerra Allaman 04 Mar '13

04 Mar '13

Lembrando que todos os testes implementados em softwares que existem são one-way. Em 2007 se não me engano propuseram uma extensão do teste de Levene para análises multi-way. (S,f,P) Allaman

1 0

Gerar mapa de distribuição por UF
by Kaue Veras 04 Mar '13

04 Mar '13

Bom dia Senhores (as), Eu tenho uma base com três colunas e mais de 100000 linhas que contém a UF o CEP e a Quantidade de clientes por cada CEP/UF, eu gostaria de saber se alguém conhece algum pacote ou meio que eu possa inserir esses dados numa distribuição por UF no mapa do Brasil. Desde já agradeço. Atenciosamente, Kauê P. Veras da Cunha Estatística - UERJ Cel: (21) 8254-9601

3 12

Usando o termo Error() com model.matrix para simular dados de Medidas Repetidas
by Augusto Ribas 04 Mar '13

04 Mar '13

Ola, normalmente para tentar entender melhor algumas analises, eu tento gerar exemplo e ver como ficam os resultados e gráficos. Normalmente eu uso o seguinte método: ################################################################### #gerando os tratamentos example<-data.frame(time=factor(rep(c(0,1,2),each=60)), trat=rep(rep(c("controle","dose2","dose3"),each=20),each=60), id=factor(rep(1:60,3))) str(example) #Usando model.matrix para fazer a matriz que multiplica os parâmetros, se quiser simular uma "two way" matrix.1<-model.matrix(~time*trat,data=example) head(matrix.1) #Gerando valores dos parâmetros param.1<-c(1:9) #Gerando os dados de resposta, que é a multiplicação da matriz pelos valores dos parâmetros response.1<-rnorm(nrow(matrix.1),matrix.1%*%param.1,1) #Assim a analise vai estimar os parâmetros usados, e eu posso ver vários parâmetros diferentes e o que acontece. summary(lm(response.1~example$time*example$trat)) #Mas eu gostaria de simular algo com medidas repetidas, tipo usando o termo Error() que da para usar com a função AOV. #Mas eu fiquei meio perdido como proceder, para fazer o mesmo esquema acima mas com medidas repetidas, tendo o error term. matrix.2<-model.matrix(~time*trat+Error(id),data=example) #Olhando como se faz a mesma coisa com a função gls do pacote nlme, suspeito que teria que fazer uma matriz de correlação, mas ai eu não sei como faria aquela multiplicação la em cima. #Alguém saberia como proceder? ##################################################################################### Os comando para fazer a medidas repetidas com nlme que me refiro são desse tutorial: http://www.ats.ucla.edu/stat/r/seminars/Repeated_Measures/repeated_measures… Olhando isso e a função aov eu vejo que tenho que separar o Error, mas não estou sabendo como. library(nlme) longg <- groupedData(pulse ~ as.numeric(exertype) * as.numeric(time) | id, data = exer) fit.cs <- gls(pulse ~ exertype * time, data = longg, corr = corCompSymm(, form= ~ 1 | id) ) summary(fit.cs) Agradeço a atenção Abraços Augusto Ribas -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

1 0

[OT] Um bom livro de estatística
by Rodolfo Timoteo da Silva 03 Mar '13

03 Mar '13

Senhores, sou graduando em Engenharia de Produção e estou interessado em comprar um bom livro impresso de estatística//que tenha os detalhes e "por quês" matemáticos. Se existir alguma publicação com a utilização do R (pt_BR ou en), mais interessante ainda. Qualquer opinião é bem vinda. Att

3 2

Dados com distribuição espaço-temporal no R
by Gilbert Queiroz 02 Mar '13

02 Mar '13

Amigos, Como eu entro com os dados no R para obter uma modelagem de dados espaço-temporais???? Tenho um arquivo onde nas colunas estão as variáveis que quero trabalhar e nas linhas os municípios. Só que eu repito estes municípios para os anos de 2002 a 2011. É assim mesmo???? Abs. Gilbert.

1 0

simulação de dados
by Sérgio Henrique almeida da silva ju 02 Mar '13

02 Mar '13

Olá Como posso simular 2000 amostras aleatórias com 2000 sujeitos cada uma? Quero modelar um desfecho binário (óbitos e não óbitos). A distribuição dos meus dados é: ### Meus Dados id <- rep(1:2000) servidor <- c(rep('Apoio',292),rep('Tecnico',776), rep('Superior',932)) participacao <- c(rep(0,600),rep(1,1400)) # 0 = Não participou do estudo e 1 = Participou do estudo obitos <- ????? Como posso distribuir meus óbitos em relação a variável participacao, sendo uma chande de 2,5 de óbitos de quem não participou em comparação a quem participou? Obrigado -- 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) 68463637 / 94429486

2 2

summary de GLMM
by Renata Carvalho 01 Mar '13

01 Mar '13

Bom dia a todos, estou fazendo um glmm com efeito randômico usando a função lmer do pacote lme4, mas quando dou o summary do modelo não consigo verificar se há "overdispersion" nos dados. Por acaso alguém sabe como verifico isso? obrigada Renata

1 0

Quadrado médio de duas características QM(y1+y2)
by Tiago Souza Marçal 01 Mar '13

01 Mar '13

Boa noite pessoal, estou querendo estimar a correlação genética entre duas variáveis (y1 e y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM (a+b). Alguém sabe como estimar este quadrado médio para as duas variáveis no R. Att. Tiago.

5 11

Re: [R-br] Resolvendo uma distribuição Uniforme SIMPLES.
by Ivan Bezerra Allaman 28 Fev '13

28 Fev '13

Eu não entendi muito bem o que vc quer realmente, mais você pode fazer algo deste tipo: http://nbcgib.uesc.br/lec/prof/ivan/sim/fdpu Para saber como fazer o que está no link abaixo acesse: http://www.rstudio.com/shiny/ (S,f,P) Allaman

2 1

Gerar números aleatórios de um dist. T mutivariada
by Carlos Pereira 28 Fev '13

28 Fev '13

Bom dia pessoal, Queria saber se tem um pacote que gera números aleatórios de uma distribuição T multivariada. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

3 3