similar to: Interaction

Hi, I have 35 data files for reading. I would like get a program for performing reading of 35 files at once. All are of the type: Dados1.raw, Dados2.raw and so on. If the files have the same number of columns, I can read with the following commands: rm(list=ls()) filenames = list.files(path="~/Silvano/Arq", pattern="Dados+.*raw") names = substr(filenames, 1, 7) for(i in

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) But, I have a data set with 1.000 columns and I need extract 70 to use (like snps in command above). This 70 snps are in a file. So I create a file to extract them with

Hi, I have 22 files (A1, A2, ..., A22) with different number of columns, totaling 10,000 columns: c1, c2, c3, ..., c10000 I have another file with a list of 100 columns that I need to extract. These 100 columns are distributed in 22 files. How to extract the 100 columns of the 22 files? I have done it "manually" with the following commands, for example: cromo1 = read.table ("~

Hi. I have 2 tables, with same dimensions (8000 x 5). Something like: tab1: V1 V2 V3 V4 V5 14.23 1.71 2.43 15.6 127 13.20 1.78 2.14 11.2 100 13.16 2.36 2.67 18.6 101 14.37 1.95 2.50 16.8 113 13.24 2.59 2.87 21.0 118 tab2: V1 V2 V3 V4 V5 1.23 1.1 2.3 1.6 17 1.20 1.8 2.4 1.2 10 1.16 2.6 2.7 1.6 11 1.37 1.5 2.0 1.8 13 1.24 2.9 2.7 2.0 18 I need generate a table of averages, the

Hi, I have to do 10,000 linear regression analysis, and the response variable (RESP) is the same for all independent variables (10,000). y ~ x[i] i = 1, ..., 10000 For each analysis must extract the p-value and put them in an orderly increasing. I thought an analysis of the type: ana = numeric(10000) for(i in 1:10000){ mod = lm(RESP~x[i] p-value[i] = summary(mod)$coe[2,4] } Could

Hi, I want to build the table of a football league with 11 teams. All play together. So will 55 games. Since there are an odd number of teams in each round a team will not play. The games will be: games = urnsamples(1:11, x = c('A','B','C','D','E','F','G','H','I','J','K'), size=2, replace=F,

Hi, I have this data set: obitoss = c( 5.8,17.4,5.9,17.6,5.8,17.5,4.7,15.8, 3.8,13.4,3.8,13.5,3.7,13.4,3.4,13.6, 4.4,17.3,4.3,17.4,4.2,17.5,4.3,17.0, 4.4,13.6,5.1,14.6,5.7,13.5,3.6,13.3, 6.5,19.6,6.4,19.4,6.3,19.5,6.0,19.7) (dados = data.frame( regiao = factor(rep(c('Norte', 'Nordeste', 'Sudeste', 'Sul', 'Centro-Oeste'), each=8)), ano =

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 Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346

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 Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346

Hi, I need construct box plot graph, but I want keep Groups order karla = data.frame( Groups = factor(rep(c('CPre','SPre','C7','S7','C14','S14','C21','S21'), 11)), Time = rep(c(0,7,14,21), 11), Resp = valor ) boxplot(Resp~Groups, order=T) doesn't work. How do this? -------------------------------------- Silvano Cesar da

Hi, I build a graph taking into account the times: 1, 2,4,6,8,10,12,15,18,21,24,28,32 and 48. Be that the scale of the X axis does not look right. It seems equidistant. (graph attached) What changes have I to do in the following commands so that the scale be correct? interaction.plot(Tempo, Trat, Valor, ylim=c(0, 2), las=1, lty=c(1,2,3,4), lwd=3, bty='l',

----- Original Message ----- From: "Silvano" <silvano at uel.br> To: <r-help at 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 > tables for each question by CODE. >

Hi, I'm building a graph (barplot) in which the X axis label disappears. I tried to use the option mgp of par() and I could not get the desired result. Note that want the axis labels horizontally. caes = c(37,20,19,16,75,103) names(caes) = c("Pinscher", "Pastor \n Alem?o", "Poodle", "Rottweiller", "SRD", "Outros") caess =

Hi, I'm building Tukey's table using qtukey function. It happens that I can't get the values of Tukey's one degree of freedom and also wanted to eliminate the first column. The program is: Trat <- c(1:30) # number of treatments gl <- c(1:30, 40, 60, 120) # degree freedom tukval <- matrix(0, nr=length(gl), nc=length(Trat)) for(i in

Hi, I have a dataset with many variables and observations. The variable Group has two levels: C and P, the Month variable has four levels: 0, 1, 2 and 3. I want to extract a subset of the variable Weight, considering only 1 and 3 levels for Months of the Group variable. I tried the command subset but it did not work. Any suggestions? Thanks, -------------------------------------- Silvano

Hi, I did the analysis of variance of a split-plot and the effect of treatment was significant. I would like compare treatment means using Tukey. I can't extract the mean square to apply HSD.test to use in agricolae package. anava = aov(ganhos ~ Blocos + Trat*Supl + Error(Blocos/Trat)) names(anava) summary(anava) require(agricolae) HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05)

Hi, I have a dichotomous variable (Q1) whose answers are Yes or No. Also I have 2 categorical explanatory variables (V1 and V2) with two levels each. I used logistic regression to determine whether there is an effect of V1, V2 or an interaction between them. I used the R and SAS, just for the conference. It happens that there is disagreement about the effect of the explanatory variables

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

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

Hi, I'm using Sweave to prepare a descriptive report. Are at least 20 tables built with xtable command of kind: <<echo=F, results=hide>>= q5 = factor(Q5, label=c("N?o", "Sim")) (q5.tab = cbind(table(q5))) @ <<echo=F, results=tex>>= xtable(q5.tab, align="l|c", caption.placement = "top", table.placement='H') @ I'm