similar to: subarray

Dear R-expert, For histogram function, can we get the table of bin and frequency like in excel, together with the histogram? Therefore, we can check the number of data included. Thank you so much for your attention and help. ____________________________________________________________________________________ [[elided Yahoo spam]]

Hello everyone, I am trying to plot an MvA plot with the following code: dat<-read.table(file="C:\\Documents and Settings\\.....txt",header=T) file.show(file="C:\\Documents and Settings\\Ow.......txt") library(sma) data(MouseArray) dat.o<-as.list(dat) dat.o$R.2<-as.matrix(dat[,c(1:5)]) dat.o$G.2<-as.matrix(dat[,c(6:10)]) dat.o$Gb.2<-as.matrix(dat[,c(11:15)])

Hello everyone, I have a question as to what code should be used if one wanted to subset from a dataframe of 7000rows by 38 columns.For example if you wanted to generate a array of 100 by 38 would you not use the following: Z<-(dat,dim=c(100,38)) where dat is the dataframe of 7000by 38 Whenever I do this and type show(Z) I get a output of 100 by 38 they all say "integer,7000" as

Hello to everyone, I am trying to compute a mean and show it to the screen and then save to later be used as a boxplot.I have used the following code: dat<-read.table(file="C:\\Documents and Settings\\Owner\\My Documents\\Yeast\\Yeast.txt",header=T,row.names=1) file.show(file="C:\\Documents and Settings\\Owner\\My Documents\\Yeast\\Yeast.txt") x<-dat[2,23:46]

Hello everyone, I am trying to plot a boxplot but am coming up with the error : "could not find function dat" I have used the following code. dat<-read.table(file="C:\\Documents and Settings\\.......txt",header=T,row.names=1,blank.lines.skip=F,na.strings="NA") file.show(file="C:\\Documents and Settings\\......txt")

Hello everyone, I am trying to use the following code: dat<-read.table(file="C:\\Documents.....txt") class2<-read.table(file="C:\\Documents.....txt") class1<-read.table(file="C:\\Documents .....txt") dat<-as.data.frame(dat) dimnames(dat)[[1]]<-paste("g",c(1:nrow(dat)),sep="")

Hello everyone, I have two problems which I am unable to solve : 1.I am trying to add the  row labels (g1-g2000) to the very left of a data table. The data table is 2000 rows  by 62 columns.I have used the following code. read.table(file="C:\\Documents and Settings\\Owner\\My Documents\\colon cancer1.txt",header=T,row.names=1) rowname(dat) <- paste("g", c(1:nrow(dat)),

Hello everyone, I am using the following code to try to calculate the mean : dat<-read.table(file="C:\\Documents and  Settings.....txt") dat<-as.numeric(dat) x1.m<mean(dat) I am getting the following error message Error in eval.with.vis(expr,envir,enclos): (list) object cannot be coerced to type"double' I do not understand what is wrong as I thought that I have changed

Hello everyone, I have a data frame in which I am wanting to eliminate the row labels and then relabel the rows with g1-g2000.I have used the following code: dat<-read.table(file="C:\\Documents and Settings\\Owner\\My Documents\\colon cancer1.txt",header=T,row.names=1) file.show(file="C:\\Documents and Settings\\Owner\\My Documents\\colon cancer1.txt") I thought that this

useR's I have a matrix from which I want to take multiple subsets from, according to a particular scheme I will now describe. The matrix below (mat) is 5x5, and I want to take 9 subsets of it, each of dimension 3x3. The best way to explain what the result should look like is with the following: > dat <- c(3,6,1,9,12,9,2,10,6,5,3,13,1,4,8,9,4,6,10,11,2,7,3,5,10) > miss <-

Hi. I'm new to R, and I'm still learning R's system for addressing subsets of data structures. I'm particularly interested in the problem of selecting subarrays based on complex criteria involving the dimnames (as opposed to the values of the cells) of the array. Here's an example of such a problem. Suppose I have an array x of unknown dimensions (it may have been passed as

Hi, When I call the *switch* function first time, it works. but when I call it at the second time, it does nothing. The version I use is R version 2.9.0 Under development (unstable) (2009-02-21 r47969) here is the output: > organism="human" > species <- switch(organism, human <- "Hs", fly <- "Dm", mouse <-

Define an array > v<-1:256 > dim(v)<-rep(4,4) Subarrays can be obtained as follows: > v[3,2,,2] [1] 71 87 103 119 > v[3,,,2] [,1] [,2] [,3] [,4] [1,] 67 83 99 115 [2,] 71 87 103 119 [3,] 75 91 107 123 [4,] 79 95 111 127 In the general case this procedure is very tedious. Given an array A, dim(A)=(dim_1,dim_2,...,dim_d) and two vectors

Dear helpeRs- I'm using a for loop to create a series of models. I'm trying to assign a name to each model created, using the loop index. The loop gets stuck at the name of the model, giving the error "target of assignment expands to non-language object". The linear model runs without error; only the name is problematic. Here is the current loop syntax. The use of dat

Hi, Few days ago I have asked about spatial join on the minimum distance between 2 sets of points with coordinates and attributes in 2 different data frames. Simon Knapp sent code to do it when calculating distance on a sphere using lat, long coordinates and I've change his code to use Euclidian distances since my data had UTM coordinates. Typically one data frame has around 30 000 points

Hi all (really probably just Deepayan): In the plot below I want to add text on either side of each violin plot that indicates the number of observations that are either positive or negative. I'm trying to do this with ltext() and I've also monkeyed about with panel.text(). The code below is generally what I want but my calls to ltext() are wrong and I'm not sure how to fix them.

Hi, I cant understand where I am going wrong.Below is my code.I would really appreciate your help. Thanks. > genfile<-read.table("c:/tina/phd/bs871/hw/genfile.txt",skip=1) > > #read in SNP data > snp.dat <- as.matrix(genfile) > snp.name <- scan("c:/tina/phd/bs871/hw/genfile.txt",nline=1,what="character") Read 100 items

I am using the "predict" function on a support vector machine (svm) object, and I don't understand why I can't predict on a dataset with more observations than the training dataset. I think this problem is a generic "predict" problem, but I'm not sure. The original svm was fit on 50 observations.

I have a zooreg object and I want to be able to generate a value for seasons and 11-day composites paste it onto my zoo data frame, along with year, month and days. Right now I have the following to work from: eg. dat.zoo.mdy <- with(month.day.year(time(dat.zoo)), cbind(dat.zoo, year, month, day, quarter = (month - 1) %/% 3 + 1, dow = as.numeric(format(time(dat.zoo), "%w")))) For

This must be a very simple question, but I can't find any information on it elsewhere, sorry. When extracting information from a list using column headers, how do I get R to interpret something as a variable rather than a string? For example: xx$"YAL002" works, but this doesn't: gene <- "YAL002" xx$gene neither do xx$parse(gene) xx$eval(gene)