Abril 2011 - R-br - DINF UFPR

(sem assunto)
by Juraci Mendes Moreira 15 Abr '11

15 Abr '11

Estou tentando determinar um tamanho otimo de parcela com base em um experimento em branco, usando o modelo plato de resposta linear, mas deu esse erro e não sei como resolver. Usando a mesma rotina em analises anteriores deu certo. mas essa não. sem mais para o momento antecipo agradecimentos. jura Erro em qr.qty(QR, resid) : NA/NaN/Inf em chamada de função externa (argumento 5) Além disso: Warning message: In Ops.factor(lhs, rhs) : - not meaningful for factors > lrp.fit

2 1

Teste de t^2 multivariado de Hotellings
by Alexandre dos Santos 15 Abr '11

15 Abr '11

Boa tarde pessoal, Estava tentando comparar dois grupos de indivíduos com 5 variáveis, usando o teste de t^2 de Hotelling, ate ai tudo bem, mais quando simulei os indivíduos e suas variáveis com pelos menos um NA no vetor, a função não conseguia realizar o cálculo, dando a mensagem Erro em solve.default(v) : sistema é computacionalmente singular: condição recíproca número = 0 Tentei inserir na função na.omit(), sem sucesso, alguem poderia me dar um help, obrigado, segue a rotina: #####################Teste de t multivariado de Hotellings############################################ ##Depois de algum teste de cluster (ex. Mahalanobis) vou observar se os dois grupos possuem diferenças ##Grupo 1 ind.pop.1<-1:30 #######Os indivíduos med.pop.1<-data.frame(mma=Mod(rnorm(30)),mmc=Mod(rnorm(30)),mmd=Mod(rnorm(30 )),mca=Mod(rnorm(30)),mca2=Mod(rnorm(30)))#######As medidas dos indivíduos dados.pop.1<-cbind(ind.pop.1,med.pop.1)#######Planilha com os indivíduos e suas medidas head(dados.pop.1) ##Grupinho Oligo 3 ind.pop.3<-1:6 #######Os indivíduos med.pop.3<-data.frame(mma=Mod(rnorm(6)),mmc=Mod(rnorm(6)),mmd=Mod(rnorm(6)), mca=Mod(rnorm(6)))#######As medidas dos indivíduos mca2<-c(1.22222,1.666666,0.4500000,NA ,2.33333,0.789000) dados.pop.3<-cbind(ind.pop.3,med.pop.3,mca2)#######Planilha com os indivíduos e suas medidas head(dados.pop.3) hotelling = function(y1, y2) { ### Função que faz o cálculo de Hotellings by http://sas-and-r.blogspot.com/2010/05/example-737-calculation-of-hotellings. html # check the appropriate dimension k = ncol(y1) k2 = ncol(y2) if (k2!=k) stop("input matrices must have the same number of columns: y1 has ", k, " while y2 has ", k2) # calculate sample size and observed means n1 = nrow(y1) n2 = nrow(y2) ybar1= apply(y1, 2, mean); ybar2= apply(y2, 2, mean) diffbar = ybar1-ybar2 # calculate the variance of the difference in means v = ((n1-1)*var(y1)+ (n2-1)*var(y2)) /(n1+n2-2) # calculate the test statistic and associated quantities t2 = n1*n2*diffbar%*%solve(v)%*%diffbar/(n1+n2) f = (n1+n2-k-1)*t2/((n1+n2-2)*k) pvalue = 1-pf(f, k, n1+n2-k-1) # return the list of results return(list(pvalue=pvalue, f=f, t2=t2, diff=diffbar)) } ## Fazendo o teste entre o grupo 1 e 3 (result= hotelling(dados.pop.1,dados.pop.3)) Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29

1 0

Gráficos de superfície!!
by Thiago De paula protásio 15 Abr '11

15 Abr '11

Pessoal; Estou tendo dificuldades em colocar no gráfico de superfície tonalidades de cor diferenciada em cada parte do gráfico, pois isso facilita a visualização e interpretação da relação de cada variável independente com a variável independente. Estou usando os seguintes comandos: TCF=seq(68,80) RG=seq(32,37) f1=function(x,y)3696.046+62.183*x-33.414*y PCS=outer(TCF,RG,f1) persp(TCF,RG,PCS,theta=-50,phi=-0,shade=2,ticktype = "simple",box=T,border=T,r = sqrt(20), d = 1,col="blue") Agradeço a ajuda! Atenciosamente, Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Data - time Zones
by eder＠leg.ufpr.br 15 Abr '11

15 Abr '11

Bom dia pessoal, Eu tenho datas no seguinte formato: DATE <- "00Z29Sep2006" # hora Z dia mês ano # Eu consegui mudar o formato da data com o comando: Sys.setlocale("LC_TIME", "USA") DATE <- strptime(paste(substr(DATE,4,12),substr(DATE,1,2),sep=''),"%d%B%Y%H") Sys.setlocale("LC_TIME","Portuguese_Brazil.1252")## Volta ao original DATE Meu objetivo final e ter a data no formato dd-mm-aaaa hh:mm:ss ('13-04-2011 08:35:00') Alguma tem uma dica, de como mudar o formato desta data direto sem ter que usar o Sys.setlocate? Atenciosamente Éder

