search for: jinghuazhao

Displaying 20 results from an estimated 22 matches for "jinghuazhao".

2019 May 31
2
use of buffers in sprintf and snprintf
..., sizeof(tempname), "%s.%d", of1name, j) >= sizeof(tempname)) Rf_error("file name is too long"); BTW: most OSes systems have a path limits that are no lower than 256 so you should allow at least as much. Cheers, Simon > On May 29, 2019, at 11:49 AM, jing hua zhao <jinghuazhao at hotmail.com> wrote: > > Dear R-developers, > > I am struggling with packaging with sprintf and snprintf() as the following WARNINGS from gcc 9.x, > > hap_c.c:380:46: warning: ?%d? directive output may be truncated writing between 1 and 10 bytes into a region of size betw...
2024 Apr 16
1
read.csv
...i.nlm.nih.gov/15214961/ https://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-1044-7 https://www.nature.com/articles/d41586-021-02211-4 https://www.theverge.com/2020/8/6/21355674/human-genes-rename-microsoft-excel-misreading-dates On Tue, Apr 16, 2024 at 3:46?AM jing hua zhao <jinghuazhao at hotmail.com> wrote: > > Dear R-developers, > > I came to a somewhat unexpected behaviour of read.csv() which is trivial but worthwhile to note -- my data involves a protein named "1433E" but to save space I drop the quote so it becomes, > > Gene,SNP,prot,log10p &gt...
2019 Jun 21
4
Calculation of e^{z^2/2} for a normal deviate z
...point representation cuts out at 1e+/-308 or so. If you want to go outside that range, you may need explicit manipulation of the log values. qnorm() itself seems quite happy with much smaller values: > qnorm(-5000, log=TRUE) [1] -99.94475 -pd > On 21 Jun 2019, at 17:11 , jing hua zhao <jinghuazhao at hotmail.com> wrote: > > Dear Rui, > > Thanks for your quick reply -- this allows me to see the bottom of this. I was hoping we could have a handle of those p in genmoics such as 1e-300 or smaller. > > Best wishes, > > > Jing Hua > > _____________________...
2010 Jul 09
4
Mysterious behavior
I had trouble with some tests for the survival suite last night that I cannot explain. Framework: Ubuntu Linux, R2.11. For testing survival I have a separate directory and Makefile. I pull everything into the local .RData, no packages, library, or namespace. (It's easier to add test modifications to a routine in a chain of calls). A test of survreg + psline would fail because
2018 May 28
5
readLines function with R >= 3.5.0
On 28.05.2018 11:07, G?bor Cs?rdi wrote: > stdin() is not the same as file("stdin"), see the note in ?stdin. In particular stdin() works in an interactive session but not when R -f / Rscript is used, since it does not wait for the user to input anything: $ R -f readLines.R R version 3.5.0 (2018-04-23) -- "Joy in Playing" Copyright (C) 2018 The R Foundation for Statistical
2024 Apr 16
5
read.csv
Dear R-developers, I came to a somewhat unexpected behaviour of read.csv() which is trivial but worthwhile to note -- my data involves a protein named "1433E" but to save space I drop the quote so it becomes, Gene,SNP,prot,log10p YWHAE,13:62129097_C_T,1433E,7.35 YWHAE,4:72617557_T_TA,1433E,7.73 Both read.cv() and readr::read_csv() consider prot(ein) name as (possibly confused by
2019 May 31
0
use of buffers in sprintf and snprintf
...the absolute pathname that counts." Admittedly, it's been several years when I last looked into it, but at the time when I wrote that I spent lots of time investigating it. /Henrik > > Cheers, > Simon > > > > > > On May 29, 2019, at 11:49 AM, jing hua zhao <jinghuazhao at hotmail.com> wrote: > > > > Dear R-developers, > > > > I am struggling with packaging with sprintf and snprintf() as the following WARNINGS from gcc 9.x, > > > > hap_c.c:380:46: warning: ?%d? directive output may be truncated writing between 1 and 10 bytes...
2020 Nov 10
0
Help on mapping memory
...hing is that it works fine with the documentation data and even with this data there were times it could be tweaked to work (therefore PAN.txt, test.rda and test.log were as intended there). I have extracted the pan.f, pan.R from the package and leave the Bash/R scripts all here, https://github.com/jinghuazhao/R/tree/master/tests, short of adding a driver program to pan.f and debug without R but before doing that any idea/insight would be greatly appreciated. Thank you so much, Jing Hua Zhao [[alternative HTML version deleted]]
2019 Jun 23
2
Calculation of e^{z^2/2} for a normal deviate z
...1e+/-308 or so. If you want to go outside that range, you may need explicit manipulation of the log values. qnorm() itself seems quite happy with much smaller values: > >> qnorm(-5000, log=TRUE) > [1] -99.94475 > > -pd > >> On 21 Jun 2019, at 17:11 , jing hua zhao <jinghuazhao at hotmail.com> wrote: >> >> Dear Rui, >> >> Thanks for your quick reply -- this allows me to see the bottom of this. I was hoping we could have a handle of those p in genmoics such as 1e-300 or smaller. >> >> Best wishes, >> >> >> Jing Hua &...
2024 Apr 16
1
read.csv
...nt character (.), optionally followed by an exponent part (an e or E character followed by an optional sign and a sequence of digits)." [Of course, there would be nothing to stop e.g. "1433E1" from being converted to numeric.] -pd > On 16 Apr 2024, at 12:46 , jing hua zhao <jinghuazhao at hotmail.com> wrote: > > Dear R-developers, > > I came to a somewhat unexpected behaviour of read.csv() which is trivial but worthwhile to note -- my data involves a protein named "1433E" but to save space I drop the quote so it becomes, > > Gene,SNP,prot,log10p...
2019 Jun 21
4
Calculation of e^{z^2/2} for a normal deviate z
Hello, Well, try it: p <- .Machine$double.eps^seq(0.5, 1, by = 0.05) z <- qnorm(p/2) pnorm(z) # [1] 7.450581e-09 1.228888e-09 2.026908e-10 3.343152e-11 5.514145e-12 # [6] 9.094947e-13 1.500107e-13 2.474254e-14 4.080996e-15 6.731134e-16 #[11] 1.110223e-16 p/2 # [1] 7.450581e-09 1.228888e-09 2.026908e-10 3.343152e-11 5.514145e-12 # [6] 9.094947e-13 1.500107e-13 2.474254e-14 4.080996e-15
2019 May 30
2
use of buffers in sprintf and snprintf
...NINGS were initially given with https://www.stats.ox.ac.uk/pub/bdr/gcc9/README.txt and those specification might as well used in --configure-args above. Besst regards, Jing Hua ________________________________ From: R-devel <r-devel-bounces at r-project.org> on behalf of jing hua zhao <jinghuazhao at hotmail.com> Sent: 29 May 2019 15:49 To: r-devel at r-project.org Subject: [Rd] use of buffers in sprintf and snprintf Dear R-developers, I am struggling with packaging with sprintf and snprintf() as the following WARNINGS from gcc 9.x, hap_c.c:380:46: warning: ?%d? directive output may...
2019 Jun 21
0
Calculation of e^{z^2/2} for a normal deviate z
...point representation cuts out at 1e+/-308 or so. If you want to go outside that range, you may need explicit manipulation of the log values. qnorm() itself seems quite happy with much smaller values: > qnorm(-5000, log=TRUE) [1] -99.94475 -pd > On 21 Jun 2019, at 17:11 , jing hua zhao <jinghuazhao at hotmail.com> wrote: > > Dear Rui, > > Thanks for your quick reply -- this allows me to see the bottom of this. I was hoping we could have a handle of those p in genmoics such as 1e-300 or smaller. > > Best wishes, > > > Jing Hua > > ___________________________...
2008 Feb 07
5
pnorm
Dear R list, I calculated a two-sided p values according to 2*(1-pnorm(8.104474)), which gives 4.440892e-16. However, it appears to be 5.30E-16 by a colleague and 5.2974E-16 from SAS. I tried to get around with mvtnorm package but it turns out to be using pnorm for univariate case. I should have missed some earlier discussions, but for the moment is there any short answer for a higher
2019 Jun 23
0
Calculation of e^{z^2/2} for a normal deviate z
...tside that range, you may need explicit manipulation of > the log values. qnorm() itself seems quite happy with much smaller values: > > > >> qnorm(-5000, log=TRUE) > > [1] -99.94475 > > > > -pd > > > >> On 21 Jun 2019, at 17:11 , jing hua zhao <jinghuazhao at hotmail.com> > wrote: > >> > >> Dear Rui, > >> > >> Thanks for your quick reply -- this allows me to see the bottom of > this. I was hoping we could have a handle of those p in genmoics such as > 1e-300 or smaller. > >> > >> Best...
2018 Jul 08
2
[R] consider running tools::compactPDF(gs_quality = "ebook")
Dear all, I run R CMD build --compact-vignettes="both" gRbase and/or R CMD build --compact-vignettes="gs+qpdf" gRbase and in the log from r-devel (on winbuilder) I get * checking sizes of PDF files under 'inst/doc' ... WARNING 'gs+qpdf' made some significant size reductions: compacted 'gRbase-arrays.pdf' from 421Kb to 115Kb consider
2019 Jun 24
2
Calculation of e^{z^2/2} for a normal deviate z
...> the log values. qnorm() itself seems quite happy with much smaller values: >> > >> >> qnorm(-5000, log=TRUE) >> > [1] -99.94475 >> > >> > -pd >> > >> >> On 21 Jun 2019, at 17:11 , jing hua zhao <jinghuazhao at hotmail.com> >> wrote: >> >> >> >> Dear Rui, >> >> >> >> Thanks for your quick reply -- this allows me to see the bottom of >> this. I was hoping we could have a handle of those p in genmoics such as >&...
2019 Jun 24
0
Calculation of e^{z^2/2} for a normal deviate z
...> the log values. qnorm() itself seems quite happy with much smaller values: >> > >> >> qnorm(-5000, log=TRUE) >> > [1] -99.94475 >> > >> > -pd >> > >> >> On 21 Jun 2019, at 17:11 , jing hua zhao <jinghuazhao at hotmail.com> >> wrote: >> >> >> >> Dear Rui, >> >> >> >> Thanks for your quick reply -- this allows me to see the bottom of >> this. I was hoping we could have a handle of those p in genmoics such as >&...
2017 Apr 20
1
Intel MKL compiling issue
On our Scientific Linux 6, there is gcc --version gcc (GCC) 4.4.7 20120313 (Red Hat 4.4.7-17) and later version (a symbolic at HOME to the system directory), export MKL=/home/jhz22/11.3.3.210/mkl export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$MKLROOT/lib/intel64 ./configure --prefix=/home/jhz22 LDFLAGS=-L/home/jhz22/lib CPPFLAGS=-I/home/jhz22/include \ --enable-R-shlib --with-lapack \
2017 Jun 04
2
read.table
Hi All, I wonder if there should be one character for quote= in read.table, i.e., > args(read.table) function (file, header = FALSE, sep = "", quote = "\"'", dec = ".", ... I have a file containing the following lines, 08248-GOTERM 3'-phosphoadenosine 5'-phosphosulfate biosynthetic process 08279-GOTERM 3'-phosphoadenosine