Setembro 2013 - R-br

Manipulação de base.
by Gerson R. Primo Jr 18 Set '13

18 Set '13

Pessoal, Estou com a seguinte situação: Tenho uma base e quero separa-las em duas amostras uma com 20% da base e outra com 80%. Só que na base com os 20% quero "balancear" ponderando por uma variável dicotônica, onde eu tenha 50% de 'sim' e 50% de 'não' dessa variável. Alguém sabe como posso fazer isso no R? *Att, Gerson R. Primo Jr** *

2 2

Fatorial desbalanceado em bloco
by Ze Henrique 16 Set '13

16 Set '13

Gostaria de saber se algum pode me ajudar e se o erre faz isto: Como obter as médias marginais de um experimento fatorial em bloco desbalanceado, porem quando não der interação. Obrigado desde já. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

3 4

Re: [R-br] Gráfico (Intervalo com números decimais)
by luciano.f.sgarbi＠gmail.com 14 Set '13

14 Set '13

Use ylim=c(0,0.6) ---------- Enviada do meu celular Nokia ------Mensagem original------ De: Eric Leandro <eric_lsp(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Date: Sexta-feira, 13 de Setembro de 2013 20h9min50s GMT-0700 Subject: [R-br] Gráfico (Intervalo com números decimais) 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. _______________________________________________ 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.

2 1

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 13 Set '13

13 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 13 Set '13

13 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