- R-br - DINF UFPR

RES.: Função runs.test! Falta de Ajuste!!!!
by Thiago De paula protásio 09 Abr '11

09 Abr '11

Walmes e demais membros da lista; Analisei aqueles dados conforme vc sugeriu.... Os dados são os mesmos que eu enviei a lista recentemente..... No meu caso eu tenho repetições......e os resultados que obtive com o teste de falta de ajuste sugere que o modelo não linear não está adequado. Por favor, veja o que vc acha! Caso vc e os demais membros da lista queiram reproduzir....basta utilizar os dados que enviei... As variáveis são as mesmas!! m0 <- lm(UEQ~factor(UR)) > summary(m0) Call: lm(formula = UEQ ~ factor(UR)) Residuals: Min 1Q Median 3Q Max -1.1275 -0.2663 0.0900 0.3284 0.6525 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 5.3113 0.1631 32.556 < 2e-16 *** factor(UR)40 -0.3262 0.2307 -1.414 0.164 factor(UR)50 1.8825 0.2307 8.159 1.09e-10 *** factor(UR)60 2.9800 0.2307 12.916 < 2e-16 *** factor(UR)70 3.4750 0.2307 15.062 < 2e-16 *** factor(UR)80 7.8688 0.2307 34.106 < 2e-16 *** factor(UR)90 16.9163 0.2307 73.321 < 2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.4614 on 49 degrees of freedom Multiple R-squared: 0.9941, Adjusted R-squared: 0.9934 F-statistic: 1372 on 6 and 49 DF, p-value: < 2.2e-16 > anova(n0,m0) # teste da falta de ajuste Analysis of Variance Table Model 1: UEQ ~ A * (1 - exp(B * UR)) Model 2: UEQ ~ factor(UR) Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F) 1 54 202.799 2 49 10.433 5 192.37 180.69 < 2.2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Muito Obrigado! 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

RES.: Função runs.test()!!!
by Thiago De paula protásio 09 Abr '11

09 Abr '11

Walmes e demais membros da lista; Com base naquelas dados que eu envie a lista ....Ajustei o modelo abaixo e gostaria de testar se o mesmo é compatível...... Também ajustei um polinômio de ordem 2 e fiz uma anova() entre os dois modelos....Vc acha que esse procedimento é válido? Dessa forma, estava pesquisando e encontrei esse teste de sequências da função runs.test(). Dê uma olhada nos comandos abaixo..... n0 <- nls(UEQ~A*(1-exp(B*UR)), data=da, start=list(A=-2, B=0.03)) summary(n0) # curva ajustada A <- coef(n0)["A"] B <- coef(n0)["B"] plot(UEQ~UR, da) curve(A*(1-exp(B*x)), add=TRUE, col=2,lty=1) title("Tratamento 1 ") SQE <- summary(n0)$sigma^2*summary(n0)$df[2] SQE SQT <- var(da$UR)*(length(da$UR)-1)# soma de quadrado total corrigida R2 <- 1 - SQE/SQT R2 erro1=residuals(n0) #pacotes necessários.... library(tseries) library(quadprog) library(zoo) runs.test(factor(sign(erro1))) # permite testar se os dados diferem da curva ajustada ########Modelo polinomial model=lm(UEQ~UR+I(UR^2)) summary(model) anova (model, n0) Muito obrigado pelas dicas ..vou olhar aqui no banco de teses da UFLA e encontrar a sua dissertação..... Valeu mesmo..... Suas dicas tem me ajudo bastante! 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

Função runs.test()!
by Thiago De paula protásio 09 Abr '11

09 Abr '11

Prezado membros da R-Br; Gostaria de perguntar uma coisa a vcs: O teste de sequência da função runs.test() é confiável para verificar se os dados diferem significativamente da curva ajustada? No meu caso uma função não linear..... Ou existe uma maneira melhor de verificar isso? Obrigado pela ajuda! 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

Correlação de Spearman
by Rafael Felice Fan Chen 09 Abr '11

09 Abr '11

Não sei se me expliquei bem, a minha intenção é verificar se existe alguma relação do preço pago pelo sexo, os dados não tem uma distribuição normal, e estão assim tabulados preco sexo 98 1 99 2 98 1 97 2 ...... Não é apropriado usar a correlação de spearman? Qual a interpretaçao do resultado do teste de Wilcox? Obrigado

1 0

Converter o vetor de slot de lista para que cada elemento do vetor seja um slot dela
by Walmes Zeviani 09 Abr '11

09 Abr '11

