- R-br - DINF UFPR

Informação arquivo xls - Detalhes
by Eder David Borges da Silva 24 Abr '12

24 Abr '12

Boa noite pessoal, Venho mais uma vez recorrer as suas sabias ajudas, em: http://www.leg.ufpr.br/~eder/teste.xls deixei um arquivo xls, ###Caso queiram carregar ele: download.file('http://www.leg.ufpr.br/~eder/teste.xls','teste.xls') file.info('teste.xls') # Minha duvida é a seguinte, se clicarmos com o direito e ir nas abas: # detalhes -> descrição # Veremos que tem a opção de ter um titulo, assunto, marca etc. # minha missão é modificar este titulo (esta como MEU TITULO) para outro nome. # via comando... # Alguma sugestão? # condições de contorno: SO Windows, R 2.15. Obrigado pessoal

2 2

instalação de pacotes
by André Figueiras Rutz 24 Abr '12

24 Abr '12

Caros, Estou tendo dificuldades com o R na instalação de pacotes alguem poderia me ajudar? Parece que o pacote vai para o micro mas de alguma forma na hora de carregar nao funciona sugestões? Grato André Rutz > install.packages("rJava") --- Please select a CRAN mirror for use in this session --- tentando a URL ' http://cran-r.c3sl.ufpr.br/bin/windows/contrib/2.15/rJava_0.9-3.zip' Content type 'application/zip' length 746108 bytes (728 Kb) URL aberta downloaded 728 Kb package ‘rJava’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\user\AppData\Local\Temp\RtmpoBqDjS\downloaded_packages > require(rJava) Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry > require("rJava") Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry

3 4

Como saber se é um 'bug' ou uma 'feature', no comportamento de 'plot'?
by Cleber N.Borges 24 Abr '12

24 Abr '12

olá a todos da lista, estava fazendo uns testes e me deparei com um comportamento do qual *acho* que pode ser incoerente. Vamos ao problema: > par(mfrow=c(2,2)) > plot(1:10) # grafico 01 > plot(1:10, cex.axis=0.5, las=1) # grafico 02 > plot(table(1:10)) # grafico 03 > plot(table(1:10), cex.axis=0.5, las=1) # grafico 04 a mudança dos parametros: 'cex.axis' e 'las' funcionam muito bem na mudança do gráfico 1 para 2, mas não há efeito na mudança de 3 para 4. Segunda pergunta: Como saber qual é a verdadeira função 'plot.ALGUMACOISAAQUI' que será chamada para fazer o grafico de um objeto da classe 'table'? dado que não há 'plot.table' como função... Grato pela atenção de todos. Cleber

2 3

Dúvida sobre ANOVA tipo III
by Olympio Neto 24 Abr '12

24 Abr '12

Pessoal, Tenho lido vários livros sobre R, muitos blogs (R-bloggers) e as excelentes páginas do Wiki-Leg da UFPR e as famosas Ridículas do Walmes e em muitas delas vejo que em muitos dos tutoriais, é necessária, para a obtenção da Anova Marginal (Tipo III), setar o contraste como "contr.sum" e outros que nada falam sobre isso, ficando setado o padrão do R, o "contr. treatment". A explicação para a alteração do contraste é que a soma deste deve ser nula, o que não acontece com o "contr.treatment". Caso não seja feita a alteração, os resultados podem estar errados. O artigo que tenta por luz neste assunto é aqui indicado: http://goanna.cs.rmit.edu.au/~fscholer/anova.php Normalmente eu uso o pacote "car" para cálculos que análises de variância do tipo III, mas nos exemplos que vejo, dados que simulo etc, não vejo diferenças. Acredito que nem poderiam haver, mesmo. A ANOVA trivial (tipo I) é invariante à parametrização, mas fico muito em dúvida se a tipo III também seria. O que vocês pensam? Já observaram diferenças importantes em suas análises? Abs e muito obrigado. Olympio

4 10

Paste
by Ricardo Fernandes Campos Junior 24 Abr '12

24 Abr '12

Olá, pessoal. Eu sei que essa pergunta é provavelmente muito simples, mas eu não consigo encontrar como fazê-lo. Eu tenho um objeto de caracteres de tamanho genérico, para fim de exemplo, vou usar 3. Eu quero criar um outro objeto com estes mesmos caracteres, mas concatenados. Porém não consigo. Tentei isso: exemplo = c("A", "B", "C") novoobjeto = paste(exemplo, " ") O que me é retornado é isso: "A" "B" "C" Mas o que eu quero é isso: "A B C" Alguém pode ajudar? -- Abraços, Ricardo Fernandes Campos Junior Departamento de Genética e Biologia Evolutiva Universidade de São Paulo

