Re: [R-br] Parcelas sub-divididas com fatorial nas parcelas

Obrigado pelas dicas Walmes. E aproveito para fazer mais duas perguntas. Para aplicar a abordagem de modelos mistos para o exemplo que enviei o efeito aleatório deveria ser Parcela/Sub-parcela? Com os artifícios presentes no tutorial que você me indicou eu consigo ajustar o efeito de cada uma da

[R-br] Parcelas sub-divididas com fatorial nas parcelas

Boa noite pessoal, estou com experimento de três fatores e gostaria da opinião de vocês para modelar uma análise de variância. De forma mais direta eu tenho dois ambientes e em cada um deles tenho um experimento em arranjo fatorial em DIC com um tratamento adicional. #CMR set.seed(102) resp <-

Re: [R-br] Método de Jacobi

Da uma olhada neste site e verifique se os códigos te atendem. http://morotalab.org/Mrode2005/mme/mme.html#section0002 Att. Tiago. Wecsley Prates escreveu: Olá a todos... Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações line

Re: [R-br] Sistema é computacionalmente singular - Matriz/Mahalanobis

Boa tarde Andre, como o Walmes afirmou seu exemplo ainda não é reproduzível. Entretanto, segue um exemplo para servir de inspiração para o seu problema. Uma dica é você buscar entendimento da lógica de obtenção da matriz de variância e covariância Σ nas documentações da função mahalanobis.dist()

Re: [R-br] Dúvida sobre escrito em xlab

expression ("t ha"^-1) Att. Tiago. Cesar Rabak escreveu: Veja no help (e no demo do R) "expression" 2015-05-15 17:08 GMT-03:00 Leonardo Monteiro : > Senhores, boa noite. > > Como modifico o formato da legenda: t/ha para t ha-1? > Sendo que o -1 deve ficar sobreescrito. > > obrigado, > > Leon

Re: [R-br] Fatorial tiplo com dados adicionais

Te indico o pacote ExpDes se você quer uma rotina simples. Att. Tiago. Alisson Lucrécio escreveu: Caro Colegas da lista r-br, Preciso de uma ajudar para montar uma análise de variância para um fatorial tripo (Solo x Solução de Extração x Concentração) e um fator adicional água. Os dados fato

Re: [R-br] Algoritimo newtonRaphson

André tente este adaptação: newtonRaphson <- function (x) { x - ((exp(-x^2) - x^3 + 2) / (-2 * exp(-x^2) - 3 * x^2))} x <- 0.5# valor inicial old <- 0# valor anterior tol <- 0.0001y <- xi <- 0 while (abs(old - x) > tol) { i <- i + 1 old <- x x <- newtonRaphson(x) # Repetir se delta

Re: [R-br] Partição de somas de quadrado usando contrastes ortogonais

Walmes desconsidere a ultima mensagem. Não tinha me atentado para a falta de ortogonalidade entre os contrastes (somatório de ai*bi =0). Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciaçã

Re: [R-br] Partição de somas de quadrado usando contrastes ortogonais

Sim Walmes esta é a saída que desejo. No entanto, observe que se eu alterar a ordem dos contrastes as somas de quadrado alteram. Exemplo que você me enviou:mat.c [,1] [,2][1,] -12[2,]0 -1[3,] 1 -1 Df Sum Sq Mean Sq F value Pr(>F)trat 2 91.5

[R-br] Partição de somas de quadrado usando contrastes ortogonais

Boa tarde pessoal, estou querendo pré-determinar contrastes ortogonais na anova para estimar suas somas de quadrado. Entretanto, a ordem com que eu insiro os contrastes está diferenciando as somas de quadrado (acredito que seja função do ajuste sequencial) Alguém sabe como determinar as somas de

Re: [R-br] [Regressao logística] - comparacao de coeficientes beta

Acho que você deseja fazer a analise de identidade de modelos. O procedimento foi descrito por Regazzi (2004, 2010) para o delineamento inteiramente ao acaso e blocos ao acaso respectivamente. Caso tenha interesse nesta metodologia em R é só entrar em contato. Att. Tiago. Manuella Lech Cantu

Re: [R-br] Gráfico/CL50

