Thank you for responding floramaria.

Yes, I am not totally sure how we will be able to cope with a whole bunch of stark raving normal types.

One of the aspects of the aducan trials that perhaps has not been emphasized enough on this thread is that 0.372 SUVR is the end of the road. 0.372 SUVR = 0! This may or may not be actually true, though as of now that is what is on the table. We can complain all we want, but there it is. We can't ex...

NF52, do you know of any online AD mouse databases? It would be best to not reinvent the wheel if this can be avoided. It was only recently that I heard of research into what might be our specific family AD gene. I suppose that the lab then starts by designing some CRISPRed mice. I actually was quit...

Thank you for replying SusanJ. Yes, I think that we now need to be aware that genetic uplift is starting. Life will never be the same. Natural reproduction is no longer plausible. Might not be that long before it is illegal. The idea of simply spinning the genetic roulette wheel is already absurd. H...

Code Chunk rough # Another stylized chunk that might help out. # Idea is to take in the rsnums from the polygenic score and find there genotypes # in the vcf file. Search through the vcf file and find genotypes. # Put all of this into a new array poly. # Poly would have everything: It would have the...

Code Chunk rough
# All that is needed to calculate the polygenic score is a scalar multiplication.
# Need to write a function that takes a set of rsnumbers from the polygenic set and then finds them in the vcf file



```{r}
beta<- c(1,2,1,4,5)
intgt[1:5]*beta
```

Owwh, I was worried that I was going to have to bluff it and fill up a page of code simply because I didn't know what I was doing. It looks like I can simply avoid all that mess with a few lines of code. I added them above. This is the current output (below). As we see when we add across rows with t...

# Code Chunk 2 # All this code does is extract the genotypes for each location as separate character strings # I am not entirely sure whether this is the best approach for this; I thought why not # put this into an array? We'll stick with this version for now; if trouble happens plan b. # First pair...

# Code Chunk 1 # Documented version of the code. # This just sets things up. # Install the vcfR package only once; use the library call every time you use the program # Sometime I will provide everyone a J1.vcf file to follow along # Will need to download R and Rstudio for the code to work. # Then o...

Thank you flora. The code is not that impressive, though I wanted to lift it off my hard drive before it became unfindable and completely incomprehensible. The goal here is to have a nice and simple program that would take a list of rsnumbers with their associated betas and output a polygenic score....