similar to: computing and showing mean

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 subset a matrix with 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") Z.matrix<-as.matrix(dat) Y<-array(dat,dim(,[46:63]) but I keep getting the error message: error in dim

hi all i have a function that uses two inputs, say xdata and ydata. An example is the following, simple1<-function(xdata,ydata) { ofit<-lm(ydata~xdata) list(ofit) } say i use arbitray number for xdata and ydata such that D = x1 x2 y 1 1 10 2 6 6 3 10 7 x<-D[,1:2] and y<-D[,3] if one uses these inputs and rund the program we get the following: >simple(xdata=x,ydata=y)

Hello all, I have been using R for about 3 weeks and I am frustrated by a problem. I have read R in a nutshell, scoured the internet for help but I either am not understanding examples or am missing something completely basic. Here is the problem: I want to plot data that contains dates on the x axis. Then I want to fit a line to the data. I have been unable to do it. This is an example of

I have a set of data in the form x1, y1, z1 x1, y2, z2 ... x1, yN, zN x2, y1, z(N+1) x2, y2, z(N+2) ... x2, yN, z(2N) ...and so on... xM, yN, val(M*N) I have been trying to figure out how to get R to use this data in a persp plot. So far the only thing that I can figure out to do is to break the data file into three different files. The first file contains the x-coordinate data: x1 x2

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

Hi! When running the following code I get the message: Error in eval(expr, envir, enclos) : object "A" not found fm <- function(p,ydata, env = parent.frame()) { #fit model y ~ (b0+b1x1+...+bpxp)*exp(g2plus*z2plus) where bi and g2plus are parameters #create design matrix for linear part fo00 <-

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

Airon, I used R2.4.0 on a Windows XP (SP2) (not Chinese :-)) and it still works: > data = read.csv("u:/airon.csv") > xdata = data$Adj..Close > modwt.la8 = modwt(xdata, "la8", n.level=6) > summary(modwt.la8) Length Class Mode d1 1467 -none- numeric d2 1467 -none- numeric d3 1467 -none- numeric d4 1467 -none- numeric d5 1467 -none- numeric d6 1467

Hello I am sorry fuch such a stupid question. Suppose I have a table of data having a lot of NAs and I want to replace those NAs by the mean of the column before NA replacement. How is it possible to do that efficiently ? Thanks in advance, Julie -- Julie Bernauer Yeast Structural Genomics http://www.genomics.eu.org

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)

Dear all: I have the following codes: Xdata<-c(2,3,8,9,10) Ydata<-1:5 gap.plot(Xdata, Ydata,gap=c(5,6),gap.axis="x",type="o") However, the type='o' seems only work on the first part of gap plot, the second half of the plot always just points, you can not add lines on that part, any help will be highly appreciated. I would like to have these two parts of

Hi: I have encountered problems with imodwt and universal.thresh.modwt and cannot find any reference in R Search. Hope someone can give me some ideas: Starting with modwt.la8 <- modwt(xdata, "la8", n.level=6) <-- this seems to work fine (1) ydata <- imodwt(modwt.la8) will always give ydata as numeric(0) (no values) instead of being a time series data with

Hi, all I am trying to make a plot with a axis break and I want the whole plot to be line, not points. However, when I execute the following command half of the graph is points and the other lines. gap.plot(Xdata, Ydata,gap=c(5,6),gap.axis="x",type="l") I think it might be a bug in plotrix. I would greatly appreciate your input. If there is another way to do it, I

Dear list members. I have problems with the usage of cv.glm from the boot package. Here are some parts of the script I wanted to use: data <- read.table("selected_2D.csv", header=TRUE, sep=",") ? glm.fitted <- glm("ydata$ y ~ 1 + density + vsurf_ID6 + vsurf_S ", data=data) error <- cv.glm(data=data, glm.fitted, K=6) ydata$y is a separate data set, where

Greg, Two thoughts: 1. It might be possible that 'vars' is a reserved word of sorts and if you change the name of your vector RF might be happier 2. A way that works for me is to call importance as follows: sel.imp <- importance(sel.rf, class=NULL, scale=TRUE, type=NULL) and then use the 'names' of the imp data frame to be absolutely clear to RF you are talking about the

Dear R-enthusiasts, What a wonderful package R is, many thanks to all you have contributed (I've just started using it)...I have a few questions that I am having some trouble finding the answer to: I am using RPgSQL to grab some data from a postgresql database (about 179 columns and a few thousand rows), and I want to do a lot of pairwise comparison of the rows...(this is yeast genomic

Airon, I don't think you have to find an English computer 'cause the following must work in your Chinese one :-) Let me explain. First of all, change your lines to xdata <- ckhdat$Adj..Close[1:1447] #names(ckhdwt.la8) <- c("w1", "w2", "w3", "w4", "w5","w6", "v6") note the # sign, i.e., DO NOT change the names

Hello, I have a y-axis label that reads: "S. schenckii yeast cells". The part that reads "S. schenckii" needs to be in italic style, the rest of the text is normal style. How can I specify the different font styles for each part of the y-axis label? Thank you, Judith

Hello there, I was wondering if anyone might be able to help me as I'm pretty new to R. I'm trying to create a boxplot from a data table returned from an sproc. I have the following code, which generates the plot as I'd like it: library("RODBC"); conn <- odbcConnect("datawarehouse"); # connect to datawarehouse results <- sqlQuery(conn, "call