Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error: $ operator is invalid for atomic vectors #12

Closed
jhurst5 opened this issue Nov 27, 2017 · 5 comments
Closed

Error: $ operator is invalid for atomic vectors #12

jhurst5 opened this issue Nov 27, 2017 · 5 comments

Comments

@jhurst5
Copy link

jhurst5 commented Nov 27, 2017

Hi Xiaolei,

I am getting the following error

[1] "Current loop: 2 out of maximum of 10"
[1] "Optimizing Pseudo QTNs..."
Error: $ operator is invalid for atomic vectors
In addition: Warning messages:
1: In ID.GP/bin.size :
 longer object length is not a multiple of shorter object length
2: In if (!is.na(inclosure.size)) { :
  the condition has length > 1 and only the first element will be used
3: In parallel::mclapply(1:m, seqQTN.optimize.parallel, mc.cores = ncpus) :
  all scheduled cores encountered errors in user code
Execution halted

I have tried changing minor things, but cannot figure out what the issue is. Any help would be appreciated.

@jhurst5
Copy link
Author

jhurst5 commented Nov 27, 2017

Here is the input:

genotype = attach.big.matrix("/path/mvp4.hmp.geno.desc")

phenotype = read.table("/path/mvp4.hmp.phe", head=TRUE)

map = read.table("/path/mvp4.hmp.map", head=TRUE)

myMVP = MVP(
phe=phenotype,
geno=genotype,
map=map,
nPC.FarmCPU=3,
perc = 1,
priority="speed",
ncpus=5,
vc.method="EMMA",
maxLoop=10,
method.bin="EMMA",
threshold=0.05,
method="FarmCPU"
)

@jhurst5
Copy link
Author

jhurst5 commented Nov 27, 2017

In addition, when I set Max Loop to 1, it runs to completion

@XiaoleiLiuBio
Copy link
Collaborator

Could you please check if there are non numeric records in second and third column in map file?

@jhurst5
Copy link
Author

jhurst5 commented Dec 4, 2017

No, there are no non-numeric records.

@XiaoleiLiuBio
Copy link
Collaborator

Please send me (xiaoleiliu@mail.hzau.edu.cn) a demo data that can repeat the error, I will test and debug it. Thank you.

@hyacz hyacz closed this as completed Sep 12, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants