- R-br - DINF UFPR

Cálculo numérico no R
by Gilbert Queiroz 14 Set '13

14 Set '13

Pessoal, Alguém sabe indicar algum material sobre Cálculo numérico usando o R? Abs

1 0

Gráfico (Intervalo com números decimais)
by Eric Leandro 14 Set '13

14 Set '13

Boa noite pessoal, Me chamo Eric e recentemente comecei a usar o R e tenho o seguinte problema: Quero plotar um scatterplot com intervalo no eixo Y de 0.0 a 0.6. No entanto com o comando que estou tentando o intervalo que aparece é de -1 a 1. O comando que uso é o seguinte. plot(ph,phpredito, xlab="pH medido",ylab="pH predito",pch=16,xlim=range(0:5),ylim=range(0:0.6)) Att Eric Leandro Silva Pereira Mestrando em Clima e Ambiente Instituto Nacional de Pesquisas da Amazônia Programa de Grande Escala da Biosfera-Atmosfera na Amazônia Av. André Araújo, 2936 - Aleixo, CEP 69060-001, Manaus - AM.

1 0

ordenar valores de uma matriz
by Aline Lipsky 13 Set '13

13 Set '13

Boa tarde. Eu preciso de uma ajuda. Tenho a seguinte matriz: 2011 2012 2013AC 59.1 30.4 69.6DF 0.4 0.7 1.4EP 33.3 54.8 23.1NPEP 7.2 14.1 5.8 Gostaria de ordenar os dados em ordem decrescente para cada variável (anos de 2011,2012 e 2013). Grata, Aline

2 2

matematica com o R
by Samuel Carvalho 12 Set '13

12 Set '13

Caros, boa tarde! Alguém poderia por gentileza auxiliar com a duvida relativa ao codigo abaixo? # RCMR y <- expression(exp(a+b/x)) D(y,'x') #primeira derivada dy.dx <- expression(-(exp(a + b/x) * (b/x^2))) #resultado da primeira derivada D(dy.dx, 'x') #segunda derivada fx <- function(a,b,x) exp(a + b/x) * (b * (2 * x)/(x^2)^2) + exp(a + b/x) * (b/x^2) * (b/x^2) #resultado da segunda derivada optim(c(1), fx, a=3.3,b=-1.2,method="BFGS") # A ideia é encontrar o ponto de inflexao para o modelo representado pelo objeto y. Por definição o ponto de inflexao é o valor de x que faz a segunda derivada da # função igual a zero. Montei este rcmr porém não estou certo se a otimização de fx está igualando a zero, até porque os resultados não estão condizentes. # Gostaria se possível também de alguma dica de algum livro que tenha dicas de cálculos matemáticos com o R pois gostaria de adaptar algumas rotinas do maple # para o R, até porque recentemente precisei isolar um termo de um modelo não linear que manualmente não é tão simples e gostaria de fazer estes calculos com o # R e nao mais com o maple. Segue um exemplo simples de um calculo que gostaria de reproduzir com o R y = b0 + b1x b0=? b0 = y - b1x # esta é a saída que gostaria de encontrar com o R. Ate cheguei a ver um pouco sobre o pacote Ryacas mas sem muito sucesso. Ja deixo aqui meus agradecimentos ==================================== Samuel P. C. Carvalho Engº Florestal [UFLA] Mestre em Ciências Florestais [UFLA] Doutor em Recursos Florestais [ESALQ/USP] =============================================

3 2

[OT] Curso geo[r-]referenciamento
by Benilton Carvalho 12 Set '13

12 Set '13

Alguém saberia indicar um bom curso no tópico no MS? Ou mesmo à distância? Grato,

4 3

Obtendo valores de x a partir de valores determinados y em uma equação de segundo grau
by Alisson Lucrecio 12 Set '13

12 Set '13

Caro Colegas da lista r-br, Bom dia. Eu estou fazendo um trabalho de fertilidade do solo e preciso obter o níveis de potássio (K) do solo que correspondem a 90, 70 e 50 % da produtividade relativa (PR), eu fiz um ajuste para de uma equação do segundo grau (PR = -0.0042 K^2 + 1.3248 K ), assim para obter o K no solo posso fazer pelo delta, mas não queria resolver por delta para encontrar os valores de K da função, assim alguém saberia me dizer se existe outro foma? Obrigado. Alisson Lucrécio da Costa

2 1

Re: [R-br] Gerar e analizar dados por TRI
by Ivan Bezerra Allaman 12 Set '13

12 Set '13

A mensagem de erro é extremamente óbvia!!! e o seu exemplo não é reproduzível pois faltou colocar "library(ltm)". O problema está na especificação do argumento "model"!!! Ao abrir a função tem-se: probs <- if (model == "grm") { gammas <- lapply(thetas, function(x) { nx <- length(x) if (IRT) cbind(plogis(x[nx] * (z - matrix(x[-nx], n, nx - 1, TRUE))), 1) else cbind(plogis(matrix(x[-nx], n, nx - 1, TRUE) - x[nx] * z), 1) }) lapply(gammas, function(x) { nc <- ncol(x) cbind(x[, 1], x[, 2:nc] - x[, 1:(nc - 1)]) }) } Aqui, este objeto é calculado de modo que contém valores negativos em suas listas. Logo, e claramente, a função "sample" irá acusar erro, uma vez que, X <- matrix(0, n, p) for (j in 1:p) { for (i in 1:n) X[i, j] <- sample(ncatg[j], 1, prob = probs[[j]][i, ]) } Este objeto "probs" é justamente usado no argumento "prob" da função "sample" e como sabemos, não existe probabilidade negativa segundo os axiomas de Kolmogorov. Detectei o erro, porém não posso lhe dar a solução pois TRI não é algo que entendo!! (s,f,p) Allaman

3 6

Modelos não-linear misto (nlme), obtenção do intervalo de confiança
by Alisson Lucrecio 12 Set '13

12 Set '13

Caro Colegas r-br, Boa tarde, Eu estou usando a função "augPred" para obter as estimativas do modelo não-linear misto, entretanto gostaria de plotar com o intervalo de confiança de 95%. Você saberiam qual função posso usar para obter o intervalo de confiança? Em lm posso usar a função predict(...., interval="confidence"), mas nlme ela não funciona. Obrigado. Alisson Lucrécio da Costa

1 0

PGLS - Phylogenetic Generalized Least Squares
by marcelo claro de souza 11 Set '13

11 Set '13

Boa tarde, Gostaria de comparar a variação no acúmulo de N foliar entre duas comunidades vegetais utilizando regressões lineares filogenéticas (PGLS ou OLS). Tal análise foi realizada por Penũelas et al 2011 (tabela 2) ao comparar o acúmulo de nutrientes floiares entre duas comunidades vegetais, uma em Borneo e a outra no Hawaii ( http://www.esajournals.org/doi/abs/10.1890/ES10-00185.1) O pacote CAPER possui a função 'pgls' capaz de rodar esse tipo de análise, entretanto, eu só encontrei exemplos sobre comparações de duas variáveis dentro da mesma comunidade: require(caper) data(shorebird) shorebird <- comparative.data(shorebird.tree, shorebird.data, Species, vcv=TRUE) mod <- pgls(log(Egg.Mass) ~ log(M.Mass), shorebird) summary(mod) Alguém poderia me dar uma luz sobre como fazer isso, mas comparando duas populações? Agradeço imensamente por qualquer ajuda. Abraço Marcelo

1 0

Transformar string em binário
by Gerson R. Primo Jr 11 Set '13

11 Set '13

Pessoal, Estou com uma base que algumas variáveis são dicotômica (SIM/NÃO) e o campo está preenchido com a palavra SIM/NÃO. Gostaria de transforma-la em binárias, ou seja, SIM = 1 e NÃO = 0. É uma base de 800k e acredito que com essa mudança posso diminuir seu tamanho e também facilita a analise. T em algum comando que possa fazer essa transformação? No aguardo. Att, *Gerson R. Primo Jr** *

3 4