On 2020-04-10 19:05 -0500, Ana Marija wrote:
> I am not sure what I am suppose to run
> from your codes. Can you just send me
> lines of codes which I should run?
> (without part where you are loading
> your data frames) (assuming my files
> are as I showed them)
Dear Ana,
try these lines:
I am not sure what I am suppose to run from your codes.
Can you just send me lines of codes which I should run? (without part
where you are loading your data frames)
(assuming my files are as I showed them)
Or the whole idea was to remove sep=" " from everywhere?
On Fri, Apr 10, 2020 at 7:01 PM R
On 2020-04-10 18:46 -0500, Ana Marija wrote:
> so if I understand correctly you would just remove sep=" " from my codes?
>
> Thank you so much for working on this.
> Is there is any chance you can change my original code (pasted bellow)
> with changes you think should work?
>
> library(SNPRelate)
so if I understand correctly you would just remove sep=" " from my codes?
Thank you so much for working on this.
Is there is any chance you can change my original code (pasted bellow)
with changes you think should work?
library(SNPRelate)
# get PLINK output
plink.genome <- read.table("plink.geno
On 2020-04-10 17:05 -0500, Ana Marija wrote:
> it didn't work unfortunately with your
> example:
>
> > plink.genome[idx,]
> character(0)
Hi! Perhaps csv formatting are better
suited for these emails ... the two
ibdlist lines I added still matches in
this example, added lookup for kinship
..
it didn't work unfortunately with your example:
> plink.genome[idx,]
character(0)
here is my whole:
plink.genome <- read.table("plink.genome", header=TRUE)
FID1 IID1 FID2 IID2 RTEZ Z0 Z1 Z2
PI_HAT PHE DST PPC RATIOIBS0IBS1IBS2 HOMHOM
HETHET
On 2020-04-10 15:38 -0500, Ana Marija wrote:
| Hi,
|
| I have this code:
Dear Ana,
none of the ID tuples in the head
outputs you provided matches, so I added
two lines in ibdlist that matches up;
perhaps if you provided more lines that
would have matched in a pastebin
somewhere ...
optio
Hi,
I have this code:
library(SNPRelate)
# get PLINK output
plink.genome <- read.table("plink.genome", header=TRUE)
> head(plink.genome)
FID1 IID1FID2 IID2 RT EZ Z0 Z1 Z2 PI_HAT PHE DST
1 fam1054 G1054 fam1054 G700 OT 0 0.0045 0.9938 0.0017 0.4986 -1 0.839150
2 fa
8 matches
Mail list logo