Use a função points() Att. Tiago. Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA escreveu: Olá, Tenho um gráfico dose-resposta e preciso marcar um ponto específico, o ponto que representa a CL50. Como insiro essa informação no gráfico? Obrigada, Flávia - Mensagem original

Re: [R-br] Sub-matrizes em R

Obrigado izi mas creio que ainda não é o que procuro, pois, eu posso chegar a ter um número indefinido de sub-matrizes. Desta forma eu teria que declarar uma a uma. Ex: a, b, c ... k. Att. Tiago. # Tiago de Souza Marçal - Graduand

[R-br] Sub-matrizes em R

Boa noite, alguém sabe como selecionar 2 ou mais sub-matrizes de uma matriz maior de forma direta sem usar um loop CMR: mat <- matrix(1:25,nc=5) z = c(1,4)lz = c(2,5)k <- numeric(0) for(j in 1:2) print( k <- mat[z[j]:lz[j],z[j]:lz[j]] ) # Gostaria de chegar a este resultado sem o loop Estou

[R-br] for dentro de while

Boa noite, estou com um código que necessita o uso de um for dentro de um while. Entretanto, o uso do for aumentou o tempo de execução. Gostaria de saber se alguém tem alguma sugestão para contornar este problema Acredito que a questão seja organizacional, portanto, não estou enviando um CMR

Re: [R-br] Matrizes de variâncias e covariâncias do pacote nlme

Agradeço mais uma vez pelas valiosas recomendações Walmes. Vou verificar os pontos que você mencionou. Mas fiquei com uma ultima dúvida. Eu consegui reproduzir as estimativas de componentes de variância da função lme4::lmer() com a função nlme::lme() usando o comando list. Mas não consegui inser

Re: [R-br] Matrizes de variâncias e covariâncias do pacote nlme

Classes) help(varClasses) À disposição. Walmes. Em 10/11/2014 21:19, "Tiago Souza Marçal" escreveu: Boa noite pessoal, gostaria de saber quais as estruturas de matrizes eu posso obter combinando os argumentos correlation e weights nas funções lme e gls?? Att. Tiago.

[R-br] Matrizes de variâncias e covariâncias do pacote nlme

Boa noite pessoal, gostaria de saber quais as estruturas de matrizes eu posso obter combinando os argumentos correlation e weights nas funções lme e gls?? Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES

Re: [R-br] Ajuste de um semivariograma

por regressão não linear com “chutes” iniciais para os #parametros* *ff<-variogramLine(f,maxdist=2,n=500,min=1.0e-6)#geração de pontos do modelo* *plot(ff,col="blue")#plotando o modelo* *points(vgm1[,2],vgm1[,3], col="red")#plotando o variograma experimental* Em 7 de

[R-br] Ajuste de um semivariograma

Boa tarde pessoal, gostaria de saber como faço para manipular os pares de pontos no semivariograma ? Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e

Re: [R-br] Ordenamento entre colunas

Tenta deixar o cabeçalho da seguinte forma F09 D07 F05 F04 D02 F03 F02 D05 D06 F01 D03 F06 D04 D08 D01 F07 D09 D10 F08. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da

Re: [R-br] Ordenamento entre colunas

