similar to: numbers as part of long character

hello all sorry for the following "none" R related question. does anyone know of a reference to calculate the following identity: |I + ABC| where I is an identity matrix and A, B,C may not have to be square matrices? you help will be greatly appreciated. H. V. Henderson; S. R. Searle SIAM Review, Vol. 23, No. 1. (Jan., 1981), pp. 53-60. provides a result to

I have a data frame named "database" with panel data, a little piece of which looks like this: Symbol Name Trial Factor1 Factor2 External 1 548140 A 1 -3.87 -0.32 0.01 2 547400 B 1 12.11 -0.68 0.40 3 547173 C 1

Since we're on the topic of book reviews, I just received Phil Spector's new R book called "Data Manipulation with R" and it is also quite a nice book. I haven't gone through it all and I won't give a detailed review but I have gotten a lot out of the first 100 pages that I have read. Note that I've been using R for almost 1.5 years so , for me, it's a

Dear R-Help, I'm new to R and struggling with weighting data when I run regression. I've tried to use search to solve my problem but haven't found anything helpful so far. I (successfully) import data from SPSS (15) and try to run a linear regression on a subset of my data file where WEIGHT is the name of my weighting variable (numeric), e.g.: library(foreign)

After many hours of debugging code, I came to the conclusion that I have a fundamental misunderstanding regarding eval, and hope that someone here can explain to me, why the following code acts as it does: foo <- function(expr) { eval(substitute(expr), envir=list(a=5), enclos=parent.frame()) } bar <- function(er) { foo(er) } > foo(a) [1] 5 > bar(a) Error in eval(expr,

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

You may want to look into using the log option to qnorm e.g., in round figures: > log(1e-300) [1] -690.7755 > qnorm(-691, log=TRUE) [1] -37.05315 > exp(37^2/2) [1] 1.881797e+297 > exp(-37^2/2) [1] 5.314068e-298 Notice that floating 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()

ok,thanks. Suh HUA SU 邮箱：suhua.tanke@gmail.com 签名由 网易邮箱大师 定制 On 09/11/2019 02:59, Cole Robinson wrote: On 8/25/19 10:21 PM, Su Hua wrote: > Hi, everyone, ask a question, which version can fully support the device type > of qemu hw/block/vhost-user-blk.c? If so, what should the format of the xml > file look like? > Hi, libvirt does not currently support vhost-user-blk - Cole

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

Hi All, I'm wondering whether there is a quick way to know the position of a specific charcater in a long character: for example frg="((D:41.04,I:41.04):45.05,(((E:2.32,((G:0.67,J:0.67):0.44,H:1.11):1.21)" and I would like to know that the 1st, 2nd, 26th, 27th, 28th ... character is "(", is there a quick way to do that? Thanks! Hua

I agree with many the sentiments about the wisdom of computing very small p-values (although the example below may win some kind of a prize: I've seen people talking about p-values of the order of 10^(-2000), but never 10^(-(10^8)) !). That said, there are a several tricks for getting more reasonable sums of very small probabilities. The first is to scale the p-values by dividing the

> >> I tried to ping the server, and it can accept all data. Is there a good way >> to check it? > > ssh -v, ssh -vv and ssh -vvv might give you some interesting information. > Yes, I tried ssh -vvv. It gave a lot of information while login, but it quit without any further information except for ?write failed, broken pipe?. >> The problem came out while

Hi All, I have a problem of putting long titles on a graph: for example, x= seq(1:100) y=seq(1:100) plot(x,y,main="p=0.05:A-B=3,C-D=10,D-E=100,A-F=2,AFR-E=3,ACE-D=1,ADEF-M=0,AED-E=10,DE-F=3,AB-J=4,AC-J=10,ED-F=1,ED-B=4,AF-B=10,CD-S=10,AM-C=4") R seems not able to print the whole title. The title content might be changing and thus I don't know how long it is beforehand. Is there a

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

In R, is there any way that I can let R not printing out the "Read 4 items " message? > a<-scan("probes.txt") Read 4 items Thanks Hua LI

Hi All, I'm working with R and want to ignore the warning messages given, is there a way to stop R from giving out warning messages any more? an example: tt = "test" as.numeric(tt) would give me the following message: [1] NA Warning message: NAs introduced by coercion I decide to ignore the warning message for now and don't want it to show any more, can someone help?

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

You need to extend docs/schemas/domaincommon.rng which is what is used for the validation virsh is reporting. Do a 'git log' on that file for other examples of patches adding new qemu command line options. If you plan to submit it to libvir-list you will also want to add tests/ - Cole On 9/11/19 2:16 AM, Su Hua wrote: > Hi, I added libvirt4.5.0 support for vhost-user-blk, which can

Dear all, I'm using the "ape" package in R and want to draw a phylogenetic tree with not only the tip labels but also some labels for the edges. e.g. Mark the edge AB as "m" in the tree ABC. Couldn't find a way to do that. Can someone help? Thanks, Hua

Dear Johnny, Yes, I have tried much larger numbers than 60 and 3 for the above two parameters respectively. And I am sure it is using ?protocol 2? because it?s uncommented in sshd_config. Is there a way to catch what?s happing before quit? I couldn?t see anything except for the line ?write failed, broken pipe?. Thanks. Hua At 2015-10-08 19:37:24, "Johnny Hughes" <johnny at