r - How to calculate rowMeans for dataframe? -

-

i trying calculate rowmeans dataframe using these command

stamrow  <- data.frame(probeid=stam[,1], means=rowmeans(stam[,c(-1,-2,-3)])) stam$means  <- rowmeans(stam[,-(1:3)])

but getting error in both case

rowmeans(stam[, -(1:3)]) : 'x' must numeric

my dataset this: 

 stam =  probeid                     chr     position        sample1 sample2 sample3 sample4 sample5 chr10fs00300029     10     3000293             0.132    0.135   0.312   0.724   0.889    chr10fs003018825        10     3018825             0.524       0.446    0.203   -0.022  0.581

thanking in advance

it works data posted. check str(stam) see if columns numeric. 

data.frame(probeid=stam[,1], means=rowmeans(stam[,-c(1:3)]))  #         probeid  means  #1  chr10fs00300029 0.4384  #2 chr10fs003018825 0.3464

update

if numeric columns start 4 15, can convert columns factor class numeric first

 stam[,4:15] <- lapply(stam[,4:15], function(x) as.numeric(as.character(x)))

and apply rowmeans

or can specify numeric in colclasses argument in read.table

stam <- read.table('file.txt', header=true, sep='',        colclasses=c(rep('character',3), rep('numeric',12)))

data

stam <- structure(list(robeid = c("chr10fs00300029", "chr10fs003018825" ), chr = c(10l, 10l), position = c(3000293l, 3018825l), sample1 = c(0.132,  0.524), sample2 = c(0.135, 0.446), sample3 = c(0.312, 0.203),  sample4 = c(0.724, -0.022), sample5 = c(0.889, 0.581)), .names = c("robeid",  "chr", "position", "sample1", "sample2", "sample3", "sample4",  "sample5"), class = "data.frame", row.names = c(na, -2l))

Comments

Post a Comment

Popular posts from this blog

matlab - "Contour not rendered for non-finite ZData" -

-
i'm trying plot frequency characteristic equation using ezplot, matlab gives following warning, "contour not rendered non-finite zdata". have used command plot frequency equations warning , plot display empty , not change axis range well. can please help. appreciated. here's code i'm using. % transfer matrix case-i, thin rotor clear all; clc; ei = 1626; l = 0.15; m = 0.44108; = 2.178*10^-4; i_p = 2.205*10^-5; itr = 0.24; i_pr = 0.479; syms p n; f = [1 l*1i l^2/(2*ei)*1i l^3/(6*ei); 0 1 l/ei -l^2/(2*ei)*1i; 0 0 1 -l*1i; 0 0 0 l]; p = [ 1 0 0 0; 0 1 0 0; 0 -it*p^2+i_p*n*p 1 0; -m*p^2 0 0 1]; p_r = [1 0 0 0; 0 1 0 0; 0 -itr*p^2+i_pr*n*p 1 0; -m*p^2 0 0 1]; = f*p*f*p*f*p*f; b = p_r*f*p*f*p*f; r = a(1,2)/a(1,4); a12_p = 0; a22_p = a(2,2)-r*a(2,4); a32_p = a(3,2)-r*a(3,4); a42_p = a(4,2)-r*a(4,4); ap(2,2) = a22_p; ap(3,2) = a32_p; ap(4,2) = a42_p; ap(4,4) = 1; c = b*ap; m = [c(3,2) c(3,4); c(4,2) c(4,...
Read more

delphi - Indy UDP Read Contents of Adata -

-
i using indy udp server read string of data. however, not know how work adata parameter. can around converting binary using bytetobin function below , convert hex using bintohex1 . feel stupid , works slow. know how can directly results without 2 conversions? thanks!! here codes: procedure tform1.idudpserver1udpread(athread: tidudplistenerthread; const adata: tidbytes; abinding: tidsockethandle); var buf: tidbytes; buffer: string; data_received: string; begin if bytestostring(adata) <> 'hello' begin buffer := hextostring('00'); setlength(buf, length(buffer)); copytidstring(buffer, buf, 0); abinding.sendto(abinding.peerip, abinding.peerport, buf); data_received := bintohex1(bytetobin(adata[1]) + bytetobin(adata[2]) + bytetobin(adata[3]) + bytetobin(adata[4]) + bytetobin(adata[5]) + bytetobin(adata[6]) + bytetobin(adata[7]) + bytetobin(adata[8]) + bytetobin(adata[9]) + bytetobin(adata[10]) + bytetobin(adata[11...
Read more

javascript - Any ideas when Firefox is likely to implement lengthAdjust and textLength? -

-
any ideas when firefox implement lengthadjust andtextlength?..... working svg application. lengthadjust , textlength working fine in chrome not working in fire fox . <svg viewbox = "0 0 500 300" version = "1.1"> <defs> <path id="s3" d="m 10,90 q 100,15 200,70" /> </defs> <g> <text font-size = "20"> <textpath xlink:href="#s3" textlength="205"> short text</textpath> </text> <use xlink:href="#s3" fill="none" stroke="black" stroke-width="1"/> </g> </svg> this implemented in firefox text elements (not textpath or tspan) in case that's ok though can move attribute text element. <svg viewbox = "0 0 500 300" version = "1.1"> <defs> ...
Read more