, Michael Dewey ()
escribió:
> Dear Ana
>
> When you installed R it gave you an option to put a shortcut on the
> desktop. If you did that then you can start the R GUI by clicking on it.
> You may need to customise it for your preferences. Specifically it
> starts R with an option cd
Hello,
I am sorry for my ignorance, but what is Rgui.exe?
El mar, 21 nov 2023 a las 12:06, Ivan Krylov ()
escribió:
> В Tue, 21 Nov 2023 10:51:59 +0100
> Ana de las Heras Molina пишет:
>
> > I uninstalled onedrive, I eliminated all the folders and then
> > reinstalled R
Hi,
I uninstalled onedrive, I eliminated all the folders and then reinstalled R
and RStudio... but it is RStudio the one creating a folder
called C:\Users\Ana\OneDrive - Universidad Complutense de Madrid
(UCM)\Documentos.
This is what I obtained with the debugging
Error in setwd(dir) : no es
Hello,
Thank you all for your responses. When I initialize R, the folder in which
it starts is:
getwd()
[1] "C:/Users/Ana/OneDrive - Universidad Complutense de Madrid
(UCM)/Documentos"
So it seems a problem with OneDrive
-Kevin: I do not really understand your question... You mean i
Hello,
I am Ana de las Heras, and I write to you because every time I open RStudio
or R directly I have the following message, before I can do anything at
all:
Error in setwd(dir) : no es posible cambiar el directorio de trabajo
At first I didn't pay much attention to it, but I am having
I have code like this:
data <- read.csv("test1.csv", stringsAsFactors=FALSE, header=TRUE)
# Graph
myplot=ggplot(data, aes(fill=condition, y=value, x=condition)) +
geom_bar(position="dodge", stat="identity", width=0.5) +
scale_fill_manual(values=c("#7b3294", "#c2a5cf", "#a6dba0", "#0088
", "18.7", "19", "6.4", "8.4", "20", "20.6"), PIE.3D.FFT = c("6",
"8.6", "5.7", "4.3", "1 .0", "1 .1", "7.6", "8.4"), PIE.3D.FFT.comm. = c
categorical data, but not continuous, numeric data which are better
> handled with box plots or strip charts.
>
> Do not use printouts of your data since it hides important information.
> Use str(a11) and dput(a11) or dput(head(a11)) to provide useful information
> about your data.
>
end= rownames(data2), beside=
TRUE,las=2,cex.axis=0.7,cex.names=0.7,ylim=c(0,80), col=c("#9e9ac8",
"#6a51a3"))
If I don't log transform my code runs.
Please advise,
Ana
[[alternative HTML version deleted]]
__
R-help@r-pr
with having an open mind is that people keep coming along and
> sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Thu, Oct 14, 2021 at 10:10 AM Ana Marija
> wrote:
>
>> Hi All,
>>
>> I hav
uct of natural log. I have to do
inverse of it.
So for example do do inverse of 0.6931472 I would do:
> 2.718281828^0.6931472
[1] 2
How do I perform this operation for every single value in this data frame?
The original data frame is this dimension:
> dim(b)
[1] 1441 18
T
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
> On Fri, Sep 17, 2021 at 12:22 PM Ana Marija
> wrote:
> >
> > Hi All,
> >
> > I plan to identify metabolite levels that differ between individuals
> > with various retinopa
code given in the tutorial so that
it is matching to mine design.
Can you please help,
Thanks
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
)
> #[1] 145092258
>
> df[14509227,] # beyond nrow(df) by 2
>
>
> Hope this helps,
>
> Rui Barradas
>
>
> Às 15:12 de 16/09/21, Ana Marija escreveu:
> > Hi All,
> >
> > I have lines in file that look like this:
> >
> >> df[1450922
N
"numeric" "integer"
> dim(df)
[1] 145092258
Tried:
> df=na.omit(df)
> dim(df)
[1] 145092258
and:
> library(tidyr)
> d=df %>% drop_na()
> dim(d)
[1] 145092258
Please advise,
Thanks
Ana
___
ec[df1]
oavec<-tausqvec*avec
buw<-sum(oavec*beta1vec)/sum(oavec)
lamvec<-cvec*sum(oavec*(beta1vec-buw)^2)/sigsq
ep<-ep+sum(wcpdf*pf(fcrit,df1,df,lamvec,lower.tail=FALSE))
}
kbpower<-ep/repn
}
kbpower<-kbpowerf()
for minfi package including making MDS plots:
https://www.bioconductor.org/packages/release/workflows/vignettes/methylationArrayAnalysis/inst/doc/methylationArrayAnalysis.html
Thanks
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,
rsplit(names, ":")[-2]]
out <- data[, .(rsid, ref_allele, eff_allele)][,
WGT := files[i]][]
}
return(out)
rm(data)
gc()
}
parallel::stopCluster(cl)
big_data <- rbindlist(lst_out, fill = TRUE)
On Wed, Dec 16, 2020 at 9:31 AM Ana Marija wrote:
>
ster(cl)
big_data <- rbindlist(lst_out)
and i got:
Error in { : task 4 failed - "$ operator is invalid for atomic vectors"
> parallel::stopCluster(cl)
I uploaded 3 RDat file here if you want to try it
https://github.com/montenegrina/sample
Thank you for looking into this
Ana
On Tue, D
- rownames(a)
data <- data.table(names, a["blup"])
nm1 <- c("rsid", "ref_allele", "eff_allele")
any idea how I can rewrite this?
On Tue, Dec 15, 2020 at 8:30 PM Jim Lemon wrote:
>
> Hi Ana,
> I would look at "data" in your second exampl
xing parallelized version.
On Tue, Dec 15, 2020 at 7:35 PM Jim Lemon wrote:
>
> Hi Ana,
> My guess is that in your second code fragment you are assigning the
> rownames of "a" and the _values_ contained in a$blup to the data.table
> "data". As I don't have
quot;, "eff_allele")
data[, (nm1) := tstrsplit(names, ":")[-2]]
return(data[, .(rsid, weight = blup, ref_allele, eff_allele)][,
WGT := files[i]][])
}
parallel::stopCluster(cl)
big_data <- rbindlist(lst_out)
I am getting this Error:
Error in { : t
lp page, ?corrplot .
>
> --
>
> David.
>
> On 11/6/20 6:08 AM, Ana Marija wrote:
>
> sorry forgot to attach the plot.
>
> On Fri, Nov 6, 2020 at 8:07 AM Ana Marija wrote:
>
> Hello
>
> I have data like this:
>
> head(my_data)
>
> subjects DIABDUR HB
sorry forgot to attach the plot.
On Fri, Nov 6, 2020 at 8:07 AM Ana Marija wrote:
>
> Hello
>
> I have data like this:
>
> > head(my_data)
> subjects DIABDUR HBA1C ESRD SEX AGE PHENO C1 C2
> 1 fam0110_G110 38 9.41 2 51 2 -
uot;correlation.pdf")
corrplot(res, type = "upper", order = "hclust",
tl.col = "black", tl.srt = 45)
dev.off()
and I am getting the plot in attach. How to make it so that my
variables are shown on the plot in the order they are in my_data data
frame?
Th
Makes sense, thank you!
On Wed, 21 Oct 2020 at 17:46, Rolf Turner wrote:
>
> On Wed, 21 Oct 2020 16:15:22 -0500
> Ana Marija wrote:
>
> > Hello,
> >
> > I have a data frame with one column:
> >
> > > remove
> >
> >
lename %in% as.character(remove$V1)]
> >
> >
> > Hope this helps,
> >
> > Rui Barradas
> >
> > Às 22:15 de 21/10/20, Ana Marija escreveu:
> >> Hello,
> >>
> >> I have a data frame with one column:
> >&g
ied doing this:
b= celFiles[!basename(celFiles) %in% as.character(remove$V1)]
but none of the 8th entries in "remove" data frame have been removed.
Please advise,
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https:
a bad idea.
> And if it's a good idea, then how much to trim.
>
>
> On Sat, Oct 10, 2020 at 5:47 AM Ana Marija
> wrote:
> >
> > Hi Bert,
> >
> > Another confrontational response from you...
> >
> > You might have noticed that I use the wor
erpret the plot. I was hoping someone here knows
that and can help me.
Ana
On Fri, Oct 9, 2020 at 11:31 AM Bert Gunter wrote:
>
> I recommend that you consult with a local statistical expert. Much of what
> you say (outliers?!?) seems to make little sense, and your statistical
> k
var sd
FQC.10090295 0.0327 0.002678 0.0517
FQC.10119363 0.0220 0.000978 0.0313
FQC.10132112 0.0275 0.002088 0.0457
FQC.10201128 0.0169 0.000289 0.0170
FQC.10208432 0.0443 0.004081 0.0639
FQC.10218466 0.0116 0.000131 0.0115
...
the distribution is not normal, it is right-skewed.
Cheers,
Ana
NP[SNP$density>400,]
and plot it again:
p <- ggplot(a, mapping = aes(x = mean, y = var))
p <- p + geom_density_2d() + geom_point() + my.theme + ggtitle("SNPS_red")
On Thu, Oct 8, 2020 at 3:52 PM Ana Marija wrote:
>
> Hello,
>
> I have a data frame like this:
>
>
e ellipses-contours?
2. how do I extract from my data frame the points which are outside of ellipses?
Thanks
Ana
snps.pdf
Description: Adobe PDF document
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mail
4:ncol(mt)) mt[,i] <- 1 + (names(mt)[i]== mt$PLATE)
Thanks!
On Tue, Sep 29, 2020 at 12:08 PM Ana Marija wrote:
>
> HI Bert,
>
> thank you for getting back to me.
> I tried this:
>
> > dat <- cbind(mc, matrix(0,ncol = 34))
> > head(dat)
> FID IID PLATE 1
A NA NA NA NA NA NA NA NA
...
so this gives me the correct columns. Now is the question of how to
replace NA with 2 id column name matches the rownname in PLATE column
with 2 otherwise it is 1.
Cheers,
Ana
On Tue, Sep 29, 2020 at 11:46 AM Bert Gunter wrote:
>
value.var to override.
Please advise,
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, m
t RHS position 1 taken as TRUE when
assigning to type 'logical' (column 6 named 'PHENO')
Please advise,
Ana
On Wed, Sep 23, 2020 at 2:48 PM Jeremie Juste wrote:
>
>
> Hello Ana Marija,
>
> I cannot reproduce your error,
>
> with a$PHENO=ifelse(a$PLASER==2 |a$
, 2020 at 11:43 AM Ana Marija wrote:
>
> Hello,
>
> I have a data frame as shown bellow.
> I want to create a new column PHENO which will be defined as follows:
> if CURRELIG==1 -> PHENO==1
> in the above subset those that have:
> PLASER==2 -> PHENO==2
> and
>
1 1 2
fam5709 G57091 1 1
fam5713 G57131 1 1
fam5715 G57151 1 1
fam5718 G57181 1 1
Please advise,
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,
HI Jim,
fantastic solution!
Thank you so much!!!
Ana
On Thu, Sep 17, 2020 at 6:01 PM Jim Lemon wrote:
>
> Hi Ana,
> Sorry it's not in ggplot, but it may help:
>
> d<-read.table(text="CHR counts name
> 1 193554 old
> 2 220816 old
> 3 174350 old
>
show counts per CHR in "old" and "new" side by side.
Please advise,
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.
% had Type 1. (my TD covariate is reference for
the type of diabetes) In the attach is the description of the data.
Cheers,
Ana
On Tue, Sep 15, 2020 at 7:59 PM David Winsemius wrote:
>
>
> On 9/15/20 8:57 AM, Ana Marija wrote:
> > Hi Abby and David,
> >
> > Thanks
uge effect or what
kind of conclusion I can draw from that?
Thanks for all your help
Ana
On Tue, Sep 15, 2020 at 1:26 AM David Winsemius wrote:
>
> There is a user-group for PLINK, easily found by looking at the page you
> cited. This is not the correct place to submit such questions
sorry not replace with NA but with empty string for a name, for example
for example this:
> geneSymbol["Ku8QhfS0n_hIOABXuE"]
Ku8QhfS0n_hIOABXuE
"MACC1"
would go when I subject it to
> geneSymbol["Ku8QhfS0n_hIOABXuE"]
Ku8QhfS0n_hIOABXuE
On Mon,
;,"odF3XHR8CVl4SAUaUQ")]
> geneSymbol[c("0lQ1XozriVZTn.PezY","uaeFiCdegrnWFijF_s","ZOluqaxSe3ndekoNng","912ny6eCHjnlY2XSCU","odF3XHR8CVl4SAUaUQ")]
0lQ1XozriVZTn.PezY uaeFiCdegrnWFijF_s ZOluqaxSe3ndekoNng
912ny6eCHjnlY2XSCU odF3XHR8CVl4SAUaUQ
better way to do this from the data I already have?
Thanks
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.h
ed for tutorials also. This list cannot
> substitute for such homework on your own.
>
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along and
> sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County&qu
It seems that "treatment" and "patient" are just vectors.
> treatment
[1] "treat" "treat" "treat" "treat" "treat" "control" "control"
[8] "control" "control" "cont
1] "treatment"
> str(treatment)
chr [1:10] "treat" "treat" "treat" "treat" "treat" "control" "control" ...
but this is not the format I need.
On Fri, Jul 31, 2020 at 9:18 PM Ana Marija wrote:
>
> Hello,
>
Hello,
I have this file:
> a=load("paired_example.Rdata")
> a
[1] "rawdata" "treatment" "patient"
I can extract "rawdata" with:
dat<-local(get(load("paired_example.Rdata")))
Can you please advise how
Thank you so much as.character(r) indeed resolved the issue!
On Sat, Jun 20, 2020 at 3:47 AM Ivan Krylov wrote:
>
> On Fri, 19 Jun 2020 19:36:41 -0500
> Ana Marija wrote:
>
> > Error in cat(x, file = file, sep = c(rep.int(sep, ncolumns - 1),
> > "\n"), :
s a foreach package but I try to avoid extras
>
> make your for statements:
>
> for ( a in rownames(f1) ) {
>
> # a will now be a row number rather than the value, so replace ' a ' in
> the paste0 with: f1[ a, 1]
>
> so
>
> ext <- paste0( "/ld/human
SNPs at the time and I need to do that for 300 pairs
On Fri, Jun 19, 2020 at 6:42 PM Rasmus Liland wrote:
>
> On 2020-06-19 14:34 -0500, Ana Marija wrote:
> >
> > I have two files (each has 300 lines)like this:
>
> The example looks quite similar to the R example in
&g
) You may have something that adds it but better to use
> something that works. So try using:
>
> library(readr)
> f1 <- read_tsv("1g.txt", col.names=F)
>
> This will give you a tibble with f1$X1 with the file in it
>
> then loop it with (a in as.list(f1[,1])
&g
ot open compressed file '1000GENOMES:phase_3:KHV.rds', probable
reason 'No such file or directory'
> lapply(r[1:4], function(x) {
+ jsonlite::fromJSON(jsonlite::toJSON(httr::content(x)))
+ })
Error in lapply(r[1:4], function(x) { : object 'r'
uot;)
>
> r <- GET(paste(server, ext, sep = ""),
> content_type("application/json"))
>
> # You presumably need to do something with 'r' at the
> moment its over written by the next loop.. were
> #
and output every line of result in list.txt
The process is illustrated in the attachment.
Please help,
Ana
lists.pdf
Description: Adobe PDF document
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/list
te in the presence of ties
as I understand high p-values here say I cannot claim statistical
support for a difference, but low p-values are not evidence of
sameness?
D should be the maximum difference between the two probability distributions?
On Wed, Jun 17, 2020 at 3:06 PM Ana Marija wrote:
>
&g
99937, num df = 1454159, denom df = 1454159, p-value = 0.7057
alternative hypothesis: true ratio of variances is not equal to 1
95 percent confidence interval:
0.9970808 1.0016750
sample estimates:
ratio of variances
0.9993739
Or some other test makes more sense?
Thanks
Ana
qq-plots-c
np="MARKER",suggestiveline
= F, genomewideline = F,main = "Gokind old",col = c("red2", "orange",
"gold1","springgreen4", "dodgerblue",
"darkorchid1","orchid1"),ylim=c(0,9),cex = 0.6)
in this qqman function.
C
og,chr="CHR",bp="POS",p="META_pval",snp="MARKER",ylim
= c(0, 10))
dev.off()
Thanks
Ana
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read th
HI Jim
thank you so much! This is amazing answer!!!
Ana
On Sat, Jun 13, 2020 at 4:09 AM Jim Lemon wrote:
>
> Right, back from shopping. Since you have fourteen rows containing NAs
> and you only want seven, we can infer that half of them must go. As
> they are neatly divided into s
iate logic to get rid of only the
> ones you don't want.
>
> Jim
>
> On Sat, Jun 13, 2020 at 12:50 PM Ana Marija
> wrote:
> >
> > Hi Rasmus,
> >
> > thank you for getting back to be, the command your provided seems to
> > add all 11 NAs to 2s
&
; > > On Sat, Jun 13, 2020 at 10:46 AM Ana Marija wrote:
> > > >
> > > > I am trying to make a new column
> > > > "pheno" so that I reduce the number
> > > > of NAs
> > >
> > > it looks like those two NA values in
>
exclude = NULL)
12
859 828 11
Am I am doing something wrong?
Thanks
Ana
On Fri, Jun 12, 2020 at 8:06 PM Jim Lemon wrote:
>
> Hi Ana,
> From your desired result, it looks like those two NA values in PLASER
> are the ones you want to drop.
> If so, try this:
>
> b
,1)
> table(b$pheno, exclude = NULL)
12
859 828 11
I would like to reduce this number of NAs to be 7
so I would like to have in "pheno column"
7 NAs
825 2s (825=691+14+70+7+43)
and the rest would be 1s (866=1698-7-825)
How can I change the above command t
e shown on the x-axis?
On Thu, Jun 11, 2020 at 4:39 PM wrote:
>
> Your dots are too big!
>
> Add
>
> geom_points(... , size = 1
>
> May need to play... 0.5 or 0.1?
>
> On 11 Jun 2020 22:26, Ana Marija wrote:
>
> I tried it,
> ggplot( data = tmp.tidy) +geom
I tried it,
ggplot( data = tmp.tidy) +geom_point( aes(y = BP,x = CHR,color=key)
,position = "jitter" )
I got the attached
On Thu, Jun 11, 2020 at 4:18 PM wrote:
>
> Try adding
> position = "jitter" to the geom_point(...
>
>
>
> On 11 Jun 2020 21:41, Ana M
Hello,
I tried your code and this is what I got
I really need two groups side by side shown per chromosome as it is here:
https://imgur.com/a/pj40c
on the image there are 4 groups I do have only two
On Thu, Jun 11, 2020 at 11:52 AM wrote:
>
> On 2020-06-11 15:59, Ana Marija wrote:
>
me, and CHR would go from 1 to 22
On Thu, Jun 11, 2020 at 9:26 AM wrote:
>
> On 2020-06-11 14:54, Ana Marija wrote:
> > Hello,
> >
> > I expected it to look like this:
> > https://imgur.com/a/pj40c
> >
>
> Ah - so all on the one plot? - so you don't w
value: num 0.952 0.475 0.529 0.255 0.295 ...
Unfortunately qqman doesn't do this kind of overlay of two plots
On Wed, Jun 10, 2020 at 11:24 PM John wrote:
>
> On Wed, 10 Jun 2020 15:36:11 -0500
> Ana Marija wrote:
>
> > Hello,
> >
> > I have a data frame lik
ssoc.logistic.jpeg")
manhattan(results_log,chr="CHR",bp="BP",p="P",snp="SNP", main =
"Manhattan plot:logistic")
dev.off()
this code does work with another data set, the jpeg() does work. I
will try now with PNG
On Wed, Jun 10, 2020 at 5:17 PM Jim Le
uot;SNP", main =
"Manhattan plot: logistic")
dev.off()
but I am getting:
Error in plot.window(...) : need finite 'ylim' values
Calls: manhattan ... do.call -> plot -> plot.default -> localWindow ->
plot.window
Execution halted
Please advise,
Thanks
Ana
Hi Rui,
thank you so much, that is exactly what I needed!
Cheers,
Ana
On Wed, Jun 3, 2020 at 11:50 AM Rui Barradas wrote:
>
> Hello,
>
> If you want to filter out rows with any of the values in a 'unwanted'
> vector, try the following.
>
> First, create a te
created just remove any row that contains any of those variables.
Thanks
Ana
On Wed, 3 Jun 2020 at 11:00, Bert Gunter wrote:
> I suggest that you forget all that fancy stuff (and this is not a use
> case for regular expressions).
> Use %in% with logical subscripting instead -- basic R func
grep('^E11', x, value = TRUE))
d1=unlist(s1)
d11=unique(d1)
> d11
[1] "E11" "E119" "E113" "E115" "E111" "E114" "E110" "E118" "E116" "E112"
[11] "E117"
I need help with c
Hi Jim,
not in this case, but thanks for asking!
Ana
On Mon, Jun 1, 2020 at 10:04 PM Jim Lemon wrote:
>
> So recombination sticks out its foot before us. Do you want to account
> for gene linkage?
>
> JIm
>
> On Tue, Jun 2, 2020 at 11:55 AM Ana Marija
> wrote:
>
r","Chr"))
nn2<-merge(nn1,ret1,by=c("Marker","Chr"))
> Marker3<-nn2$Marker
> length(Marker3)
[1] 3742962
> Marker4<-nn1$Marker
> length(Marker4)
[1] 373
On Mon, Jun 1, 2020 at 8:50 PM Ana Marija wrote:
>
> Hi David,
>
> that is a
43494 9
How would I rewrite this code so that is merging by Chr and Marker
column? It seems that a Marker can be under a few Chr.
On Mon, Jun 1, 2020 at 8:41 PM David Winsemius wrote:
>
>
> On 6/1/20 5:40 PM, Ana Marija wrote:
> > Hi Jim,
> >
> > thank you so much
.
Cheers
Ana
On Mon, Jun 1, 2020 at 7:31 PM Jim Lemon wrote:
> Hi Ana,
> Not too hard, but your example has all the "marker" fields in common.
> So using a sample that will show the expected result:
>
> neu1<-read.table(text="Chr BP Marker MAF A1 A2 Directio
991:G:A 0.0982587 A G - 0.272936
1608
2234642 1 14726 1:14726:G:A 0.1340170 A G - 0.594538
1608
Is there is a way to create another 3 data frames, say neu2, nep2, ret2
which would only contain lines that have the same entries in Marker column
for all 3 data frames?
Thanks
Ana
4156 G A + 0.484813 1641
2 10 10645 10:10645:A:C 0.216027 C A + 0.73597 1641
Can you please tell me what colClasses=colClassvec suppose to do?
Thanks
Ana
On Mon, Jun 1, 2020 at 4:13 PM David Winsemius
wrote:
>
> On 6/1/20 1:37 PM, Ana Marija wrote:
&g
coming along and
> sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Mon, Jun 1, 2020 at 1:38 PM Ana Marija
> wrote:
>
>> Hello,
>>
>> I have a dataframe like this:
>>
>> Chr
ble("gokind.neuropathy.fin", header=T,stringsAsFactors=FALSE)
and every column is numeric. How to say have all numeric ones stay numeric
like: Chr, BP, MAF, pValue, N
Thanks
Ana
[[alternative HTML version deleted]]
__
R-help@r-projec
ous(breaks
= unique(d))
ed
where:
> head(e)
ESRD pheno
11 1
21 1
31 2
41 1
51 1
> sapply(e,class)
ESRD pheno
"integer" "factor"
On Fri, May 22, 2020 at 1:52 PM Ana Marija wrote:
>
> Hello,
>
> I made t
Hello,
I made the plot in attach via:
ed<-ggplot(e) +
geom_bar(aes(x = ESRD, fill =
factor(pheno,labels=c("control","case"+scale_fill_manual(values=c("#56B4E9","#E7B800"))+labs(fill="pheno")
amp;consent=&temp=1
On Fri, May 22, 2020 at 12:18 AM Eric Berger wrote:
>
> Hi Ana,
> This is a very common question about ggplot.
> A quick search turns up lots of hits that answer your question. Here
> are a couple
> https://community.rstudio.com/t/trouble-scaling-
ou the whole dataset if you would like to try with it
On Thu, May 21, 2020 at 11:14 PM Jim Lemon wrote:
>
> Hi Ana,
> Just noticed a typo from a hasty cut-paste. Two lines should read:
>
> casehist<-table(cut(aafd$HBAIC[aafd$pheno=="case"],breaks=0:15))
> cont
the result would basically look something like this on in attach or
the overlay of those two plots
On Thu, May 21, 2020 at 5:23 PM Ana Marija wrote:
>
> Hello,
>
> I have a data frame like this:
> > head(a)
> FID IID FLASER PLASER DIABDUR HBA1C ESRD pheno
dth=.5,
position="dodge")
there is 848 who have "case" in pheno column and 892 who have
"control" in pheno column.
I would like to have on y-axis shown percentage of individuals which
have either "case" or "control" in pheno instead of count.
Please advis
5,annotatePval = 0.0001)
but I am getting this error:
Error in textxy(topSNPs$pos, -log10(topSNPs$P), offset = 0.625, labs =
topSNPs$SNP, :
formal argument "cex" matched by multiple actual arguments
Do you by any chance know how to do this?
Cheers
Ana
On Wed, May 20, 2020 at 4:1
Hi Michael,
can you please send me code how that would be done?
Thanks
Ana
On Tue, May 19, 2020 at 11:18 AM Michael Dewey wrote:
>
> Dear Ana
>
> Perhaps paste together SNP and GENE using paste() and then supply that
> as the snp parameter.
>
> Michael
>
> On 19/05/2
SNP=="rs4081570",]
SNPP CHR BP GENE
1 rs4081570 6.564447e-05 19 15918647 UCA1
> a[a$SNP=="rs11867934",]
SNPP CHR BP GENE
1021 rs11867934 6.738066e-06 17 16933404 FLCN
Right now my plot only has SNP name for those 2
Hi Rui,
thank you so much that is exactly what I needed!
Cheers,
Ana
On Fri, May 15, 2020 at 5:12 PM Rui Barradas wrote:
>
> Hello,
>
> I have tried several options and with large dataframes this one was the
> fastest (in my tests, of the ones I have tried).
>
>
> s1 &
ch start with E10.
Thanks
Ana
On Fri, May 15, 2020 at 12:13 PM Jeff Newmiller
wrote:
>
> Read about regular expressions... they are extremely useful.
>
> df1 <- tot %>% filter_all(any_vars(grepl( '^E10', .)))
>
> It is bad form not to put spaces around the <
of strings that start with E10...
it would look something like this:
[1] "E102" "E109" "E108" "E103" "E104" "E105" "E101" "E106" "E107"
Thanks
Ana
__
R-help@r
uot;Populations in millions (log scale)")+
+ylab("Total number of murders (log scale)")+
+ggtitle("US Gun Murders in US 2010")+
+scale_color_discrete(name="Region")+
+theme_economist()
Is this an issue with my d
>
> t.test(x, mu = beta0)$statistic
>
>
> But in this case the estimator is the estimator for the mean value.
>
> Hope this helps,
>
> Rui Barradas
>
> Às 15:54 de 06/05/20, Ana Marija escreveu:
> > Thanks Patrick, so in conclusion this is fine?
> > z-sco
Thanks Patrick, so in conclusion this is fine?
z-score=Beta/StdErr
On Wed, May 6, 2020 at 9:52 AM Patrick (Malone Quantitative)
wrote:
>
> Guessing for Ana, but no, that's a different meaning. Beta/StdErr is a
> z statistic--a test statistic against (usually) the tails of the
gt; Hope this helps,
>
> Rui Barradas
>
> Às 15:31 de 06/05/20, Ana Marija escreveu:
> > Hi Rui,
> >
> > Thank you for getting back to me. Is there is a better way to
> > calculate Z scores if I have p values, SE and Beta?
> >
> > Thanks
> &g
1 - 100 of 342 matches
Mail list logo