5 10

Dúvida sobre função
by Victor Eduardo 15 Abr '11

15 Abr '11

Pessoal, há a possibilidade de criar uma única função que tenha varias váriaveis sendo que depois dessa função eu queria que cada variavel analisada da função gerasse uma coluna para cada uma, por exemplo, os comandos seriam mais ou menos esses: # Diretório wrk.dir <- "/home/nathalia/Área de Trabalho/synop/" dat.dir <- "/home/nathalia/Área de Trabalho/synop/" setwd(wrk.dir) # Leitura dos dados coluna<-c("sid", "ano", "mes", "dia", "hora", "hh", "lon", "lat", "height", "hls", "ps", "pmsl", "dd", "ff", "t", "td") dados<-read.table('201103.txt',head=T,col.names=coluna) attach(dados) # Criar uma nova coluna para a planilha dados #dados<-transform(dados, temp=t-273) fx<-function(var.temp, var.pre){ ifelse(var.temp> -25 & var.temp <40 | var.pre>930 & var.pre<1050, 1, 0) } dados<-transform(dados, CQ_temp= fx(t-273)) attach(dados) No caso, está errado a minha função, pois estou usando o | e isso não vai ficar legal, pois acabam que tornam elas dependentes, e quero analisar de forma independente. Agora na hora do transform, eu queria ver se dava para pegar cada variavel da minha função e ir adionando os 0 e 1 em coluna das respectivas variaveis. É possível isso? Abraços!

2 1

Arquivo Binario
by eder＠leg.ufpr.br 14 Abr '11

14 Abr '11

Pessoal, Minha duvida agora e com arquivos binarios, contextualizando: Quero abrir arquivos de saida de modelos climáticos (Eta), estes modelos soltam arquivo no formato GRIB, em em palavras R-anas uma array onde cada variável meteorológica seria um z dento do array (dados[x,y,z]), no arquivo ctl (indica como o binario é) indica qual o código falha (NA), em arquivo que o código era -9999 a função readBin funcionou, em outro arquivo que o código era 1.0E20 ela não funcionou ou seja estou suspeitando que ela não esta entendo este código falha, exemplo: formato <- readLines('http://www.leg.ufpr.br/~eder/ll20gfsa2006120100+2006120100.ctl') formato cone <- file('http://www.leg.ufpr.br/~eder/ll20gfsa2006120100+2006120100.dat','rb') dados <- readBin(cone,numeric(),195*200*58,4) length(dados) head(dados) close(cone) pelas informações em formato era para ter no arquivo binário 195*200*58= 2262000 números, porem ela carrega apenas 1911000, faltando dados, desta forma. Alguma dica? Atenciosamente

3 6

Ajustar escala em plot
by Marcelo Claro de Souza 14 Abr '11

14 Abr '11

Estou plotando alguns gráficos de regressões lineares da seguinte forma: cor=read.csv2('correlações.csv',h=TRUE) cor=as.data.frame(cor) simbolos=c(19,19,19,19,19,19,19,2,2,2,2,2,2,2,3,3,3,3,3,3,3) N.P.cor=lm(N~P,data=cor) plot(N~P,data=cor,pch=simbolos,ylab=expression(N~(g~Kg^{-1})),xlab=expression(P~(g~Kg^{-1}))) abline(N.P.cor) legend("topleft",c("Ep","Es","Et"),pch=c(19,2,3)) Como eu faço para delimitar a escala dos eixos x e y deixando por exemplo x(0:2) e y(0:30) Muito obrigado Marcelo Claro de Souza Biólogo, Doutorando em Biologia Vegetal, Depto. de botânica, IBRC- UNESP.16 - 91914780 http://lattes.cnpq.br/5691388412841376

3 2

Res: configuracoes do RStudio
by Samuel Carvalho 14 Abr '11

14 Abr '11

Obrigado Lucas. O problema era mesmo a versão do RStudio que estava utilizando, tanto que foi discutido isso no blog do software Benilton o problema não era como sugerido na lista Valeu ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 1

Sugestões
by Alexandre Galvão Patriota 14 Abr '11

14 Abr '11

Caros, Gostaria de parabenizar aos organizadores pela nova lista. Só tenho alguma sugestão: Existem muitas informações, o que é importante para os leitores. Porém, algumas informações importantes poderiam ser mais diretas e de fácil acesso sem precisar ler todo o conteúdo. Por exemplo, a senha de acesso. Está no meio do texto e uma batida de olho rápida eu não consegui acha-la. Só depois de ler o conteúdo quase inteiro do email eu percebi qual era a minha senha. Poderia ser uma coisa direta: Login: XXX Senha: YYY Quem não quiser ler o texto todo não precisará. Isso se aplica a outras informações importantes também. Fica a dica, Abraço a todos, Alexandre.

1 0

Intercept
by Edson Lira 14 Abr '11

14 Abr '11

Boa tarde a todos. Como faço para retirar o intercep do modelo de regressão, não consegui ver isso na ajuda online. Edson Lira Estatístico Manaus-Amazonas

3 2