生物信息学：基因芯片数据分析

按：股票市场经常这个仓，那个仓的，今天去baidu了一下，blog以志之。以下内容均来自百度知道，特此声明。

 

买入股票叫进仓，卖出全部的股票叫平仓。买入股票后把股票放在帐户内不卖出叫持仓。
买入股票后股价下跌，然后在低价再买入股票以滩低每股平均成本价叫补仓。
把帐户内所有的钱全部买入股票叫重仓，用一半左右的钱买入股票叫半仓，用一小部份钱买入股票叫轻仓。

那么什么叫做老鼠仓？

 

老鼠仓是指庄家在用公有资金在拉升股价之前，先用自己个人（机构负责人，操盘手及其亲属，关系户）的资金在低位建仓，待用公有资金拉升到高位后个人仓位率先卖出获利，当然，最后亏损的是公家资金。大家都知道，中国股市的特色就是无庄不成股，而老鼠仓就存在于这些大大小小的庄股当中。券商是庄股中的主力队员，利用自身具有融资的天然优势，从社会各方面融入大量资金坐庄拉升股票。坐庄本来是为了赚钱盈利的，但券商坐庄很少有真正赚钱的，原因就在于券商把股票拉升后，大量底部埋仓的老鼠仓蜂拥出货，券商又在高位接盘。这样的结果就是券商亏损累累，老鼠仓赚个钵满盆满。这便是当今券商被掏空的主要形式。

按：今天第一次听说了“蒙代尔一克鲁格曼不可能三角形”，觉得这个东西挺有意思。所以转一篇文章过来与浏览此blog的各位亲朋好友共同分享。文章原文来自百度知道，其中有一幅图已经失落了，各位看官请多多包涵。

 

关于汇率政策与货币政策的协调问题，在20世纪60年代初，弗莱明和蒙代尔认为：资本的自由流动、货币完全独立、汇率稳定三项目标中，一国政府最多只能同时实现两项。这一结论被后人称为“蒙代尔三角”。克鲁格曼在199B年初发表的一篇讨论亚洲金融危机的文章中，也谈到“蒙代尔三角”的问题，他称之为“The eternal triangle”(永恒的三角形)；国内学者易纲总结为“蒙代尔一克鲁格曼不可能三角形”，（如下图）并且提出X+Y+M=2定理。

http://www.forumcn.com/inews/upload/snatch/2004-05- 08/img_new.asp_id_71668 （图）

图中三角形坐标系三角形的三个角点表示三种不同的汇率制度。中国占住两条边，即可以独立制定货币政策和保持汇率稳定，但要实行外汇管制，资本不能自由流动。X+Y+M=2主要用来分析我国未来的货币政策和汇率政策的走势。按照定理，X、Y、M三者之和一定等于2。货币政策必须保持独立，因此M不能变，必须接近1或等于1；而随着对外开放不断扩大，资本自由流动不断增加，Y一定要增加；变化的唯一可能性就是X值要缩小，这意味着汇率要朝着有弹性的浮动汇率方向发展。这个模型常被引用来说明关于我国为什么应该实施浮动汇率制，或者说明实施浮动汇率制是发展的方向。

美国财长斯诺访华也传达了同样的观点。他在接受采访时说，中国要扩大对外开放，从中国自身的利益角度出发，应该实施由市场供求决定的灵活的汇率制度。政府不要干预过多，否则不可能存在比较独立的货币政策。关于央行的此次紧缩银根的行动，国内有看法认为，央行最终选择提高准备金率来控制过快的国内信贷增长，这种选择显然受到我国盯住美元的汇率制度的约束，汇率稳定优先的考虑使得央行在一定程度上丧失了货币政策的独立性。

