提供两种思路：
思路1：排除非0元素，然后计算平均值：

方差类似。
思路2：将“0“替换成NA，然后就可以用mean 或者var 中的na.rm=T来计算。

万分感谢~~~{:3_46:}

大神，麻烦再帮我看看我这个代码问题在哪里？运行不出来~~~
read.table("D:/Program Files/R/R/R-2.15.2/bin/i386/test.csv",header=T,sep=",")->mydata
mydata->A          #A是一个矩阵
Ex=function(x){
mean(x)
}                          #Ex这里是一个函数
En0=function(x){
x1=x-Ex                 #所以x-Ex是不是就不能直接计算
(1/sqrt(2*pi))*(1/n)*sum(abs(x1))
}
He0=function(x){
var(x)
sqrt(abs(var-En0^2))
}
Ex=function(x){
mean(x[x!=0],na.rm=T)
}
En=function(x){
En0(x[x!=0],na.rm=T)
}
He=function(x){
He0(x[x!=0],na.rm=T)
}
apply(A,2,Ex)
apply(A,2,En)
apply(A,2,He)

你可以把你的test.csv文件用xls传上来，你没必要拿一串代码来测试我对代码的熟悉程序，别人的时间也很宝贵，只有运行才能比较快地发现你的代码问题之所在。
还有，apply(A,2,Ex)这个可以运行吧？
你要问问题，要有数据，有代码，还要有报错。
只有这样，别人才能尽快地给你解答，而不是一个问题，需要来回问3、4次。

好吧~懂了  刚来论坛，不懂规矩~受教了