search for: pp1

I'm trying to obtain numerical derivative of a probability computed with mvtnorm with respect to its parameters using grad() and jacobian() from NumDeriv. To simplify the matter, here is an example: PP1 <- function(p){ thetac <- p thetae <- 0.323340333 thetab <- -0.280970036 thetao <- 0.770768082 ssigma <- diag(4) ssigma[1,2] <- 0.229502120 ssigma[1,3] <- 0.677949335 ssigma[1,4] <- 0.552907745 ssigma[2,3] <- 0.784263100 ssigma[2,4]...

...# in the file list, I have all the filenames. n=length(list[[1]]) for (i in 1:n){ aux <- paste("p",i,sep="") assign(aux, as.matrix(read.table(list[[1]][i]))) } R creates the matrices p1,p2,...,pn. I want to manipulate the matrices p1,...,pn and create new matrices pp1,...,ppn. The dimensions of matrices ppi must be 100x100. So, I do the following: for (i in 1:n){ aux2 <- paste("pp",i,sep="") assign(aux2, array(0,dim=c(100,100))) } My question is the following: I want to assign to the matrix ppi[1,1] the value of pi[1,1] for i=...

Dear all, I'm looking for primary help at aggregating table() results and at writing a loop (if useful) My dataset ( http://goo.gl/gEPKW ) is composed of 23k rows, each one representing a point in the space of which we know the land cover over 10 years (column y01 to y10). I need to analyse it with a temporal sliding window of 5 years (y01 to y05, y02 to y06 and so forth) For each period