事实上，目前我国货币政策并不具有完全的独立性，目前我国M并不等于1。货币政策通过利率影响生产和消费的渠道并不畅通，没有形成市场供求决定意义上的利率。按照“先利率，后汇率”的开放原则，利率的市场化比汇率放开更紧要。因此在我国建立一个有效的货币政策体系前，必须保持汇率稳定，以减少风险和波动。正如美国经济学家罗伯特·蒙代尔所言，中国货币政策最大的成就就是保持了外汇稳定。外汇稳定已成为中国实施货币政策的重要基础和指导方针。中国应当坚持稳定的汇率制度。如果人民币汇率出现大幅度波动，将给中国经济带来许多不稳定因素，并很可能使中国重蹈20世纪七八十年代日本经济衰退的覆辙。

 

央行前不久加息27个基点，有人预计中国可能由此开启了加息通道，年内可能不止一次继续加息。那么加息究竟能不能解决流动性过剩以及资产价格膨胀问题呢？

一般认为，加息会使企业融资成本增加，限制投资过热；另外，加息还使人们更愿意去储蓄，社会消费需求较少，从而促使物价下降；还有就是提高利率促使贷款资金真正流入到盈利高的行业，从而能优化资源配置。

 

但是中国的市场经济不够健全，企业行为受利率影响也有着独特的中国国情。郎咸平在2005年的财经节目上讲到过，国企受利率影响小，很多国企贷款之后根本没有考虑到要还的问题，因此不管利率是多少都照借不误。这些借出去的钱所投资的项目受市场影响也小，因此造成了很多重复建设，过剩投资的行为。而真正对市场敏感的民营企业，利率高了，他们就不敢借款了，好多该投资的项目就没法投资。这样提高利率反而造成了该降的不讲，不该降的反而降了。当然，即使郎监管的话是纯净的真理，也只是反映了2005年的情况，2007年是什么样还要具体分析。

 

论述未完，请在明天这个时候继续关注本blog。

     只可惜beyond的主将黄家驹于九三年就已经与世长辞，听他的歌更让人有种恍如隔世的感觉。鲜活的生命仍然跳跃在耳边，让人去爱，去恨，去宣泄，去放纵，可是斯人却已经一去不返了。“羊公碑尚在，读罢泪沾襟”！听其歌曲，想见其为人，同样的感慨怎能不油然而生？

 

 

 

 

在bioconductor的字典里，marray是microarray的简写，它是一个程序包，用来给cDNA芯片作标准化(normalization)。这里说一说marray中的一个class，叫做marrayRaw，它是开发者自己定义的一种数据类型。

 

marrayRaw，可能是MicroARRAY Raw data的简写，也就是芯片原始数据的意思。与marrayRaw相对应，marray里还有一个class叫做marrayNorm，这应该是MicroARRAY NORMalized data的简写，即标准化之后的芯片数据。

 

marrayRaw类型实际上是一个list,它里面包含了一些子数据元素。想看它包含了哪些数据元素，可以使用如下R命令：

slotNames("marrayRaw")

 

R这时便会给出答案如下：

[1] "maRf" "maGf" "maRb" "maGb" "maW" "maLayout" "maGnames" "maTargets" "maNotes"

 

把R给出的答案翻译成汉语就是：marrayRaw这个class里面包含有maRf（前景红光强度，也就是芯片上那些小点点的红色荧光强度），maGf（前景绿光强度），maRb（背景红光强度）,maGb（背景绿光强度）,maLayout（芯片布局）,maGnames（基因名，ID等信息）,maTargets（芯片实验设计信息）,maNotes（相当于一个附注，用途不大）。

 

芯片的原始数据往往是GenePix软件或Spot软件生成的，比如GenePix生成的数据会含有一大堆数字都标明是前（背）景红（绿）光强度,那么marrayRaw里面的maRf是哪个呢？我看了看，发现marray软件包取的是Median一项，其他各项好像都被舍弃了。

 

创建一个marrayRaw对象很简单，最常用的办法是使用marray提供的read.GenePix函数等，至于read.GenePix怎么使用，则不是这里所要讨论的内容。当然，如果不觉得麻烦，也可以通过new函数手工创建。

 

如果有了一个marrayRaw对象，假设叫做data，这个时候可以通过@操作符来获取其中的子元素，比如想把data里面的maRf的信息拿出来，可以使用如下命令：

 

Rf <- data@maRf