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

2014-05-27 Thread francesca casalino
May 2014, francesca casalino wrote: I am having this problem because I need to run a meta-analysis and to align all the variants between the different studies included in the meta-analysis I need to know the effect allele used to get the beta (so that I can flip the beta if the effect allele

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 david.du...@qimr.edu.au: francesca casalino francy.casal...@gmail.com asked: Does anyone know how to find the reference allele used for genetic

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 francy.casal...@gmail.com: 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 david.du...@qimr.edu.au: francesca casalino

[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, but

[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] ggplot2 and facet_wrap help

2013-02-18 Thread francesca casalino
, Feb 18, 2013 at 6:04 AM, francesca casalino francy.casal...@gmail.com wrote: 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

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 francy.casal...@gmail.com: Dear R experts, I have the logarithms of 2 values: log

[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

[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

[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-

[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] Mean or mode imputation fro missing values

2011-10-11 Thread francesca casalino
not be sensible to simply fill missing values by mean or mode as multiple imputation becomes the norm this day. For your specific question, na.roughfix in randomForest package would do the work. Weidong Gu On Tue, Oct 11, 2011 at 8:11 AM, francesca casalino francy.casal...@gmail.com wrote

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

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

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

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

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