Can anyone explain what I'm missing here? max(pp1 <- package_version(c("0.99999911.3","1.0.4","1.0.5"))) ## [1] ?1.0.4? max(pp2 <- package_version(c("1.0.3","1.0.4","1.0.5"))) ## [1] ?1.0.5? I've looked at ?package_version , to no avail. Since max() goes to .Primitive("m...

Hi I looked, but I didn't find it: I need an array [0 to 101] where each element is a list (the result of Kest in spatstat). I.e. I want to do: A[2] <- Kest(pp1) A[3] <- Kest(pp2) ... A[101] <- Kest(pp100) How can I create A ? Rainer -- Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT) Department of Conservation Ecology University of Stellenbosch Matieland 7602 South Africa

...y1 mat1 <- cbind(c(index(ts01)[1], index(ts01), index(ts01)[length(index(ts01))]), c(ymin, coredata(ts01), ymin)) ? #vertices for area under y2 mat2 <- cbind(c(index(ts01)[1],index(ts01), index(ts01)[length(index(ts01))]), c(ymin, y2, ymin)) ? pp2 <- as(mat2, "gpc.poly") pp1 <- as(mat1, "gpc.poly") ? m <-max( (max(coredata(ts01)) - bl)/3, abs(min(coredata(ts01)) - bl)/3) m1 <- bl+m m2 <- bl+2*m ? ##### cutting the blue poly above the baseline # poly above the baseline s1 <- setdiff(pp1, pp2) y2 <- rep(m1, length(index(ts01))) mat2a...

Hi I have two processes which take with a certain probability (p1 and p2) x number of years to complete (age1 and age2). As soon as thge first process is completed, the second one begins. I want to calculate the time it takes for the both processes to be completed. I have the following script which gives me the answer, butI think there must be a more elegant way of doing the calculations

Hi, Thanks for the reply. Thats where the problem starts. If I set ppp0 as the default gw the internet doesnt work anymore. This is how im doing it... route del default route add default gw 217.32.81.74 dev ppp0 if I put it back to ppp1... route del default route add default gw 217.32.68.73 dev ppp1 It works fine again. Whats up with that? Cheers, Chris ----- Original Message ----- From: "andybr" <andybr@bol.com.br> To: <cbolton@hirstanddanson.co.uk> Cc: <lartc@mailman.ds9a.nl> Sent: Thur...

...PR#1964) Could you us a hint what the problems are? And an example to reproduce them? On Fri, 30 Aug 2002 stakb@nus.edu.sg wrote: > Full_Name: Kaushik Bhattacharyya > Version: 1.5.1 > OS: Solaris > Submission from: (NULL) (137.132.3.10) > > > Main R-function used: > > pp1 > > function(X) > > { > > x_as.matrix(X) > > trial_function(a) > > { > > clusproj(x,a) > > } > > test.nlm_nlm(trial,rep(1,ncol(X))) > > theta_test.nlm$estimate > > theta_theta/sqrt(sum(theta^2)) > > cluster.index_round(1/test.nlm$min...

...39;m running into a few problems. Its got two adsl modems connected to it, both connected to seperate 512kbs lines. Now I''ve followed the simple source routing in the advanced routing howto to the letter but it doesnt work. I''ve got it autoconnecting on startup and redhat puts ppp1 as the default gateway, this is then setup for masquerading for the entire network. Therefore I''ve tried setting up ppp0 as the deafult gateway for only one computer (10.0.0.11), as it says at http://lartc.org/howto/lartc.rpdb.html#LARTC.RPDB.SIMPLE I''ve done everything it says t...

(The previous sent message had some issues...) ---------------------------------------------- Hello to everybody, I just want make a quick summary: OBJECTIVE: discuss for finding a way to represent logical (derived from source language abstractions) address space also in the IR to be able to exploit this information for optimizations. CURRENT STATE: The partial workaround probably used IMO

...s share some features with other posters on this thread: * different sized address regs * most address regs are larger than int regs * address regs only support a sub-set of operations * some address spaces are disjoint, some may overlap. so far, I declared a data layout with "-p0:xx:yy-pp1:zz:aa-p2:... " etc. ... implemented getPointerWidthV() (and getPointerRegClass(), which never seems to get called). I see that global data items are correctly sized according to my pointer definitions - but that all operations are carried out as if pointers were sized for address space 0. W...

Dear R users, I have a function (below) which encompasses several tests. However, when I run it, only the output of the last test is displayed. How can I ensure that the function root(var) will run and display the output from all tests, and not just the last one? Thank you, b. root <- function(var) { #---Phillips-Perron PP.test(var, lshort = TRUE) PP.test(var, lshort = FALSE)

Hello to everybody, I would like to start a discussion about a possible extension of address space concept in LLVM. The idea was born starting from this discussion in the clang mailing list (first msg: http://lists.cs.uiuc.edu/pipermail/cfe-commits/Week-of-Mon-20130715/084011.html - interesting point: http://lists.cs.uiuc.edu/pipermail/cfe-commits/Week-of-Mon-20130722/084499.html) where

...: API unit is 2^-32 Joules, 5 fixed counters, 655360 ms ovfl timer [ 0.667411] RAPL PMU: hw unit of domain pp0-core 2^-14 Joules [ 0.667412] RAPL PMU: hw unit of domain package 2^-14 Joules [ 0.667412] RAPL PMU: hw unit of domain dram 2^-14 Joules [ 0.667413] RAPL PMU: hw unit of domain pp1-gpu 2^-14 Joules [ 0.667413] RAPL PMU: hw unit of domain psys 2^-14 Joules [ 0.667772] clocksource: tsc: mask: 0xffffffffffffffff max_cycles: 0x2879c5f06f2, max_idle_ns: 440795220049 ns [ 0.667991] simple-framebuffer simple-framebuffer.0: framebuffer at 0xb0000000, 0x1d5000 bytes, mapped t...

...27000] (64MB) mapped at [ffff880081c27000-ffff880085c26fff] [ 1.741665] RAPL PMU detected, API unit is 2^-32 Joules, 3 fixed counters 163840 ms ovfl timer [ 1.741674] hw unit of domain pp0-core 2^-16 Joules [ 1.741679] hw unit of domain package 2^-16 Joules [ 1.741684] hw unit of domain pp1-gpu 2^-16 Joules [ 1.745047] futex hash table entries: 2048 (order: 6, 262144 bytes) [ 1.745190] audit: initializing netlink subsys (disabled) [ 1.745258] audit: type=2000 audit(1449332805.708:1): initialized [ 1.746590] Initialise system trusted keyring [ 1.746841] HugeTLB registere...

...: API unit is 2^-32 Joules, 4 fixed counters, 655360 ms ovfl timer [ 8.811025] RAPL PMU: hw unit of domain pp0-core 2^-14 Joules [ 8.811026] RAPL PMU: hw unit of domain package 2^-14 Joules [ 8.811027] RAPL PMU: hw unit of domain dram 2^-14 Joules [ 8.811028] RAPL PMU: hw unit of domain pp1-gpu 2^-14 Joules [ 8.850141] usb 1-1.1: New USB device found, idVendor=0451, idProduct=8142 [ 8.850144] usb 1-1.1: New USB device strings: Mfr=0, Product=0, SerialNumber=1 [ 8.850146] usb 1-1.1: SerialNumber: 1F0008593B9B [ 8.851590] hub 1-1.1:1.0: USB hub found [ 8.851654] hub 1-1.1...

Hello, I am encountering a 30 second hang during boot, in the Xorg process just before the display manager comes up. I have bisected the problem to the following commit: commit 2b700825e7a7702fb862edba1262c98040dc1bf6 Author: Ben Skeggs <bskeggs at redhat.com> Date: Thu Aug 20 14:54:22 2015 +1000 drm/nouveau/mc: move device irq handling to platform-specific code The hang only

...89000] (64MB) mapped at [ffff8800b9a89000-ffff8800bda88fff] [ 0.482861] RAPL PMU detected, API unit is 2^-32 Joules, 3 fixed counters 163840 ms ovfl timer [ 0.482864] hw unit of domain pp0-core 2^-16 Joules [ 0.482866] hw unit of domain package 2^-16 Joules [ 0.482867] hw unit of domain pp1-gpu 2^-16 Joules [ 0.483074] microcode: CPU0 sig=0x306a9, pf=0x10, revision=0x1b [ 0.483087] microcode: CPU1 sig=0x306a9, pf=0x10, revision=0x1b [ 0.483103] microcode: CPU2 sig=0x306a9, pf=0x10, revision=0x1b [ 0.483115] microcode: CPU3 sig=0x306a9, pf=0x10, revision=0x1b [ 0.483208]...

Hello, following up on myself, it was suggested on IRC that I better attach a dmesg output. Here's the output of a clean boot & echo 0f > /sys/.../pstate cycle. I can't spot a message that relates to the reclock action, and there's only one weird "nouveau 0000:04:00.0: clk: base: 7 MHz, boost: 7 MHz" message. On the other hand: # cat