similar to: hist(1) results in a segmentation fault (at least on my machine) (PR#86)

I use knncat to make a predictive model and get misclass rate > knncat.m<-knncat(training.new,k=c(10,20),classcol=5) > knncat.m Training set misclass rate: 36.88% then I try to calculate prediction accuracy by the following: > pr.knncat.train <- predict (knncat.m,training.new,training.new,train.classcol=5,newdata.classcol=5) > tb.knncat.train <-table (pr.knncat.train,

Dear all, I have done further modifications on the histogram functions that I reported earlier this year, and I hope this can be of use and perhaps included in the distribution. I have been using this stuff a couple of months myself, and while it is nothing sophisticated, it has it's applications. :-) I did a few small modifications today to make it a bit more compact. I have modified the

/usr/src/sys/modules/hifn/Makefile appears to need an update to make it same as HEAD. Same with ..modules/ubsec/Makefile. Patches: Index: Makefile =================================================================== RCS file: /home/ncvs/src/sys/modules/ubsec/Makefile,v retrieving revision 1.2.2.1 diff -u -r1.2.2.1 Makefile --- Makefile 21 Nov 2002 23:38:47 -0000 1.2.2.1 +++ Makefile 5 Jun 2003

--- filters/delay/nbdkit-delay-filter.pod | 8 ++++++++ filters/delay/delay.c | 26 ++++++++++++++++++++++++++ 2 files changed, 34 insertions(+) diff --git a/filters/delay/nbdkit-delay-filter.pod b/filters/delay/nbdkit-delay-filter.pod index c2eb172..2e2ac74 100644 --- a/filters/delay/nbdkit-delay-filter.pod +++ b/filters/delay/nbdkit-delay-filter.pod @@ -11,6 +11,7 @@

Hi Folks, I'd like to know how hist() decides how many cells to use when it ignores my "suggestion" to use say 'hist(...,breaks=50)'. More specifically, I have the results of 10000 simulations, each returning an 8-vector, therefore 8 variables each with 10000 values. Some of these 8 have somewhat skew distributions. Say one of these 8 variables is X. I ask for H <-

I discovered a bug in hist(). Try the following: x<-c(-5,-4,-4,-4,-3,-3,-3,-3,-2,-2,0,0,0,0,1,1,1,3,3,5,6) # Note that sum(x)<0: sum(x) # [1] -13 hist(x) # looks ok hist(x,freq=F) # negative bars !! # and finally this gives not 1: sum(hist(x,plot=F)$rel.freqs) # [1] -0.8076923 The reason is, that "sum(x)" is used instead of "length(x)" in the following line near

Hello: I am trying to work with a couple of microarray data sets, using platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.1 year 2003 month 11 day 21 language R In the shortcut for invoking R I have set

I have question about log2 transformation and performing mean on log2 data. I am doing analysis for ELISA data. the OD values and the concentration values for the standards were log2 transformed before performing the lm. the OD values for samples were log2 transformed and coefficients of lm were applied to get the log2 concentration values. I then backtransformed these log2 concentrations and the

At the risk of being beaten about the face and body, can somebody explain why the middle example: log2(2^3); floor(log2(2^3)) is different than examples 1 and 3? > log2(2^2); floor(log2(2^2)) [1] 2 [1] 2 > log2(2^3); floor(log2(2^3)) [1] 3 [1] 2 > log2(2^4); floor(log2(2^4)) [1] 4 [1] 4 > DrC [[alternative HTML version deleted]]

-------- --Multipart_Thu_Apr__8_13:33:51_1999-1 Content-Type: text/plain; charset=US-ASCII Hello guys, I always wondered why smbclient keeps timestamp in `put'ting a file but does not keep it in `get'ting a file. Alternative to keep the time stamp in `get'ting file is to use -Tc option or smbtar script, but they are a bit hairy for interactive use, and it can only dump files with an

Hi list, Can somebody explain why there are these warning messages?? I just don't get it.? I'm using R 2.15.1 on WinXP. Thanks! Tao > x [1] -2.143510 -1.157450 -1.315581? 1.033562 -1.225440 -1.179909 >? ifelse(x>0, log2(x), -log2(-x)) [1] -1.099975 -0.210950 -0.395700? 0.047625 -0.293300 -0.238675 Warning messages: 1: In ifelse(x > 0, log2(x), -log2(-x)) : NaNs produced

I figured out something new that I would like to see if I can do this more easy with R then Excel. I have these huge files with data. For example: DataFile.csv ID Name log2 1 Fantasy 5.651 2 New 7.60518 3 Finding 8.9532 4 Looeka -0.248652 5 Vani 0.3548 With like header1: ID, header 2: Name, header 3: log2 Now I need to get the $ID out who have a &log2 value higher then 7. I

Hello. My name is Akashah. i work at metabolic laboratory. From my study, i found that volcano plot can help a lot in my section.  i already studied about the volcano plot and get the coding to run in R software, unfortunately, there is may be something wrong with the coding. This is because  no graph appear, but no error (blue color text) was shown on the R console. Below is the coding for

Just noticed nobody has commented on this, so I took a look, and I can't see the problem (2.1.0 looks the same as 2.1.1 to me). What is supposed to be wrong? Duncan Murdoch On 5/23/2005 2:35 PM, pelzc at ohsu.edu wrote: > Full_Name: Carl Pelz > Version: 2.1.0 & 2.0.1 > OS: Windows 2000 > Submission from: (NULL) (137.53.64.158) > > > # Paste this all into an R

Hello everybody! I have 2 >issues< concerning methods applied to missing data. I think they're bugs, but who knows. 1. var(NA) returns Error in var(NA) : missing observations in cov/cor instead of NA. I expanded the summary-function to my.summary including SDev, in order to use it with tapply, which crashes in case of groups with no valid data. 2. is a similar problem. I use

HELLO, anybody... could you help me to check the below coding for volcano. what is the mistake? what the plot could not display? # volcano_plot.r # # Author: Amsha Nahid, Jairus Bowne, Gerard Murray # Purpose: Produces a volcano plot # # Input: Data matrix as specified in Data-matrix-format.pdf # Output: Plots log2(fold change) vs log10(t-test P-value) # #

I submitted a series of patches to convert all uses of log2 values to non-log2 values (That was harder than I thought because the types of the two are the same. They only different in meaning. So it was not easy to distinguish them. I split up patches so that anyone can verify correctness of the conversion that I've done.) >From now on, please always use non-log2 alignment values

Dear list, I have a data frame like: > log2.ratios[1:3,1:4] Clone a1 a2 a3 1 GS1-232B23 -0.0207500 0.17553833 0.21939333 2 RP11-82D16 -0.1896667 0.02645167 -0.03112333 3 RP11-62M23 -0.1761700 0.08214500 -0.04877000 how to make it to look like: > log2.ratios[1:3,1:4] a1 a2

This is very minor, but shouldn't log2(quote(1:10)) throw an error,the same as log() and other math functions do? It looks like log2 and log10 evaluate a call object instead of throwing a non-numeric-argument error. They do object to non-call language objects, like expressions. > log2(quote(1:10)) [1] 0.000000 1.000000 1.584963 2.000000 2.321928 2.584963 [7] 2.807355 3.000000

Hello, I'd like to fit a user defined function to a data set, but I have problems to find my problem. The user defined function is a combination of two rectangular functions, and the listing below gives an example for what I want to do. The problem is, that I get the error message for fit1 and fit2 "Error in nlsModel(formula, mf, start, wts) : singular gradient matrix at initial