Uma maneira eficiente é usando o comando sort(). CMR: data = structure(c(-1.66349302275793, -0.403539151411636, 0.516649058010755, -0.865397435428028, -1.63106115127737, 0.447384559126393, 0.358753627652744, 1.37063744488393, -2.05266509010824, 0.275401420277465, 0.108863869878301, -2.1638459759

Re: [R-br] Multivariada no R

A bibliotecas 'candisc', 'MASS'. E es próprias funções da interface R (dos pacotes básicos) manova(), hclust(), dist(), heatmap() entre outras. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista d

Re: [R-br] ADAPTAR UMA MATRIZ NO R

Obrigado pelas considerações Walmes. A tal matriz que enviei é só parte de uma matriz maior que codifica efeitos dos cruzamentos, blocos, capacidade de combinação geral, capacidade de combinação especifica, efeito materno e efeito recíproco. A matriz de incidência para os efeitos recíprocos codi

[R-br] ADAPTAR UMA MATRIZ NO R

Boa noite pessoal,estou trabalhando com o desenvolvimento matricial dos efeitos dialélicos dos modelos de Griffing no R.No entanto, estou tendo dificuldade para adaptar a matriz CR (mostrada abaixo) para a matriz CR1 (mostrada abaixo) que codifica os efeitos recíprocos. Qualquer sugestão é bem v

Re: [R-br] Ajuda em análise de trilha (path analysis)

O pacote agricolae realiza a análise de trilha. Uma boa abordagem teórica você vai encontrar no livro "Modelos Biométricos Aplicados ao melhoramento" Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES

Re: [R-br] Teste de Tukey em experimento fatorial

Uma opção mais simples seria utilizar o package "ExpDes". # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ###

Re: [R-br] Histograma e boxplot juntos

Você pode fazer isso usando o comando par() para dividir a janela gráfica e depois usar as funções boxplot() e hist(). Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica

Re: [R-br] Parcela subdividida: problema na declaração do erro associado a parcela

Mesmo com o aviso use o comando summary(m0) e o quadro será formado. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ###

Re: [R-br] Plotar média de dois grupos

Talvez a função lattice::barchart() seja o que você esta procurando. De uma olhada no link: http://ridiculas.wordpress.com/2011/05/21/fatorial-com-desdobramento-da-interacao-grafico-e-tabela-de-medias/ Att. Tiago. # Tiago de

[R-br] Algoritmo de Fisher's Scoring para o procedimento REML

Bom dia pessoal, Estou trabalhando com a estimação de componentes de variância via modelos mistos usando o procedimento de máxima verossimilhança restrita (REML), no entanto, em minhas pesquisas pela internet encontrei o algoritmo de Fisher's Scoring para solução do seguinte modelo: y = Xa +

Re: [R-br] Gráfico de pontos

Nos envie um CMR para que possamos tentar ajudá-lo. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #

Re: [R-br] ScottKnott

Caso você queira conferir tente usar a função do pacote EXPDES.pt::dic(trata,respo,quali=c(TRUE),mcomp="sk"). Tente executar seu comando da seguinte forma: sk1<-with(dad,SK(y=respo, model="respo~trata",which="trata")) Att. Tiago. ##

Re: [R-br] Selecionar as variáveis físicas e químicas que melhor correlacionam com uma variável visual em experimento em fatorial 2x3

De uma olhada na análise de correspondência, ela esta implementada no pacote MASS. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento

[R-br] Customização de gráficos com a função plot() [Resolvido]

Consegui aqui. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #

[R-br] Customização de gráficos com a função plot()

Boa noite pessoal, gostaria de fazer um gráfico de dispersão com o CMR abaixo, entretanto, gostaria de variar o símbolo no parâmetro pch a cada tratamento e se possível a cor do símbolo. dados<-structure(list(trat = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 4L,

[R-br] Interface em R

Boa tarde pessoal, Andei dando uma olhada no pacote tclk e suas funcionalidades na criação de interfaces. No entanto, gostaria de saber se existe algum plugin que possa ser instalado no computador para que a interface pudesse ser construida via netbeans. Att. Tiago.

Re: [R-br] Substituir NA em data.frame 1 por dados data.frame 2

m data.frame 1 por dados data.frame 2 Segue recorte das data.frame 1 (normais.clima) e 2 (inmet.mensal). Obrigado. Alisson Lucrécio da Costa On Thursday, October 17, 2013 8:54 PM, Tiago Souza Marçal wrote: Envie um CM

Re: [R-br] Substituir NA em data.frame 1 por dados data.frame 2

Envie um CMR para podermos ajuda-lo. att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ###

[R-br] Remover NULL de um output da função lapply

Bom dia pessoal, gostaria de saber como remover os NULL do CMR abaixo. x=c("b","a","c") u=lapply(1:3,function(i){ if(x[i]=="a"){ 1} } ) att. Tiago. # Tiago de Souza Marçal - Graduando em Agro

Re: [R-br] Produção de Gráficos: PCA

Vinícius, já usei o pacote scatterplot3d para plotar gráficos de variáveis canônicas em 3D. Segue um CMR: dat <- data.frame(x=1:10, y=1:10, z=1:10, label=letters[1:10]) library(scatterplot3d) s3d <- scatterplot3d(dat$x, dat$y, dat$z) text(s3d$xyz.convert(dat$x, dat$y, dat$z), labels=dat$label,

[R-br] Subscrito em textos

Boa noite, gostaria de saber se ao escrever um texto em R é possível representar 2^ns na forma 2ns. CMR cat("\n 2^ns \n") Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de

Re: [R-br] Fatorial desbalanceado em bloco

Walmes, ao analisar um experimento similar me deparei com a seguinte situação a função popMeans() ajusta as médias pelo número de repetições existentes e na hora da obtenção das médias marginais ela usa a média aritmética, já o pacote TukeyC utiliza como estimador das médias marginais a média p

Re: [R-br] matematica com o R

Samuel, não utilizo a função optim portanto, não sei o procedimento que ela usa. Uma maneira seria as tentativas baseadas em aproximações gráficas embora ao se igualar a segunda derivada a zero a solução seria analítica (x = -b^2/2*b). Outra sugestão é você dar uma olhada na função while() e

Re: [R-br] Obtendo valores de x a partir de valores determinados y em uma equação de segundo grau

Alisson, se você não quer usar a solução algébrica você pode encontrar a solução por tentativas. Ex: fx=function(x){-0.0042*x^2 + 1.3248*x} x=seq(99,99.1,0.01) y=fx(x) plot(y~x) abline(h=90) abline(v=99.02) Att. Tiago. ###

Re: [R-br] ExpDes

60 628.20219 Erro em shapiro.test(anava$residuals) : sample size must be between 3 and 5000 Em 14 de agosto de 2013 14:37, Tiago Souza Marçal escreveu: Natalia ente nos códigos da função e edite esta parte retirando o teste de Shapiro wilk. Dep

Re: [R-br] ExpDes

Natalia ente nos códigos da função e edite esta parte retirando o teste de Shapiro wilk. Depois e só colar a função no R e utilizar a função que você tinha mostrado. fat3.rbd(hibrido, populacao, zona, bloco, Prod, sigT = 0.1, sigF = 0.1) # A propósito o alpha que você esta utilizando é de 10%

Re: [R-br] Função de busca em um arquivo de dados

-br] Função de busca em um arquivo de dados Vê se resolve. Sorteio = c(9, 3, 5) dados$caracu[Sorteio ] Em 13/08/2013 03:09, "Tiago Souza Marçal" escreveu: Prezados, estou com uma dificuldade de criar uma função de busca em um data.frame(). A situação é a seguinte e

[R-br] Função de busca em um arquivo de dados

Prezados, estou com uma dificuldade de criar uma função de busca em um data.frame(). A situação é a seguinte eu tenho um arquivo de dados com rótulos identificados em ordem crescente de 1 a n. carac id 1 A1 1 2 A2 2 3 A3 3 4 A4 4 5 A5 5 6 B2 6 7 B3 7 8

[R-br] Função str()

Prezados, andei dando uma olhada na saída do str() do exemplo da função aov() e fiquei com uma dúvida. Como posso definir a saída de dados para que eles apareçam de forma resumida no formato Named num como mostrado abaixo: List of 13 $ coefficients : Named num [1:13] 54.88 1.71 1.68 -1.82 -1

Re: [R-br] Problema com decimais

Sergio não estou conseguindo reproduzir o seu CMR. Mas caso você não tenha tentado da uma olhada na função round(). Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área

Re: [R-br] Popular um data.frame em dois loops

Tente usar o for(). Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas

Re: [R-br] R em pesquisa operacional

Rodolfo, da uma olhada na documentação do pacote "boot". No link abaixo você vai encontrar exemplos resolvidos pelo método simplex no pacote "boot".

Re: [R-br] Modelo não linear

x-1)),1,12, add=T,lty=2,col="red") ###Mas é um bocado de parâmetros para 12 pontos. Bem, espero ter ajudado. Em 30 de julho de 2013 19:39, Tiago Souza Marçal escreveu: Boa noite pessoal, e

[R-br] Modelo não linear

Boa noite pessoal, estou querendo modelar uma curva de embebição de sementes através do seguinte modelo não linear: y=A*(1-B*exp(-C*tempo)) + exp(D + E*(tempo-t0)) t0 = 6 Gostaria de saber se existe um procedimento parecido com o do pacote manipulate do Rstudio no R para ajustar este mode

[R-br] Eliminar o nome das colunas de um data.frame()

Boa noite pessoal, gostaria de saber se é possível dar nomes em colunas alternadas de um data.frame() ou apagar o nome das colunas. CMR: x<-c("a","a","a","a") y<-c(3,3,3,3) data.frame(x,y) x y 1 a 3 2 a 3 3 a 3 4 a 3 # Este é o output normal do R Gostaria de sa

[R-br] Matriz

Boa tarde pessoal estou tentando construir uma matriz com letras e números, entretanto, gostaria de retirar as aspas de ambos que aparecem na matriz. Alguém tem alguma sugestão? Segue abaixo um CMR: x<-c("a","a","a","a") y<-c(2,2,2,2) matrix(c(y,x),nc=2) Att. Tiago.

Re: [R-br] Fatorial duplo com tratamento adicional

Aparentemente não vejo nenhum problema mas envie um CMR que podemos tentar ajuda-la. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento

[R-br] Arranjo com 3 fatores

Bom dia pessoal, estou tentando montar um arranjo com os fatores inoculo, genótipo e tempo. Acho que não posso analisar por uma análise de fatorial triplo pelo fato das avaliações serem realizadas tempos diferentes. Alguma sugestão? Agradeço desde já pela ajuda. Att. Tiago. ##

Re: [R-br] erro bar

Geovane o ultimo valor do seu arquivo (63) só tem um representante, portanto, não é possível calcular o erro padrão ep=sd/sqrt(n). Tente utilizar o algoritmo abaixo: dados=structure(list(idade = c(35L, 35L, 35L, 36L, 36L, 36L, 36L, 37L, 37L, 40L, 40L), peso = c(59L, 64L, 87L, 52L, 90L, 45L, 87

Re: [R-br] Plotagem 3D

gronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas # From: tito.co...@gmail.com Date: Mon, 3 Jun 2013 10:56:33 -0300 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Plotagem 3D da uma olhada na função cloud do lattice Tito Conte 20

[R-br] Plotagem 3D

Bom dia pessoal, estou querendo plotar um gráfico tridimensional com os escores canônicos que estão abaixo: dados<-structure(list(V1 = c(-18.7175, -20.1915, -21.1557, -17.5965, -18.142, -13.857, -18.0086, -14.0436, -12.571), V2 = c(38.3004, 36.4955, 34.4964, 39.9303, 40.4819, 38.3763, 34.986

Re: [R-br] Anova de fatorial duplo em DBC (um fator quantitativo)

Caso continue dando problemas na função de um dput()nos seus dados para que possamos usar como CMR. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genétic

Re: [R-br] codigo de leitura de dados

Será que o erro pode estar associado com o operador % acompanhado dos números na terceira coluna. Att. Tiago. # Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética

Re: [R-br] Aumentar letras do plot glm

labels, mainand sub-titles, respectively De qualquer forma, obrigada. Renata 2013/5/27 Tiago Souza Marçal Normalmente quando eu trabalho com o plot()antes eu chamo a função par()e dentro dela eu uso o comando cex (ex: cex=2). Att. Tiago. Date: Sun, 26 May 2013 18:40:41 -0300 From

Re: [R-br] Aumentar letras do plot glm

Normalmente quando eu trabalho com o plot()antes eu chamo a função par()e dentro dela eu uso o comando cex (ex: cex=2). Att. Tiago. Date: Sun, 26 May 2013 18:40:41 -0300 From: renataecol...@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Aumentar letras do plot glm Boa noite a todos,

Re: [R-br] Anova de fatorial duplo em DBC (um fator quantitativo)

Você não enviou um CMR mas desconfio que a ordem de entrada dos seus dados esta incorreta. Tente usar o seguinte comando: fat2.dbc(fator1,fator2,bloco,resposta,quali=c(T,F),mcomp="tukey", fac.names=c("fator1","fator2"), sigT=0.05,sigF=0.05) No exemplo acima eu considerei que o seu fator 2 é q

Re: [R-br] Classificar fosforo e textura do solo

Hélio tente o seguinte comandos: recP<-function(argila,Psolo){ t<-numeric(0) if(100>=argila && argila>=60 && 2.7>=Psolo && Psolo>=0){ t<-c(t,"Multo Baixo") } if(100>=argila && argila>=60 && 5.4>=Psolo && Psolo>2.7){ t<-c(t,"Baixo") } if(100>=argila && argila>=60 && 8>=Psolo && Psolo>5.

Re: [R-br] transformação raiz quadrada do arco-seno

Se os seus valores estiverem em % use o seguinte comando: y <- asin(sqrt(x/100))*180/pi Att. Tiago. Date: Tue, 14 May 2013 14:52:22 -0700 From: isabelsousaamo...@yahoo.com.br To: R-br@listas.c3sl.ufpr.br Subject: [R-br] transformação raiz quadrada do arco-seno Boa noite pessoal, Gostaria de faze

Re: [R-br] Matriz X (Y = XB + e) [Resolvido]

Obrigado a todos! Att. Tiago. CC: r-br@listas.c3sl.ufpr.br From: edinhoes...@yahoo.com.br Date: Thu, 25 Apr 2013 12:00:17 -0400 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Matriz X (Y = XB + e) ?solve [. ]'s.Edson LiraEstatísticoMa-Am Em 23/04/2013, às 15:19, Tiago Souza Marçal esc

[R-br] Matrix X (Y = XB + e)

Prezados, gostaria de saber se existe uma função que simplica a obtenção desta matriz? m<-matrix(c(rep(1,15),rep(1,5),rep(0,10),rep(0,5),rep(1,5),rep(0,5),rep(0,10),rep(1,5)),nr=15) Att. Tiago. ___ R-br

[R-br] Matriz X (Y = XB + e)

Prezados, gostaria de saber se existe uma função que simplica a obtensão desta matriz? m<-matrix(c(rep(1,15),rep(1,5),rep(0,10),rep(0,5),rep(1,5),rep(0,5),rep(0,10),rep(1,5)),nr=15) Att. Tiago. ___ R-br maili

Re: [R-br] Matrix Hessiana - duvida

nao deveria ter retornado um array 3x3x3? 2013/4/13 Tiago Souza Marçal Carlos tente os seguintes comandos: d3 <- deriv3(~x1/((1+((x1/x2)**(x3*t)-1))),c("x1","x2","x3"), function(t,x1,x2,x3){ NULL }) t = rnorm(10,5,3) x1=8500 x2=0.0156 x3=0.26 hess

Re: [R-br] Erro no calculo da matriz distancia

O seu exemplo não é reproduzivel. Envie um CMR. Att. Tiago. Date: Sat, 13 Apr 2013 15:10:20 +0100 From: bernardosss...@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Erro no calculo da matriz distancia Olá a todos Estou a tentar calcular a distancia euclidiana e da me o seguinte erro di

Re: [R-br] Matrix Hessiana - duvida

Carlos tente os seguintes comandos: d3 <- deriv3(~x1/((1+((x1/x2)**(x3*t)-1))),c("x1","x2","x3"), function(t,x1,x2,x3){ NULL }) t = rnorm(10,5,3) x1=8500 x2=0.0156 x3=0.26 hessiana<-d3(t,x1,x2,x3) Espero ter contribuido Att. Tiago. From: ccpsilv...@hotmail.com To: r-br@listas.c3sl.ufpr.br

Re: [R-br] duvidas sobre o r

Para fazer uma anova em dql é só usar os seguntes comantos: trat<-factor(trat) linha<-factor(linha) coluna<-factor(coluna) av<-aov(resp~linha+coluna+trat) # Como pode ver a linha e a coluna são somadas aos tratamentos. Caso você queira fazer uma analise completa de um experimento dql eu te indic

Re: [R-br] (sem assunto)

Obrigado Henrique Att. Tiago. From: www...@gmail.com Date: Wed, 20 Mar 2013 20:57:26 -0300 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] (sem assunto) Tente assim: plot(0, main=expression(KNO[3]~L[-1])) 2013/3/20 Tiago Souza Marçal Boa noite pessoal, gostaria de saber como

Re: [R-br] (sem assunto)

Obrigado Fernando. Att. Tiago. > Date: Wed, 20 Mar 2013 21:00:58 -0300 > From: fernandoma...@gmail.com > To: r-br@listas.c3sl.ufpr.br > Subject: Re: [R-br] (sem assunto) > > Ou > > plot(0, main=expression(KNO[3] * phantom(x) * L[-1])) > > > > > --- > Fernando de Pol Mayer > Doutorando em Est

Re: [R-br] (sem assunto)

Obrigado Walmes. Att. Tiago Date: Wed, 20 Mar 2013 20:58:04 -0300 From: walmeszevi...@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] (sem assunto) Espaços nas expressões são obtidas com o tio (~), veja plot(0, main=expression(KNO[3]~~~L[-1])) À disposição. Walmes.

[R-br] (sem assunto)

Boa noite pessoal, gostaria de saber como substituir o sinal de vezes por um espaço? CMR plot(0, main=expression(KNO[3]%*%L[-1])) Att. Tiago. ___ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.u

Re: [R-br] Pacote drc

LSE,type="average",xlim=c(0.15,5),ylim=c(0,1),log="x") axis(1,conc) y=seq(0,1,0.1) axis(2,y) lines(tratVec, Li, lty = 2) lines(tratVec, Ls, lty = 2) Stephan

Re: [R-br] Pacote drc

,xlim=c(0.15,5),ylim=c(0,1),log="x") axis(1,conc) y=seq(0,1,0.1) axis(2,y) lines(tratVec, Li, lty = 2) lines(tratVec, Ls, lty = 2) Stephan

Re: [R-br] Pacote drc

er nenhuma correção como proposto por Abbott. Além disso é preciso verificar qual tipo de distribuição (e parâmetros) que você está usando, pois existe uma grande diversidade no drc quando comparado por exemplo com o Polo. Att,

[R-br] Pacote drc

Bom dia pessoal, gostaria de saber se alguem da lista utiliza o pacote drc para estimar a CL-50. Caso tenha gostaria que me esclarecesse as seguintes dúvidas: Ele calcula a mortalidade corrigida pelo critério de Abbott? Porque as estimativas geradas por este pacote são diferentes das geradas

Re: [R-br] ANOVA PARA MODELOS NÃO LINEARES

Ok Walmes, entendi perfeitamente quando o modelo é de intercepto. Agradeço mais uma vez pela ajuda. Att. Tiago. Date: Thu, 7 Mar 2013 16:37:54 -0300 From: walmeszevi...@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] ANOVA PARA MODELOS NÃO LINEARES Sim, a lógica é essa. As vezes i

Re: [R-br] ANOVA PARA MODELOS NÃO LINEARES

Ok Walmes, muito obrigado. Então sempre que eu tiver um modelo e aplicar estas restrições, ou seja, arbitrar que um certo numero de constantes é 0 e observar uma redução do modelo para uma constante ou para um produto de uma constante por uma variável como no seu exemplo (A*t), eu terei um m

[R-br] ANOVA PARA MODELOS NÃO LINEARES

Boa tarde pessoal, eu dei uma olhada no post que o Walmes fez no dicas ridículas usando o modelo y = A*t/(V+t) + D*t e ele justificou que se V e D fossem zero este se tornaria um modelo de intercepto, ou seja, y = A*t . Entretanto, estou trabalhando com um modelo logístico y = A/(1+B*exp(-C*x))

Re: [R-br] 'subsetting' de um dataframe, como fazer uma função ou loop?

No exemplo que você enviou, o objetivo seria calcular a media para todas as variáveis numéricas em função de year? Caso positivo tente usar o for() ou lapply(). Caso tenha problemas envie um CMR para tentarmos ajuda-lo. Caso preferir use um dput() ao enviar o CMR. Att. Tiago. Date: Thu, 7 Mar 2

Re: [R-br] Quadrado médio de duas características QM(y1+y2)

t;- anova(lm(resp1 ~ genotipo, data = dados)) # SQ da resposta 1 ANOVA_resp2 <- anova(lm(resp2 ~ genotipo, data = dados)) # SQ da resposta 2 ANOVA_soma <- anova(lm(I(resp1 + resp2) ~ genotipo, data = dados)) # SQ da soma Espero ter ajudado. att,FH 2013/2/27 Tiago Souza Marçal Fernan

Re: [R-br] PLOT

Oi Cezar, eu dei uma alterada nos seus comandos. Da uma olhada se isso resolve: plot(Deslocamento,direita,type="b",pch=2,ylim=c(500,6000),axes=F) lines(Deslocamento,esquerda,type="b",pch=17,col="black") axis(1,at=c(20,50,70,100,300,500,700,1000)) axis(2,at=seq(500,6000,500)) Att. Tiago From: cez

Re: [R-br] Quadrado médio de duas características QM(y1+y2)

apenas um chute. att,FH 2013/2/27 Tiago Souza Marçal Boa noite pessoal, estou querendo estimar a correlação genética entre duas variáveis (y1 e y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM (a+b). Alguém sabe como estimar este quadrado médio para as duas variá

[R-br] Quadrado médio de duas características QM(y1+y2)

Boa noite pessoal, estou querendo estimar a correlação genética entre duas variáveis (y1 e y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM (a+b). Alguém sabe como estimar este quadrado médio para as duas variáveis no R. Att. Tiago.

Re: [R-br] GRANDES RELATÓRIOS NO R

então, ver um-a-um. Você pode também imprimir o objeto para um arquivo de texto e abri-lo num bloco de notas. Se não for esse o caso, talvez você possa reduzir o tamanho da font de texto que é mostrada na tela. Outros usuários podem ter ideias melhores. Daniel 2013/2/23 Tiago Souza Marçal

[R-br] GRANDES RELATÓRIOS NO R

Bom dia pessoal, eu fiz uma analise no r que gerou um relatório muito grande e eu não estou conseguindo todas as analises. Como faço para ver todas as análises? Att. Tiago.___ R-br mailing list R-br@listas.c3

[R-br] Contrastes multivariados

Boa noite pessoal, estou utilizando a MANOVA para analisar um experimento fatorial com duas características: f1<-c(rep("A",20),rep("B",20)) f2<-rep(c(gl(4,5)),2) res1<-c(0,0, 20, 0, 20, 20, 0, 0, 20, 40, 30, 40, 40, 20, 40, 30, 20, 30, 30, 50, 40, 40, 40, 70, 60, 40, 40, 40, 60, 20, 30, 30

Re: [R-br] Método de Newton Rhapson...

ntar o método de Gauss newton e os resultados baterem com a saida do R, então, funcionou o metodo de gauss. Att. André Em 08/02/2013 23:35, Tiago Souza Marçal < tiagosouzamar...@hotmail.com > escreveu: Só mais uma duvida a equação de densidade de probabilidade possu

Re: [R-br] Método de Newton Rhapson...

o metodo de gauss. Att. André Em 08/02/2013 23:35, Tiago Souza Marçal < tiagosouzamar...@hotmail.com > escreveu: Só mais uma duvida a equação de densidade de probabilidade possui parâmetros populacionais porque quando são feitas as derivadas parciais estes parâmetros se

Re: [R-br] Método de Newton Rhapson...

liado 43 anos seguidos, cujos valores são referentes ao indice máximo pluviometrico. Os anos foram de 1956 à 1999. Att. André Em 08/02/2013 21:01, Tiago Souza Marçal < tiagosouzamar...@hotmail.com > escreveu: O método de Gaus Newton é um método de busca numérica muito utilizado

Re: [R-br] Método de Newton Rhapson...

(GVE), bem como suas derivadas parcais em relação a cada parâmetro. https://www.transferbigfiles.com/f1220ccb-8cb3-46c6-8a68-3c95fc24ebc4?rid=CnKlj7VjIyqbgcMXXXqFfg2 desde já agradeço! Att. André Em 08/02/2013 07:01, Tiago Souza Marçal < tiagosouzamar...@hotmail.com > escreveu: An

Re: [R-br] Método de Newton Rhapson...

André se você quiser a resolução especificamente por este método da uma olhado no dicas ridículas que o Walmes postou algo a respeito. Mas eu estou trabalhando com o método de GAUSS-NEWTON e vou te passar o script para caso você queira utiliza-lo. No site abaixo tem um material explicando o mét

  1   2   >