Saudações a todos, Estou montando um procedimento em que numa certa altura do código eu preciso de uma lista. Eu possuo um vetor nomeado que quero incluir nessa lista, mas que cada elemento do vetor seja um slot da lista, sendo o primeiro slot é x NULL e o ultimo o corpo de uma função. Abaixo tento via código esclarecer o meu objetivo. #----------------------------------------------------------------------------- # uma função fn <- function(x){ sen(x) } #----------------------------------------------------------------------------- # um vetor vec <- c(a=1, b=2, c=3, d=4) #----------------------------------------------------------------------------- # a lista que sei até agora obter(preciso converte-la) list(x=NULL, vec, body(fn)) $x NULL [[2]] a b c d 1 2 3 4 [[3]] { sen(x) } #----------------------------------------------------------------------------- # como faço para obter essa lista aqui, em que cada elemento do vetor # é um slot da lista, ficando assim, mas convertendo a lista de cima? # essa abaixo é a lista no formato de que preciso list(x=NULL, a=1, b=2, c=3, d=4, body(fn)) $x NULL $a [1] 1 $b [1] 2 $c [1] 3 $d [1] 4 [[6]] { sen(x) } #----------------------------------------------------------------------------- Grato à todos. 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 5

RES.: Função sigmóide! Função nls.lm
by Ivan Bezerra Allaman 09 Abr '11

09 Abr '11

Walmes, Não generalise. Eu falei especificamente da função "nls.lm". Acredito que pelos comentários vc não leu sobre tal função. Assim como vc, também trabalhei com modelos não-lineares durante toda minha tese. Com tal função, vc não chuta qualquer coisa e não obtém qualquer coisa. Muito pelo contrário! Vc não precisa chutar nada. O algorítmo proposto pelos autores é extremamente confiável e converge qualquer coisa. Então, pode surgir a pergunta: Mas se ela converge qualquer coisa, é confiável? Sim! Pois se o modelo ajustado pelo função não descreve realmente o comportamento dos dados, no momento que utilizar outras funções para fazer inferências (confit, summary, entre outras), a função retornará um erro, o que indica que a função "apenas" converge. Sabemos muito bem que para encontrar estimativas dos estimadores nenhuma pressuposição estatística precisa ser feito. Logo, com esta função em mãos, é completamente desnecessário brincar de encontrar chutes iniciais. Portanto, fiquem tranquilos, por mais que a função encontre valores sem sentido, no momento de realizar inferência ela lhe retornará um erro, lhe indicando que aquelo modelo não é o correto. Ainda utilizo a função "nls" como já comentei, em casos no qual preciso utilizar o argumento "weights", do contrário, é muita perca de tempo, pois as 30 horas que gastastes para encontrar os chutes iniciais para utilizar a função "nls" (já perdi este tempo no paleozóico), a função "nls.lm" encontra em poucos segundos com a mesma precisão da função "nls". Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Correlação de Spearman
by Rafael Felice Fan Chen 09 Abr '11

09 Abr '11

Bom dia Estou tentando fazer uma correlação de spearman entre o preço ( variavel quantitativa continua) e o sexo do animal (variavel qualitativa nominal, nesse caso usei 1 para macho e 2 para femea) e estão pareados. Usei o seguinte comando cor(dados$preco, dados$sexo, method="spearman") E o R me deu um valor, a minha dúvida se está certo o jeito que eu fiz e como gerar um p-value? Rafael Felice Fan Chen TEL - 041 9231-9056 MSN rafaelchen(a)hotmail.com SKYPE- rafaelffchen

3 2

Re:Função sigmóide! Dados!
by Ivan Bezerra Allaman 08 Abr '11

08 Abr '11

Guilherme! Não quebre muita a cabeça com a função "nls". A pouco tempo, o Gustavo se eu não me engano, passou a função "nls.lm" da library "minpack.lm". Essa função é fantástica, pois não precisamos ficar quebrando a cabeça com chutes iniciais. Eu só utilizo a função "nls" nos casos em que preciso utilizar o argumento "weights", pois na função "nls.lm" ainda não encontrei algo parecido. Teste esta função, tenho certeza que irá gostar. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

RES.: Função Sigmóide! Obrigado!
by Thiago De paula protásio 08 Abr '11

08 Abr '11

Walmes; Muito obrigado pela imensa ajuda..... Sua didática foi fenomenal....agora compreendi! Muito obrigado! Abraço!

1 0

RES.: Função sigmóide! Dados!
by Thiago De paula protásio 08 Abr '11

08 Abr '11

Walmes e demais amigos da lista; Estou enviando os dados (link abaixo) conforme o Walmes sugeriu... Está em txt e a variável UR (umidade relativa) e a variável UEQ (umidade de equilíbrio)....A ideia é tentar o ajuste de um modelo do tipo sigmóide entre a UEQ~UR. http://www.datafilehost.com/download-a7e94384.html Estou tentando ajustar o seguinte modelo.....conforme comando abaixo, mas estou tendo dificuldades em determinar os parâmetros iniciais. n0 <- nls(UME~SSlogis(UR, Asym, scal, xmid),start=list(Asym=5,scal=1,xmid=0.5),data=umi) Erro em nlsModel(formula, mf, start, wts) : matriz gradiente singular nas estimativas iniciais dos parâmetros Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf Talvez, com esses dados que enviei vcs possam me ajudar...... Muito obrigado ao Walmes e demais membros da lista!

3 2