similar to: how to add these "axis" label？

Dear all, I''m trying to write a function, that will take as an argument, some aligned genome sequences, and using a sliding window, do pairwise comparisons of sequence similarity. Coding the sliding window I think I can manage but what I''m trying to get to grips with is getting it so as every pairwise comparison is made, no matter how many genomes are added, from 3 to N. So if

Hi I am using the following script to average a set of data 0f 62 columns into 31 colums. The data consists of values of ln(0.01) or -4.60517 instead of NA's. These need to be averaged for each row (i.e 2 values being averaged). What I would I need to change for me to meet the conditions: 1. If each run of the sample has a value, the average is given 2. If only one run of the sample has a

Dear R users, I have a data generated as the following, dat <- data.frame(matrix(sample(24), nrow=4)) dimnames(dat) <-list(rownames=c('g1','g2','g3','g4'), colnames=c("A_CH1","A_CH2","B_CH1","B_CH2","C_CH3","C_CH3")) » dat A_CH1 A_CH2 B_CH1 B_CH2 C_CH3 C_CH3 g1 16 24 7 9 14

On how to even search for this "feature" as I have no idea on what it can be. I've got a meridian linked to * (by EuroISDN) which is linked to a ISDN30. I can make calls from the meridian, and receive calls into the meridian. Great stuff. However, if someone dials an invalid number, then instead of hearing a "three tone", the line just drops and goes dead. The console

Hola Javier, Intenta los siguientes cambios g4 <- function(codigo_llega, n_caracteres){ codigo_llega <- as.character(codigo_llega) if(n_caracteres == 6) res <- substr(codigo_llega, start=0, stop=4) else res <- codigo_llega res } g4 <- Vectorize(g4) x <- c('Jorge Velez','Javier Marcuzzi','Daniel Merino','Susana deus Alvarez', 'Carlos

There was a BioConductor thread today where the poster wanted to find pairwise difference between columns of a matrix. I suggested the slow solution below, hoping that someone might suggest a faster and/or more elegant solution, but no other response. I tried unsuccessfully with the apply() family. Searching the mailing list was not very fruitful either. The closest I got to was a cryptic chunk

Carlos: Te agradezco mucho tu rápida respuesta y mucho me apena haber planteado tan mal el problema. Porque la matriz en realidad es: g1 g2 g3 g4 g5 g6 g7 g1 0 18 13 16 11 12 15 g2 18 0 25 13 22 16 10 g3 13 25 0 28 23 13 25 g4 16 13 28 0 6 7 3 g5 11 22 23 6 0 18 17 g6 12 16 13 7 18 0 8 g7 15 10 25 3 17 8 0 Entonces cada cantidad debe conservar la pertenencia al grupo

Hello, Were having issues with macs writing to our gluster system. Gluster vol info at end. On a mac, if I make a file in the shell I get the following message: smoke:hunter david$ echo hello > test -bash: test: Operation not permitted And the file is made but is zero size. smoke:hunter david$ ls -l test -rw-r--r-- 1 david realise 0 Mar 3 08:44 test glusterfs/nfslog logs thus:

Hi all, I am having a trouble with this function I wrote ################################################### p26=function(x,alpha){ # dummy variable j=1 ci=matrix(ncol=2,nrow=3) while (j<4){ if (j==2) {x=x+c(-1,1)*0.5} ci[j,]= x+qnorm(1-alpha/2)^2/2+ c(-1,1)*qnorm(1-alpha/2)* sqrt(x+qnorm(1-alpha/2)^2/4) j=j+1 if (j==3) { # exact x=x-c(-1,1)*0.5

Hi All, I’m trying to evaluate the frequency of different strings in each row of a data.frame : INPUT: ID G1 G2 G3 G4 … GN 1 AA BB AB AB … 2 BB AB AB AA … 3 AC CC AC AA … 4 BB BB BB BB… The number of different strings can vary in each row. My solution has been: for (i in 1:length(INPUT[,1])){ b=as.data.frame(table(t((INPUT[i,2:5])))) <some operations using the string values and

Hi, I'm new to R. I was wondering if anyone has an OS X native installation (not Carbon) of R on a G3 Macintosh that is able to run linear regression? I installed Jan de Leeuw's R binary for OS 10.2 Jaguar on my laptop, which is still a G3 model. When I try to run a simple regression, using lm, I get an "Illegal Instruction" message and R quits. Jan thinks it may be the case

Chris Mauritz wrote: > > The 32bit PPC release is for the older G3/G4 macs, no? I've got a > pile of unused G4 400-800mhz powermacs just sitting around The ppc32 distro will work on any NewWorld G3, G4 based machine ( and 32bit ibm power, like the older rs6k's ) - it _also_ works fine on the MacMini and G4 Xserve's. With some common sense tweaks, a G4 450Mhz is a very

I'm trying to sort a DATAFRAME by a column "ID" that contains alphanumeric data. Specifically,"ID" contains integers all preceeded by the character "g" as in: g1, g6, g3, g19, g100, g2, g39 I am using the following code: DATAFRAME=DATAFRAME[order(DATAFRAME1$ID),] and was hoping it would sort the dataframe by ID in the following manner g1, g2, g3, g6, g19,

Dear R users, I have noted a difference in the merge distances given by hclust using centroid method. For the following data: x<-c(1009.9,1012.5,1011.1,1011.8,1009.3,1010.6) and using Euclidean distance, hclust using centroid method gives the following results: > x.dist<-dist(x) > x.aah<-hclust(x.dist,method="centroid") > x.aah$merge [,1] [,2] [1,] -3 -6

Hi We deployed centos 4.8 on HP DL380 G2, G3 and G4 servers (same number of cpus). The same software is running on all servers. Only on G2 servers we see a low cpu utilization combined with high load-average. For example, during the same time of day and same workload, on a G2 I see a l-avg of 11.31 and on a G4 I see 2.68. The server response time seems unchanged and we don't have any problem

I have a really simple question I have a data frame of 8 variables (the first column is the subjects' id): SubID G1 G2 G3 G4 W1 W2 W3 W4 1 6 5 6 2 6 2 2 4 2 6 4 7 2 6 6 2 3 3 5 5 5 5 5 5 4 5

Dear R helpers,   Following is a R script I am using to run the Fast Fourier Transform. The csv files has 10 columns with titles m1, m2, m3 .....m10.     When I use the following commands, I am getting the required results. The probelm is if there are 100 columns, it is not wise to define 100 commands as fk <- ONS$mk and so on. Thus, I need some guidance to write the loop for the STEP A and

Hi, I have a glm gives summary as follows, Estimate Std. Error z value Pr(>|z|) (Intercept) -2.03693352 1.449574526 -1.405194 0.159963578 A 0.01093048 0.006446256 1.695633 0.089955471 N 0.41060119 0.224860819 1.826024 0.067846690 S -0.20651005 0.067698863 -3.050421 0.002285206 then I use confint(k.glm)

Dear R Gang: I have a Macintosh G4 powerbook, running OS 9.2. Last week I ran a statistical analysis program that I wrote in R for two days continuously. The program does repeated intensive computations on a large data set. The program completed its task successfully, but during the process I noted that the computer seemed to be running unusually hot. Shortly after the program stopped,

We were looking at conducting a pilot program with Dell to virtualize our lab environments. We are specifically looking at Dell/Citrix provisioning servers. The OS in this setup is streamed from the server on every boot. I was curious if the R Project had any licensing stipulations in this type of setup. This is in the beginning phases and we wanted to answer some questions about licensing.