3 3

Notação
by Edson Lira 24 Abr '12

24 Abr '12

Tenho esse script abaixo: require(PASWR) trat<-rep(LETTERS[1:4],5) med<-c(25,31,22,33,26,25,26,29,20,28,28,31,23,27,25,34,21,24,29,28) exem1_4<-data.frame(med,trat) exem1_4 m1<-lm(exem1_4$med~exem1_4$trat) anova(m1) summary(m1) erro<-resid(m1) erro est<-predict(m1) med-est Quando gero o erro ele me dá essa saída. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 2.00e+00 4.00e+00 -4.00e+00 2.00e+00 3.00e+00 -2.00e+00 -3.33e-16 -2.00e+00 -3.00e+00 1.00e+00 2.00e+00 -1.11e-16 8.88e-16 1.11e-16 -1.00e+00 3.00e+00 -2.00e+00 18 19 20 -3.00e+00 3.00e+00 -3.00e+00 Gostaria de transformar estes números para a notação comum, tipo 2, 3, -4 .... Me lembro que certa vez alguém postou esse problema, mas faz muito tempo. Alguém tem alguma sugestão? Edson Lira Estatístico Manaus-Amazonas

3 2

[Dúvida] Números aleatórios.
by Pedro Rafael 23 Abr '12

23 Abr '12

Formula em LaTeX: F_{T}(t;\alpha,\beta) = \Phi(\frac{1}{\alpha}(\frac{t}{\beta}^(1/2) - \frac{\beta}{t}^(1/2))) Pessoal com gerar números aleatórios com distribuição Birbaum-Saunder em C? Tenho a distribuição acumulada da Birbaum-Saunder que é escrita como acumulada de uma normal e tenho a densidade da Birbaum. A formula da acumulada está escrita acima para t e beta maiores que zero. Pensei primeiramente em pegar números aleatórios normais pela biblioteca GSL e fazer ser igual à t na formula acumulada acima. Daria certo? Se for pela transformada inversa eu teria que gerar números uniformes[0,1] e fazer F(t) = U. Depois eu tinha que ver qual o t em F me daria U, sendo U o número uniforme gerado. Outra complicação é que tenho uma acumulada da normal que é uma integral de -inf a t. Existe outra forma mais fácil? -- Saudações, Pedro Rafael Diniz Marinho.

2 1

Ajuda para evitar um loop com elegância 2
by Gustavo Marcatti 23 Abr '12

23 Abr '12

Olá pessoal, tenho que executar uma função parecida com essa no R. Mas ela normalmente demora longos 500 segundos para executar. Fiz uma versão que chama um código C, e essa é bem eficiente (demora cerca de 0,5 segundos), porém gostaria de uma função totalmente em R. Caso alguém queira ajudar, desde já agradeço. set.seed(1) vetor1 <- rep(1:4, each = 5) var1 <- abs(rnorm(4,3)) * vetor1 cond1 <-abs(rnorm(4,10)) * vetor1 var2 <- matrix(0, nrow=length(vetor1)) cond2 <- matrix(0, nrow=length(vetor1)) for(i in 1:(length(vetor1) - 1)){ if(vetor1[i] == vetor1[i+1]){ var2[i] = var1[i+1] cond2[i] = cond1[i+1] }else { var2[i] = NA cond2[i] = NA } } var2[length(vetor1)] <- NA cond2[length(vetor1)] <- NA dados <- na.exclude(data.frame(var2, cond2)) Gustavo Marcatti Eng. Florestal

2 2

R servidor Windows?
by Vinicius Brito Rocha 23 Abr '12

23 Abr '12

Senhores, é possível instalar o R em um servidor Windows para a utilização de multiplus usuários? Att. -- 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 2

bootstrap
by Maria Papa 23 Abr '12

23 Abr '12

Olá Pessoal, boa noite Estou tentando, sem sucesso, realizar uma reamostragem bootstrap para médias de 3 amostras estão organizados em uma matriz na qual, que sendo cada coluna representa uma amostra. Porém, a cada reamostragem, eu preciso que seja feita para as três amostras, ou seja, quando um dado for selecionado para formar a amostra bootstrap, deve ser retirada a linha toda. O resultado final deve ser 3 médias reamostradas, que foram obtidas em um único processo de reamostragem. Alguém sabe como fazer isso no R? Eu consegui somente quando eu realizo 3 procedimentos de reamostragem distintos. Desde já agradeço a todos célia

3 2