Re: [R] Selecting elements

2021-08-25 Thread Silvano Cesar da Costa
Wow, That's exactly what I want. But, if possible, that a list was created with the selected elements (variable and value). Is it possible to add in the output file? Thank you very much. Prof. Dr. Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas Departamento de

Re: [R] Selecting elements

2021-08-23 Thread Silvano Cesar da Costa
, decreasing=TRUE), ]) # The end result is: (Selected.data.2 = Order[c(1,2,3,4,5,6,7,8,9,11), ]) # Scheme: 3A - 4B - 2C - 1D sort(Selected.data.2$Var.1) How to make the selection of the 10 elements automatically? Thank you very much. Prof. Dr. Silvano Cesar da Costa Universidade Estadual de

Re: [R] Selecting elements

2021-08-20 Thread Silvano Cesar da Costa
ons using R? Thanks, Prof. Dr. Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas Departamento de Estatística Fone: (43) 3371-4346 Em sex., 20 de ago. de 2021 às 03:28, PIKAL Petr escreveu: > Hallo > > I am confused, maybe others know what do you wa

[R] Selecting elements

2021-08-19 Thread Silvano Cesar da Costa
Dr. Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas Departamento de Estatística Fone: (43) 3371-4346 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see ht

[R] Common elements

2019-03-08 Thread Silvano Cesar da Costa
,] 1186389 [25,] 1190425 [26,] 1192235 How can I do this? Thanks. Prof. Dr. Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas Departamento de Estatística Fone: (43) 3371-4346 [[alternative HTML version deleted

[R] Reading access file

2015-05-14 Thread silvano
Hello everybody. I have a access file to read in R but I can’t to do this. I used Hmisc package, but it doesn’t work. Someone has the commands to read this kind of file? I attached the access file. Thanks. Silvano. __ R-help@r-project.org mailing

[R] Differences between MTDFReml and kinship2

2014-11-05 Thread silvano
variance = 0.767187 residual variance = 5.6712 Both of them use REML method. Could someone tell why the difference? Thanks a lot, Silvano. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] Problems with table

2014-01-14 Thread silvano
, format.args=list(big.mark = ., decimal.mark = ,), caption.placement='top', table.placement='H') cat(\n\n\n) } Infra - function(){ for (i in 8:10){ x - paste(Q, i, sep=) cria.tabela(x, i) } } Infra() How can I solve this problem? Thanks, Silvano. --- Este email está limpo de

[R] Function inside Sweave

2014-01-08 Thread Silvano Cesar da Costa
Veterinária/Dados/A2013/A2013.rec') head(dados) tail(dados) str(dados) attach(dados) serie_1 = subset(dados, SERIE=='1') head(serie_1) tail(serie_1) attach(serie_1) How can I do this? Thanks a lot, - Silvano Cesar da Costa Universidade Estadual de

[R] ftable and data.frame

2013-12-20 Thread silvano
logistic regression and I have many tables. Thanks, Silvano. --- Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Problems with xtable?

2013-12-13 Thread Silvano Cesar da Costa
{tabular} \end{table} and I don't want this output Margins computed over dimensions in the following order: 1: DISCIPLINA 2: Q1 How can I get out it? Thanks, - Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas

[R] Randomization

2013-08-23 Thread Silvano
Hi, I have a set of 80 animals and their respective weights. I would like create 4 groups of 20 animals so that the groups have means and variances with values ??very close. How can I make this randomization in R? Thanks, -- Silvano Cesar da Costa

[R] timereg

2013-08-05 Thread Silvano
0.46100 0.004 part of them... Thanks a lot, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https

[R] means in tables

