小编bde*_*vic的帖子

如何在julia中将数据从一个处理器移动到另一个处理器？

说我有一个阵列

a = [1:10]

或者其他一些数据结构.将它放在所有其他可用处理器上的正确方法是什么,以便它们可以作为相同的变量名在这些处理器上使用？

我有Rails 3

Carrierwave 0.5.4

//app/uploaders/fasta_uploader.rb

class FastaUploader < CarrierWave::Uploader::Base
  storage :file
  def store_dir
      'public/data/01_fasta'
  end
end

//移民

class AddFileUpToCvits < ActiveRecord::Migration
  def self.up
    add_column :cvits, :fasta, :string
  end

  def self.down
    remove_column :cvits, :fasta
  end
end

//app/models/cvit.rb

class Cvit < ActiveRecord::Base
    attr_accessible :fasta
    mount_uploader :fasta, FastaUploader
end

//形成

<%= form_for(@cvit, :html => {:multipart => true, :onsubmit => "return ray.ajax()" }) do |f| %>
  ...
  ...
  <%= f.file_field :fasta %><br></br>
  <div class="actions">
    <%= f.submit "Submit"%>
  </div>
<% end %>

我得到这个错误:未初始化的常量Cvit :: FastaUploader …

我喜欢的的parallel包在R和它是多么容易和直观做并行版本apply,sapply等等.

是否有类似的并行功能replicate？

我似乎无法在任何地方找到这个问题的简单明确的答案!一切似乎都已过时或不完整!

我只是希望用户能够点击链接或按钮并下载文件(即公共文件夹中的某个位置)

我试过这个:

#view
<%= link_to "Raw blast output" ,:action => :download, :file_name => "public/data/02_blastout/#{@bl_file}" %>
#controller
def download
    send_file "#{RAILS_ROOT}/#{params[:file_name]}"
end

但我得到这个错误:

No route matches {:action=>"download", :file_name=>"public/data/02_blastout/input0.fa_x_Glyma1aaunq.bl", :controller=>"cvits"}

谢谢您的帮助!!

我在R中有一个大数据矩阵.我使用包'R.matlab'将数据转换为matlab数据,如下所示:

writeMat(con="...filepath", x=data)

我没有使用matlab的经验所以请耐心等待我:

当我将数据加载到matlab时,它说我有1x1结构.

我想把它变成矩阵形式.我试过了:

data=struct2cell(x) 

但这看起来不太合适.数据是十进制数字btw.

我想基于某些行索引在一对图中着色点.这是我用于绘制1个变量与另一个变量的代码.

cases<-which(rownames(data_no_na) %in% colnames(tumor_data))
controls<-which(rownames(data_no_na) %in% colnames(control_data))

plot(y=range(pca[,1]),x=range(pca[,2]),type='n',xlab="Principle Component 2",ylab="Principle Component 1", main="Iterative Thresholding Sparse PCA")

points(y=pca[cases,1], x=pca[cases,2], col = 'red' )
points(y=pca[controls,1], x=pca[controls,2], col = 'blue' );

一个简单的对图是这样的:

pairs(pca[,1:3])

编辑:示例:

cases<-1:10
controls<-11:20

pca<-matrix(c(rnorm(3*10,0,1),rnorm(3*10,5,1)),nrow=20,ncol=3)

我有两个大文件(27k行和450k行).他们看起来像:

File1:
1 2 A 5
3 2 B 7
6 3 C 8
...

File2:
4 2 C 5
7 2 B 7
6 8 B 8
7 7 F 9
... 

我想要两个文件中的第三列都在两个文件中的行(排除了带有A和F的注释行):

OUTPUT:
3 2 B 7
6 3 C 8
4 2 C 5
7 2 B 7
6 8 B 8

什么是最好的方式？

我有以下功能:

f1<-function(x){
  iih_data<-...stuff...
  ...more stuff...

  cl <- makeCluster(mc <- getOption("cl.cores", 6))
  clusterExport(cl, c("iih_data"))
  clusterEvalQ(cl, require(lme4))

  Tstar<-parCapply(cl, ystar, function(x){
     ostar=glmer(x ~ GENO + RACE + (1|GROUP), family="binomial",data=iih_data,nAGQ=1)
     fixef(ostar)[2]/sqrt(vcov(ostar)[2,2])
  })

  stopCluster(cl)

  ...more stuff...
}

但我得到这个错误:

Error in get(name, envir = envir) : object 'iih_data' not found

我猜这与我试图在函数内运行并行应用的事实有关.你能帮帮我解决这个问题吗？谢谢

是否可以在 julia 中交错两个数组？

例如，如果a=[1:10]和b=[11:20]我希望能够返回

20-element Array{Int64,1}:
  1
  11
  2
  12
  3
  13
  4
  14
  .
  .
  .

有点类似于 ruby​​ 可以做的在 Ruby 中合并和交错两个数组

我有一个旧的S-plus脚本,我想在R中重现结果.我唯一的问题是随机种子.我知道他们使用不同的算法进行伪随机数生成.在S-plus文件中,使用以下方法设置种子:

set.seed(337)

我能找到的相关信息是

S-plus设定种子

R设定种子

从这个文档看起来,S-plus使用"Super-duper"算法进行伪随机数生成.R有这个选项,但是它表示它需要2个整数,而S-plus只需要0到1000之间的1个整数.此外,R doc说

The two seeds are the Tausworthe and congruence long integers, respectively. 
A one-to-one mapping to S's .Random.seed[1:12] is possible but we will not
publish one, not least as this generator is not exactly the same as that in
recent versions of S-PLUS.

我不太清楚这意味着什么.那么有谁知道是否有可能复制结果？

一个旧的文章将R邮件列表上试图获得S/R/SPLUS相同的结果.