search for: zhenning

>From: Sundar Dorai-Raj <sundar.dorai-raj at PDF.COM> >Reply-To: sundar.dorai-raj at PDF.COM >To: Zhen Pang <nusbj at hotmail.com> >CC: r-help at stat.math.ethz.ch >Subject: Vectorising and loop (was Re: [R] optim "a log-likelihood >function") >Date: Wed, 29 Sep 2004 18:21:17 -0700 > > > >Zhen Pang wrote: > >> >>I also use

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz<-sapply(ma, myfunction) To my surprise, It almost costs me the same time to get my results if I directly use a loop

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz<-sapply(ma, myfunction) To my surprise, It almost costs me the same time to get my results if I directly use a loop

Dear all, Using R function "dwt", it seems that I cannot specify the wavelet transformation like Haar. What's the R code for wavelet decomposition which allows me to specify Haar wavelet transformation? Of course, if it can include "db2", that is even better. In general, I want an R function like matlab code "dwt". Thanks in advance! Zhen Li

Yes. It should have something to do with read/write permissions, but it is not clear how it happens. I can write file to C drive using R. I usually write my results matrix to a txt file in C drive. For Rprof(), the boot.out file can be created, but only with one line sample.interval=20000 The situation is the same even if I specify the directory to the D drive,where I have the full

Dear all, I have made my testing program to run successfully under unix in the background. However, my simulation work does not work. I read the foo.results file, I found it only have part of my code and not any output I want. Is there any line limit? My code is nearly 400 line. I can cut some of them, but I want to know whether there is any limit or exactly the number of limit is. Thanks.

Dear all, Our department uses the linux system and we are not allowed to submit job directly. We must make a batch to submit through "qmon". so, I make a foo.sh file, which only contains one line: nohup R --vanilla < foo.txt > foo.results foo is all my codes. It is a simulation of 200 times. I set the seed at the beginning. It is to estimate the success probability, which is

Dear all, I would like to sort the values of "member" for each group, and obtain a variable to indicate its rank within the group. For example, we have original dataset as follows: df <- data.frame(group = c(rep(1, 3), rep(5, 3)), member = c(30, 10, 22, 21, 44, 15)) group member 1 1 30 2 1 10 3 1 22 4 5 21 5 5 44 6 5 15 After

Dear LLVM developers, Recently I want to try the pass '-loop-unroll-and-jam' to observe how the IR is transformed, but I don't see the IR changed after doing the following steps. Here are the details for preparing my experiment and I have tried the LLVM 6, 7 and 8: 1) Simple 2D array source code (loop.c) is given ``` #define M 32768 #define N 32768 double a[M][N], b[M][N],

https://bugzilla.netfilter.org/show_bug.cgi?id=1685 Bug ID: 1685 Summary: Calling the nftnl_set_free function may trigger the "double free" problem. Product: libnftnl Version: unspecified Hardware: All OS: All Status: NEW Severity: critical Priority: P5

Full_Name: Zhen Liu Version: 1.5.1 OS: Window XP Home/Window NT 4.0 Submission from: (NULL) (205.142.0.75) I have installed R-1.5.1 using SetupR.exe on both my home laptop (Sony Vaio running Window XP Home) and in office (Compaq AP550 running Window NT 4.0), but they shared the same problem. Both worked for a while, but then Rgui.exe crashed on the startup, complains: The exception

I have problem evaluating the expression h = exp(x)/(exp(exp(x))-1) for large negative x. This expression is actually the probability that y = 1 when y is a Poisson random variable truncated at 0, hence must satisfy 0 <= h <= 1. However, when x < -18, I may get an h value that is larger than 1 while the true value should be a quantity that is smaller than but very close to 1. For

I''ve implemented code to send out bogus gratuitous ARP packet from vna in order to fix CR 6701114. The webrev is at: http://jurassic.eng/net/consulte.prc/export/build/xvm-6701114/webrev Some backgroud information (see CR 6701114 for more detail): During live migration of an xVM domain from one dom0 to another, the VNIC will be moved from one switch port to another. But, the

Dear professor, I'm using package nlme, but I can't find a way to specify two uncorrelated random effects. For example, a random intercept and a random slope. In package lme4, we can specify&nbsp;x + (x ll g) to realize, but how in nlme? Thanks! ???????????????????????? Zhen Wang Graduate student, Department of Medical Statistics, School of Public Health, Sun Yat-sen

# HG changeset patch # User max.zhen@sun.com # Node ID db465ecb164e6d89290ae3d4fcc029e29305ae7d # Parent c682ae7a17224dc5576ba3c5db47766b2a7dc159 Shouldn''t call domain_create_ring() for a domU which is shutting down. There is a race condition between domain_create_ring() and the removal of the do mU entry from xenstore. If the removal of the domU entry from xenstore comes lat er than the

# HG changeset patch # User max.zhen@sun.com # Date 1171024588 28800 # Node ID a2d15838ea72ee049d8a657e9ecd406b5c6cb53e # Parent d609929f140809ec1c9a059f2eeaffb7ce56c891 Xenconsoled should ignore spurious watch event. Otherwise, xenconsoled could bind to a dying domU during its suspending and cause complaint from Xen: (XEN) event_channel.c:177:d0 EVTCHNOP failure: domain 0, error -22, line 177

Dear all, I wonder if there are some function or package in R which can do the iterative proportional fitting. In the exponential model f(y1,...yn)=exp(a'yi+b'(yi*yj)+.....+c'(y1*...*yn)+constance), instead of the canonical parameters, I use maginal probability instead of a and log odds ratio instead of b. and for the order higher than 3, I use the canonical way or even assume

Hello everybody, When I modify the dovecot imap source ,in "src/imap/" dir, I meet the Error below: Error: child (imap) kill (core dump disabled). In fact, I want to record the imap operation through dovecot. So, I modify every source file in "src/imap/cmd-IMAP_OPERATION" ,when the imap command function is executed, it will write to a file cur time and operation. Could you

In current Crossbow''s implementation, vlan tag seems to be unconditionally retained for promisc mode listener even when MAC_OPEN_FLAGS_TAG_DISABLE is not specified. I saw comments in mac_rx_deliver() saying that this is deliberately designed like this. I''m wondering why we design it like this (choose not to respect MAC_OPEN_FLAGS_TAG_DISABLE flag for promisc mode listener)? Or

Hi, I am doing a project involving checking a called specific function’s argument. Suppose that the function is int f(const char* str). When I am analyzing such a snippet: f("hello") , then compiled by Clang, I will have the “hello” as a *Constant Array* in the IR code. My goal is to call APIs of LLVM to get the “hello” from IR code. Now suppose the I got the llvm::Function* fn from