Re: [R] snpStats reference allele used in genetic associations?

2014-05-27 Thread francesca casalino
requencies of the same allele found in the 1000 Genomes, I get concordance up to frequency= 0.5, then a shift in direction happens and I get discordance up to 1 for the reference frequency. Thank you very much for any suggestions you may have, Francesca 2014-05-27 2:07 GMT+01:00 David Duffy

Re: [R] snpStats reference allele used in genetic associations?

2014-05-26 Thread francesca casalino
). Thanks for your help. Francesca 2014-05-26 11:13 GMT+01:00 francesca casalino : > But then how do you know which allele is the reference and which the risk > allele (between A/T/C/G)? > > > 2014-05-26 1:41 GMT+01:00 David Duffy : > > francesca casalino asked: >>

Re: [R] snpStats reference allele used in genetic associations?

2014-05-26 Thread francesca casalino
But then how do you know which allele is the reference and which the risk allele (between A/T/C/G)? 2014-05-26 1:41 GMT+01:00 David Duffy : > francesca casalino asked: > > >> Does anyone know how to find the reference allele used for genetic >> associations ran in snpS

[R] snpStats reference allele used in genetic associations?

2014-05-23 Thread francesca casalino
Hi, Does anyone know how to find the reference allele used for genetic associations ran in snpStats? I have ran several associations using snp.rhs.tests, but I cannot tell which allele was used as the "effect allele". Is it the one coded as "Al1" in the SNP.support file? I can find the RAF (risk

[R] Unicode symbols not working with ggplot in R

2014-02-13 Thread francesca casalino
Hi, I am trying to produce a ggplot graph using specific characters in the labels, but ggplots doesn't seem to support certain symbols. For example, when I type: print("\u25E9") it shows a square which is half black, but when I try to use it in ggplot it doesn't print. I am using facet_wrap, bu

Re: [R] ggplot2 and facet_wrap help

2013-02-18 Thread francesca casalino
oes give the desired display. > > Best, > Ista > > On Mon, Feb 18, 2013 at 6:04 AM, francesca casalino > wrote: >> Dear R experts, >> >> I am trying to arrange multiple plots, creating one graph for each >> size1 factor variable in my data frame, and

[R] ggplot2 and facet_wrap help

2013-02-18 Thread francesca casalino
Dear R experts, I am trying to arrange multiple plots, creating one graph for each size1 factor variable in my data frame, and each plot has the median price on the y-axis and the size2 on the x-axis grouped by clarity: library(ggplot2) df <- data.frame(price=matrix(sample(1:1000, 100, replace =

Re: [R] Exponentiate very large numbers

2013-02-05 Thread francesca casalino
I am sorry I have confused you, the logs are all base e: ln(a) = 1347 ln(b) = 1351 And I am trying to solve this expression: exp( ln(a) ) - exp( ln(0.1) + ln(b) ) Thank you. 2013/2/4 francesca casalino : > Dear R experts, > > I have the logarithms of 2 values: > > log(a)

[R] Exponentiate very large numbers

2013-02-04 Thread francesca casalino
Dear R experts, I have the logarithms of 2 values: log(a) = 1347 log(b) = 1351 And I am trying to solve this expression: exp( ln(a) ) - exp( ln(0.1) + ln(b) ) But of course every time I try to exponentiate the log(a) or log(b) values I get Inf. Are there any tricks I can use to get a real resu

[R] Submit a new script after all parallel jobs have completed

2012-09-26 Thread francesca casalino
Dear R experts, I have an R script that creates multiple scripts and submits these simultaneously to a computer cluster, and after all of the multiple scripts have completed and the output has been written in the respective folders, I would like to automatically launch another R script that works

[R] Order a data frame based on the order of another data frame

2012-03-05 Thread francesca casalino
Hi, I am trying to match the order of the rownames of a dataframe with the rownames of another dataframe (I can't simply sort both sets because I would have to change the order of many other connected datasets if I did that): Also, the second dataset (snp.matrix$fam) is a snp matrix slot: so for e

[R] Creating data frame with residuals of a data frame

2011-10-24 Thread francesca casalino
Dear experts, I am trying to create a data frame from the residuals I get after having applied a linear regression to each column of a data frame, but I don't know how to create this data frame from the resulting list since the list has differing numbers of rows. So for example: age<- c(5,6,10,14

Re: [R] Mean or mode imputation fro missing values

2011-10-11 Thread francesca casalino
without giving any error…Does anybody have any idea what is going on? Thank you very much for all the great help! -f 2011/10/11 Weidong Gu : > In your case, it may not be sensible to simply fill missing values by > mean or mode as multiple imputation becomes the norm this day. For > your sp

[R] Mean or mode imputation fro missing values

2011-10-11 Thread francesca casalino
Dear R experts, I have a large database made up of mixed data types (numeric, character, factor, ordinal factor) with missing values, and I am looking for a package that would help me impute the missing values using  either the mean if numerical or the mode if character/factor. I maybe could use

Re: [R] Permutation or Bootstrap to obtain p-value for one sample

2011-10-09 Thread francesca casalino
Dear Peter and Tim, Thank you very much for taking the time to explain this to me! It is much more clear now. And sorry for using the space here maybe inappropriately, I really hope this is OK and gets posted, I think it is really important that non-statisticians like myself get a good idea of the

Re: [R] Permutation or Bootstrap to obtain p-value for one sample

2011-10-09 Thread francesca casalino
Thank you very much to both Ken and Peter for the very helpful explanations. Just to understand this better (sorry for repeating but I am also new in statistics…so please correct me where I am wrong): Ken' method: Random sampling of the mean, and then using these means to construct a distribution

Re: [R] Merge two data frames and find common values and non-matching values

2011-10-04 Thread francesca casalino
Sorry---I thought it worked but I think I am actually definitely doing something wrong... The problem might be that there are NA's and there are also duplicated values...My fault. I can't figure out what is going wrong... I'll be more thorough and modify the two df to mirror more what I have to ex

Re: [R] Merge two data frames and find common values and non-matching values

2011-10-04 Thread francesca casalino
Yes, your code did exactly what I needed. Thank you!! -f [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-gu

Re: [R] Import in R with White Spaces

2011-10-04 Thread francesca casalino
Ok I added quoting and it did work...Not sure why, but thank you for both your replies! -f [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide ht