2013-04-10 Thread Silvano Cesar da Costa
(tab1)){ for(j in 1:ncol(tab1)){ for(k in 1:nrow(tab2)){ for(l in 1:ncol(tab2)){ Médias = tab1$i[j] Médias but it does't work. I don't know programming. How can I do this? Thanks, - Silvano Cesar da Costa Universidade

[R] SNPRelate problems

2012-12-18 Thread Silvano Cesar da Costa
Hi, I'm trying convert plink files to gds (SNPRelate). I have the files: baep.fam, baep.bed and baep.bim and the program to convert data is: setwd('C:/Silvano/Incor/Baependi/Dados Baependi/') library(gdsfmt)# versão 0.9.10 library(SNPRelate) # versão 0.9.8 # Arquivos PLINK BED

[R] Kinship2 and GenABEL

2012-11-09 Thread Silvano Cesar da Costa
= polygenic(altura ~ SEX + idade, data=Dados, kin=kmat)) Erro em intI(i, n = d[1], dn[[1]], give.dn = FALSE) : invalid character indexing How can I calculate heritabilty using GenABEL with kinship matrix of kinship2 package? Thanks, - Silvano Cesar da

[R] Extracting columns

2012-11-08 Thread Silvano Cesar da Costa
. Is there a way to schedule these operations? - Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas Departamento de Estatística Fone: (43) 3371-4346 __ R-help@r-project.org mailing list

[R] Interaction

2012-10-29 Thread Silvano Cesar da Costa
Hi, I'm fitting a model with 3 variables: A, B and SNP. The response variable is Y. I would like fit the following model, in this order: Y ~ A + B + A*B + SNP In general, the R sets of the form: Y ~ A + B + SNP + A*B How do this? - Silvano Cesar

Re: [R] factor or character

2012-10-24 Thread Silvano Cesar da Costa
Hi, still doesn't work. Hello, I've just seen the error, you are _not_ searching for colnames in mod5.sig$snps. Corrected: Selec = todos[ , colnames(todos) %in% mod5.sig$snps] Hope this helps, Rui Barradas Em 23-10-2012 21:17, Silvano Cesar da Costa escreveu: Hi Rui, it doesn't

[R] factor or character

2012-10-23 Thread Silvano Cesar da Costa
... But it doesn't work. The output is: Selec = todos[ , colnames(todos) %in% snps] head(Selec) data frame with 0 columns and 6 rows What's is wrong? Thanks a lot, - Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas

Re: [R] factor or character

2012-10-23 Thread Silvano Cesar da Costa
, like str shows). Hope this helps, Rui Barradas Em 23-10-2012 19:02, Silvano Cesar da Costa escreveu: Hi, The program below work very well. (snps = c('rs621782_G', 'rs8087639_G', 'rs8094221_T', 'rs7227515_A', 'rs537202_C')) Selec = todos[ , colnames(todos) %in% snps] head(Selec

[R] Reading multiple files

2012-09-26 Thread Silvano Cesar da Costa
=~/Silvano/Arq, pattern=Dados+.*raw) names = substr(filenames, 1, 7) for(i in names){ filepath = file.path(~/Silvano/Dados, paste(i, .raw, sep=)) assign(i, read.delim(filepath, colClasses=c(rep(character, 5), rep(numeric, 5)), sep = )) } It happens

Re: [R] Reading multiple files

2012-09-26 Thread Silvano Cesar da Costa
Hi, I didn't notice problems with this. Will your data be read in correctly if you do away with the colClasses argument to read.delim (or read.table)? Jean Silvano Cesar da Costa silv...@uel.br wrote on 09/26/2012 09:11:33 AM: Hi, I have 35 data files for reading. I would like get

Re: [R] Reading multiple files

2012-09-26 Thread Silvano Cesar da Costa
the following code should work. for(i in names){ filepath = file.path(~/Silvano/Dados, paste(i, .raw, sep=)) assign(i, read.table(filepath)) } If not, then my guess is the error (what is the error message???) is due to the file referencing or naming. For example, it's not clear

Re: [R] Common elements in columns

2012-09-03 Thread Silvano Cesar da Costa
-data.frame(C1=sample(LETTERS[1:10],10,replace=FALSE),B1=rnorm(10,3)) set.seed(5) df4-data.frame(C1=sample(LETTERS[1:15],10,replace=FALSE),A2=rnorm(10,15)) df1$C1[df1$C1%in%df2$C1[df2$C1%in%df3$C1[df3$C1%in%df4$C1]]] #[1] G E H J A.K. - Original Message - From: Silvano Cesar da Costa silv

Re: [R] Common elements in columns

2012-09-03 Thread Silvano Cesar da Costa
) set.seed(3) df3-data.frame(C1=sample(LETTERS[1:10],10,replace=FALSE),B1=rnorm(10,3)) set.seed(5) df4-data.frame(C1=sample(LETTERS[1:15],10,replace=FALSE),A2=rnorm(10,15)) x - list(df1$C1, df2$C1, df3$C1, df4$C1) Reduce(intersect, x) [1] G E H J On Mon, Sep 3, 2012 at 7:01 AM, Silvano

[R] Common elements in columns

2012-09-02 Thread Silvano Cesar da Costa
check which common elements in the columns of C1 4 files? Thanks, - Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas Departamento de Estatística Fone: (43) 3371-4346 __ R

[R] Pedigrees

2012-08-16 Thread Silvano Cesar da Costa
, famid, .pdf)) plot(pedig[famid]) dev.off() } Error em `[.pedigreeList`(pedig, famid) : Familiy famid not found but is not working. How to solve the problem? Thanks, - Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências

[R] NULL column

2012-08-08 Thread Silvano Cesar da Costa
1 rs46 163.3730 3.304276 3.279925e-06 1 rs63 162.7924 2.084678 5.023893e-06 1 rs24 162.4252 1.837208 5.509042e-06 1 How can I do this? Thanks, - Silvano Cesar da Costa Universidade Estadual de

Re: [R] NULL column

2012-08-08 Thread Silvano Cesar da Costa
Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of arun Sent: Wednesday, August 08, 2012 3:06 PM To: Silvano Cesar da Costa Cc: R help Subject: Re: [R] NULL column HI, Try this: dat1-read.table(text=     beta0  beta1

[R] regression analysis

2012-07-25 Thread Silvano Cesar da Costa
= numeric(1) for(i in 1:1){ mod = lm(RESP~x[i] p-value[i] = summary(mod)$coe[2,4] } Could someone suggest a reading material or any suggestions, I thank you. - Silvano Cesar da Costa Universidade Estadual de Londrina Centro de Ciências Exatas

[R] France Model

2012-05-04 Thread Silvano
faltante ou infinito quando avaliando o modelo Somebody knows some package? Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org

[R] Challenge

2012-02-22 Thread Silvano
this situation? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] Longitudinal data

2011-12-27 Thread Silvano
. What would be the equivalent in R? Thank you. -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Longitudinal data

2011-12-27 Thread Silvano
Hi Uwe, was a great suggestion. Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: Uwe Ligges lig...@statistik.tu

[R] Notation

2011-11-29 Thread Silvano
Hi, what's mean / in command: betareg(inf~Grupo/Sexo, data=dados) it's a effect nested? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r

[R] Problema with Excel files

2011-11-07 Thread Silvano
, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] Fit Gompertz' curve'

2011-08-16 Thread Silvano
(Tempo, Trat, Valor, ylim=c(0, 2), las=1, lty=c(1,2,3,4), lwd=3, bty='l', col=c('red','blue','magenta','green'), ylab=Média de Gases, xlab=Tempo (h), trace.label=Doses) Thanks, -- Silvano Cesar da

[R] Using Function

2011-08-02 Thread Silvano
), desc(CVLDL), desc(TRIG), desc(URICO), desc(SAQRS), desc(SOKOLOW_LYON), desc(CORNELL), desc(QRS_dur), desc(Interv_QT) ) Thanks a lot, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346

Re: [R] Tables and merge

2011-07-05 Thread Silvano
- Original Message - From: Silvano silv...@uel.br To: r-help@r-project.org Sent: Thursday, June 30, 2011 9:07 AM Subject: Tables and merge Hi, I have 21 files which is common variable CODE. Each file refers to a question. I would like to join the 21 files into one, to construct

[R] Pie chart

2011-05-19 Thread Silvano
'), clockwise=F, radius=.7) Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Boxplot in order

2011-05-05 Thread Silvano
this? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] Hartley's table

2011-02-16 Thread Silvano
14.31403 15.04529 15.73600 16.39159 1415 2 865.53893 948.24795 3 82.18823 83.40749 4 58.02112 61.10224 5 19.09581 34.64536 6 22.58582 23.46946 7 17.01639 17.61392 What's wrong? Thanks. -- Silvano Cesar da Costa Departamento de

Re: [R] HLM Model

2011-01-28 Thread Silvano
+school, test, correlation=corSymm(form = ~ 1 |id), weights=varIdent(form = ~ 1|team), method=REML) summary(unstruct) -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346

[R] subset factor?

2011-01-14 Thread Silvano
, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] subset

