Re: [R] Selecting elements

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 Departamen

Re: [R] Selecting elements

, 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

options 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 y

[R] Selecting elements

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

,] 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

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

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

j] tab1.txt = xtable(tab1.prop, align="l|r", label=Questao, caption=paste(capt)) print(tab1.txt, format.args=list(big.mark = ".", decimal.mark = ","), caption.placement='top', table.placement='H') cat("\n\n\n") } Infra &l

[R] Function inside Sweave

p="") cria.tabela(x) } } geral() I need to change the variable DISCIPLINA created in the function by the variable DISCIPLINA in my dataset, but not working. Reading data was taken with: require(foreign) require(xtable) dados = read.epiinfo('C:/Colegiado de Veterinária/Dados

[R] ftable and data.frame

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?

& 24,07 & 100,00 \\ 6HIT011 & 0,00 & 0,00 & 3,70 & 27,78 & 68,52 & 100,00 \\ 6MOR012 & 0,00 & 0,00 & 41,82 & 47,27 & 10,91 & 100,00 \\ 6PAT013 & 0,00 & 5,56 & 29,63 & 44,44 & 20,37 & 100,00 \\ 6SOC016 & 3,6

[R] Randomization

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

chf0.14800 0.158 vf 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] means in tables

(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 Un

[R] SNPRelate problems

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 # Arq

[R] Kinship2 and GenABEL

p(pedig) > (mod1 = 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, --------

[R] Extracting columns

who need to apply the commands above 22 times. 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] Interaction

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

Re: [R] factor or character

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,

Re: [R] factor or character

ou are trying to find > colnames in a variable coded as integers (a factor, 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('rs

[R] factor or character

vels "rs10058955_A",..: 89 59 88 73 40 35 97 55 87 204 ... 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,

Re: [R] Reading multiple files

t, then I think 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 erro

Re: [R] Reading multiple files

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" wrote on 09/26/2012 09:11:33 AM: >> >> Hi, &

[R] Reading multiple files

="~/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("charact

Re: [R] Common elements in columns

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 Cesar da Costa > wrote: >> Hi Arun, >> >> it's exactly what I wanted. >>

Re: [R] Common elements in columns

E H J > A.K. > > > > - Original Message - > From: Silvano Cesar da Costa > To: r-help@r-project.org > Cc: > Sent: Sunday, September 2, 2012 7:05 PM > Subject: [R] Common elements in columns > > Hi, > > I have 4 files with 1 individuals in each

[R] Common elements in columns

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

] pdf(paste("plot", 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

Re: [R] NULL column

2.4252 1.837208 5.509042e-061 > > Bill Dunlap > Spotfire, TIBCO Software > wdunlap tibco.com > > >> -Original Message- >> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] >> On Behalf >> Of arun >> Sent: Wednesday,

[R] NULL column

1351 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

[R] regression analysis

= 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

valor 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-projec

[R] Challenge

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 the

Re: [R] Longitudinal data

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" To: "Silv

[R] Longitudinal data

use R. 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/ma

[R] Notation

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] Problema with Excel files

hat's happening. Can anyone help me? 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://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] Fit Gompertz' curve'

(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=&

[R] Using Function

CTOTAL), desc(CHDL), desc(CLDL), 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 Estadu

Re: [R] Tables and merge

- Original Message - From: "Silvano" To: 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 tables for each question b

[R] Pie chart

htYellow", "lightgreen", 'lightblue', 'white'), clockwise=F, radius=.7) Thanks, -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346

[R] Boxplot in order

, 11), Resp = valor ) boxplot(Resp~Groups, order=T) doesn't work. How do this? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list https

[R] Hartley's table

28 7 12.6 13.535307 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 Co

Re: [R] HLM Model

ade+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 Fo

[R] subset factor?

umeric? 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 the posting guide http://www.R-proj

Re: [R] subset

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&

[R] subset

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

[R] Tukey's table

ot;) tukval require(xtable) xtable(tukval) Some suggest? -- 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/lis

[R] step

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 De

Re: [R] Games

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

Re: [R] Games

sample(55,5)) (round11 = sample(55,5)) (round = rbind( round1, round2, round3, round4, round5, round6, round7, round8, round9, round10, round11 )) -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -

[R] Games

ample(11,2) sample(11,2) sample(11,2) sample(11,2) sample(11,2) ) but, isn't work very well. Some suggestion? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R

[R] xtable

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

[R] xtable with Sweave

quot;top", table.placement='H') @ 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

Re: [R] Table and Sweave

. -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: "Ista Zahn" To: Cc: "Silvano" Sent: Tuesday, May 11, 2010 5:48 PM Subject: Re: [R] Table a

[R] Table and Sweave

) medias.idades = tapply(IDADE, SEXO, mean) @ <>= 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 _

Re: [R] Extracting a subset

52.0 6 6 2 53.0 6 6 2 96.0 -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 -- - Original Message - From: "Ivan Calandra" To: Sent: Monday, May 10,

[R] Extracting a subset

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] Latex and Stangle()

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()

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 Cesa

[R] About logistic regression

0.6160.538 grau 0.562717 0.433615 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 Universida

[R] Mosaic

d like to make a mosaic to represent the numeric variable depending on 3 variables. Does anyone know how to do? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.org mailing list

[R] Hartley's table

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

[R] Normality in split-plot design

? 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 the

[R] Split Plot and Tukey

/Trat)) 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

[R] About scales in graphics

#x27;Meses', ylab='Número de atendimentos', col='LightYellow',ylim=c(400,450), yaxt='n', space=0.6) axis(2, at=seq(400,450, by=10), las=1) abline(h=400, col='black', lwd=1) Thanks, -- Silvano Cesar da Costa Departamento de

Re: [R] About graphics

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" To: &qu

[R] About graphics

quot;, bty='l', col="LightYellow", horiz=T, xlim=c(0,120)) abline(v=0) Suggestions? -- Silvano Cesar da Costa Departamento de Estatística Universidade Estadual de Londrina Fone: 3371-4346 __ R-help@r-project.or

[R] Bartlett' test

hanks, -- 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

Re: [R] Tests in Anova

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

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

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" To: Cc: "

[R] boxplot, data frame and html

iteImage('caixa.jpg', p, br=TRUE) hwrite('',p, 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 ___