2010-12-15 Thread Silvano
63.60 C 33 med3 52.17 C 34 med3 75.10 C 35 med3 75.38 C 36 med3 55.78 -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371

Re: [R] subset

2010-12-15 Thread Silvano
Phil, this is exactly what do I want. Thanks a lot. -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: Phil Spector spec

[R] Tukey's table

2010-11-03 Thread Silvano
? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] step

2010-08-27 Thread Silvano
Hi, how can I change the significance level in test F to select variable in step command? I used step(model0, ~x1+x2+x3+x4, direction=c(forward), test='F', alpha=.05) but it does't work. -- Silvano Cesar da Costa Departamento de Estatística

Re: [R] Games

2010-08-16 Thread Silvano
Hi, I want to thank all the suggestions sent, especially that of Hans. Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From

[R] Games

2010-08-13 Thread Silvano
? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Games

2010-08-13 Thread Silvano
, round8, round9, round10, round11 )) -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: peter dalgaard pda...@gmail.com To: Michael

[R] xtable

2010-07-21 Thread Silvano
Hi, How do I build a table from a regression model adjusted using xtable? Commands are: modelo1 = lm(Y~X1 + X2) influencia = influence.measures(modelo1) require(xtable) xtable(influencia) but it isn't work. Thanks, -- Silvano Cesar da Costa

[R] xtable with Sweave

2010-06-14 Thread Silvano
') @ I'm getting the following message: Too many unprocessed floats in Latex file. How to avoid these messages appearing? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346

Re: [R] Table and Sweave

2010-05-12 Thread Silvano
on. -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: Ista Zahn istaz...@gmail.com To: r-help@r-project.org Cc: Silvano silv...@uel.br Sent: Tuesday, May 11, 2010 5:48

[R] Table and Sweave

2010-05-11 Thread Silvano
, results=verbatim= Idade_freq @ or = xtable(Idade_freq) @ but they are not the way I want. Any suggestions? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R

[R] Extracting a subset

2010-05-10 Thread Silvano
subset but it did not work. Any suggestions? Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] Extracting a subset

2010-05-10 Thread Silvano
-- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: Ivan Calandra ivan.calan...@uni-hamburg.de To: r-help@r-project.org Sent: Monday, May 10, 2010 11

Re: [R] Latex and Stangle()

2010-05-06 Thread Silvano
It worked very well. Thanks Frank and Ruihong, both solutions were great. -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From

[R] Latex and Stangle()

2010-05-05 Thread Silvano
Hi, I'm using the Sweave and I would like include codes of the R in my LaTeX file. I extracts the R code with Stangle (), whose name is Relatorio.R but I can't include it in the Latex file as an appendix. Suggests? Thanks, -- Silvano Cesar da Costa

[R] About logistic regression

2010-04-01 Thread Silvano
1.2980.194 genero -0.355358 0.402650 -0.8830.377 grau:genero 0.007052 0.580837 0.0120.990 What am I doing wrong? Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346

[R] Mosaic

2010-03-24 Thread Silvano
? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] Hartley's table

2010-02-26 Thread Silvano
Hi, Does anyone know how to generate Hartley's table in R? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https

[R] Normality in split-plot design

2010-02-22 Thread Silvano
? Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] Split Plot and Tukey

2010-02-17 Thread Silvano
)) names(anava) summary(anava) require(agricolae) HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05) Thanks -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help

[R] About scales in graphics

2010-02-08 Thread Silvano
, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] About graphics

2010-02-05 Thread Silvano
? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] About graphics

2010-02-05 Thread Silvano
It was exactly what I needed, thank you. -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: Duncan Murdoch murd...@stats.uwo.ca

[R] Bartlett' test

2009-11-30 Thread Silvano
, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Tests in Anova

2009-11-25 Thread Silvano
Sarah and Dennis, thank you for help and page reference. Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https

[R] Tests in Anova

2009-11-24 Thread Silvano
Hi, how can I make tests like Dunnett and Duncan using R? Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https

Re: [R] boxplot, data frame and html

2009-10-29 Thread Silvano
Gregoire, it worked very well, thanks -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: Gregoire Pau gregoire@embl.de To: r

[R] boxplot, data frame and html

2009-10-28 Thread Silvano
, br=TRUE) closePage(p) but isn't working. What's wrong? Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https