search for: p3

...-theta[6])) ## define the diffusion function S <- function(t,x,theta) theta[5]*(x^theta[6]) ## define the transformed drift function and it's first six derivatives m0 <- function(t,x,theta) { b <- theta[5]*(theta[6]-1) p1 <- (theta[6]+1)/(theta[6]-1) p2 <- theta[6]/(theta[6]-1) p3 <- (2-theta[6])/(1-theta[6]) (-1/theta[5])*(theta[1]*((b*x)^p1) - theta[2]*((b*x)^p2) + theta[3]*b*x - theta[4]*((b*x)^p3) - 0.5*theta[6]*(theta[5]^2)*((b*x)^(-1))) } m1 <- function(t,x,theta) { b <- theta[5]*(theta[6]-1) p1 <- (theta[6]+1)/(theta[6]-1) p2 <- theta[6]/(theta[6]-1) p3...

I have 3 vectors: p1, p2, and p3. I would like each vector to be any possible value between 0 and 1 and p1 + p2 + p3 = 1. I want to graph these and I've thought about using scatterplot3d(). Here's what I have so far. library(scatterplot3d) p1 <- c(1,0,0,.5,.5,0,.5,.25,.25,.34,.33,.33,.8,.1,.1,.9,.05,.05) p2 <- c(0,1...

Suppose I have three numbers p1, p2, p3 with 0 <= p1,p2,p3 <= 1 and p1+p2+p3=1, and a function f=f(p1,p2,p3) = f(p1,p2,1-p1-p2). How to draw a contour plot of f() on the p1+p2+p3=1 plane, that is, an equilateral triangle? Functions triplot(), triangle.plot(), and ternaryplot() give only scatterplots, AFAICS -- Robin...

Hi, I am wondering if anyone can propose a simple/best way to do the following: Let's say I have a data frame dat <- cbind(expand.grid(mode=c('right','left'),time=0:3,id=c('p1','p2','p3')),y=c(3,5,rep(4,6),6,2,rep(3,6),4,4,rep(2,6))) dat mode time id y 1 right 0 p1 3 2 left 0 p1 5 3 right 1 p1 4 4 left 1 p1 4 5 right 2 p1 4 6 left 2 p1 4 7 right 0 p2 6 8 left 0 p2 2 9 right 1 p2 3 10 left 1 p2 3 11 right 2 p2 3 12 left 2...

Dear All, Is there a way I can turn off the following warning message for using multi-argument returns? multi-argument returns are deprecated in: return(p1, p2, p3, p4) Steve. ************************************************************************************** Steve Su (s.su@qut.edu.au) Postdoctoral fellow Faculty of Business Queensland University of Technology Postal Address: Steve Su, School of Accountancy, QUT, PO Box 2434, Brisbane,...

Dear all, I have two functions (f1, f2) and 4 unknown parameters (p1, p2, p3, p4). Both f1 and f2 are functions of p1, p2, and p3, denoted by f1(p1, p2, p3) and f2(p1,p2,p3) respectively. The goal is to maximize f1(p1, p2, p3) subject to two constraints: (1) c = k1*p4/(k1*p4+(1-k1)*f1(p1,p2,p3)), where c and k1 are some known constants (2) p4 = f2(p1, p2, p3) In add...

Dear R-experts, I have got a dataframe: data ID disease V1 V2 1 p1 1 2 p1 3 3 p3 3 4 p3 5 5 p5 1 From which I extract a usefull table: affect affect 1 3 5 p1 1 1 0 p3 0 1 1 p5 1 0 0 I want to merge this with anotherdataframe: age p1 23 p2 24 p3 23 p4 11 p5 45 If have tried as.data.frame(affect) and other solutions to get the following comment going:...

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120213/d2a5afa6/attachment.pl>

...done. My data looks something like: Participant Stimulus Measurement p1 s`1 5 p1 s`2 6.1 p1 s`3 7 p2 s`1 4.8 p2 s`2 6 p2 s`3 6.5 p3 s`1 4 p3 s`2 7 p3 s`3 6 As a first step I would imagine that I have to rearrange my data into a frame more like this Stimulus p1 p2 p3 s1 5 4.8 4 s2 6.1 6 7 s3 7...

Among others, datam contains the columns: logconc, tm, dose, subj, bilirubin. None of these are factor variables. The following compartment models work (the first still has not converged after 100 interations): res1 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))* (exp(-exp(p2)*tm)-exp(-exp(p1)*tm))),start=list(fixed=c(5,-2,-0.1)), fixed=list(p1+p2+p3~1),control=list(maxIter=100), groups=~subj,data=datam,verbose=T,method="ML") res2 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))* (exp(-exp(p2)*tm)-exp(-exp(p1)*t...

Hi, I installed centos 4 on a P4, configured it (and it booted fine here) then sent it to the DC to be put in a P3 but it won't boot, displays GRUB and halts. The DC guys tried grub-install in rescue mode but still the same. Any clues? Is there some difference between P3 and P4? thanks tom

...gfix Advisory 2011:1836 Upstream details at : https://rhn.redhat.com/errata/RHBA-2011-1836.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 9ea016735ff112bfae8ef2384d9cb0749ba34575c100e71a51befcbf1c9cfac2 bind-9.7.3-8.P3.el6_2.1.i686.rpm 272f07955c58464ca5d16ca62d6ec9a19a4ca0f43652af8a45954a8ccc243b1a bind-chroot-9.7.3-8.P3.el6_2.1.i686.rpm 2d6933d09473791b540b023ef798ff0607dbf09db7cac740f17f93470c5c5e9d bind-devel-9.7.3-8.P3.el6_2.1.i686.rpm ed825c6686057e755e6189bee596edb4f1ec7712d4dfbb5da655ea85a8ebda7f bind-...

...gfix Advisory 2012:0009 Upstream details at : https://rhn.redhat.com/errata/RHBA-2012-0009.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 5e336743367a5824ee075ee94cd55f7e84e106131cff72eb8549bd52bab580d8 bind-9.7.3-8.P3.el6_2.2.i686.rpm 0a2e1e34f8fff5e0eee35f1959d599e9d5523ceade4ec147552d63417e0b2b8b bind-chroot-9.7.3-8.P3.el6_2.2.i686.rpm 1551d76b2d3a564df94e2e368e7ce9c1cf7e47245f9cba70bdf0ebd7288fbe9f bind-devel-9.7.3-8.P3.el6_2.2.i686.rpm b91bf00604644abd34e832ac2865918f502bceccb451b413b69fedd4af5bdb7d bind-...

...sory 2012:0716 Important Upstream details at : https://rhn.redhat.com/errata/RHSA-2012-0716.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 081b3e564fa77e8042355bf02848a1384b8321391e19d11d1348a8c3b220a72d bind-9.7.3-8.P3.el6_2.3.i686.rpm bac48abad13d32c664a68e8730558e65a2ed280d732fc6db37109efa502a7a0f bind-chroot-9.7.3-8.P3.el6_2.3.i686.rpm 2d2debcbb81feac2ca95aed1dd8d0d138ac1ec65cc752b4a3c020f17e6601938 bind-devel-9.7.3-8.P3.el6_2.3.i686.rpm b450fbefc48d5c6fe9e624fb54f5df46b03e02a9eb69da2cc019e6d00e160600 bind-...

The following dummy data frame has factor Q (with 2 levels) nesting factor P (with levels p1 and p2 nested under q1, and p3 and p4 nested under q2), but both crossing the random variate s, which has 8 levels. The dependent measure is dv. > # The data frame: > testnest dv s P Q 1 1 s1 p1 q1 2 2 s2 p1 q1 3 1 s3 p1 q1 4 2 s4 p1 q1 5 1 s5 p1 q1 6 3 s6 p1 q1 7 3 s7 p1 q1 8 4 s8 p1 q1 9 2...

...gfix Advisory 2012:0009 Upstream details at : https://rhn.redhat.com/errata/RHBA-2012-0009.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 5e336743367a5824ee075ee94cd55f7e84e106131cff72eb8549bd52bab580d8 bind-9.7.3-8.P3.el6_2.2.i686.rpm 0a2e1e34f8fff5e0eee35f1959d599e9d5523ceade4ec147552d63417e0b2b8b bind-chroot-9.7.3-8.P3.el6_2.2.i686.rpm 1551d76b2d3a564df94e2e368e7ce9c1cf7e47245f9cba70bdf0ebd7288fbe9f bind-devel-9.7.3-8.P3.el6_2.2.i686.rpm b91bf00604644abd34e832ac2865918f502bceccb451b413b69fedd4af5bdb7d bind-...

Dear R user I follow the steps defined in Modern applied statistics page(453) to use optim. However, when I run the following code the parameters seems way off and the third parameter(p3) stayed as the initial value. below is the code: ## data da=c(418,401,416,360,411,425,537,379,484,388,486,380,394,363,405,383,392,363,398,526) ### initial values pars=c(392.25, 507.25, 0.80) p1=pars[1];p2=pars[2];p3=pars[3] ### Objective function mix.obj = function(p,x) { e =p[3]*(pnorm((x+1/...

Hi everyone, I need to recode multiple columns in a dataframe into a single column in a variety of different ways. Often the values will be TRUE/FALSE and I want a list of the columns that are true as in the Result column below: P1 P2 P3 P4 Result 1 0 0 1 1 P3,P4 2 0 1 1 1 P2,P3,P4 3 1 0 0 0 P1 4 0 0 0 0 NA 5 1 1 1 1 P1,P2,P3,P4 I'm still relatively new to R and tend to think in for loops. I'm sure there exists a more concise and elegant solution....

...(ÈËÃñ±Ò) 1.cpu InteI p4 2G/1.8G/1.7G(Socket478) £¤ 500/420/390Ôª InteI p4 1.6G/1.5G/1.4G(Socket478) £¤ 380/350/320Ôª InteI ÈüÑï4 1.7G ºÐ/É¢(Socket478) £¤ 250/220Ôª InteI ÈüÑï1G/950/900/850 £¤ 210/200/190/180Ôª InteI ÈüÑï8OO/733/ £¤ 170/150Ôª InteI ÈüÑïTuaIatin 1.3G/1.2G/1G £¤ 250/210/200Ôª InteI P3 933EB/866EB £¤ 500/390Ôª AMD AthIon XP 1800+/1700+/1600+ £¤ 350/320/300Ôª AMD À×Äñ 1.33G/900/850 £¤ 300/280/250Ôª AMD Duron 1.2G/1G/900/850/800 £¤ 250/210/180/120/110Ôª 2.Ö÷°å InteI i845/i845D ASUS P4B/ABIT BL7/B266-C/ABIT BD7 £¤ 500/410/520/390Ôª ¼¼¼Î GA 81DX/81DX3/81RXP/81PX £¤ 340/410/510/420Ôª...

...data frame. Plus, I cannot query names(mydata) to find the column index which header should contain "~" or "." > mydata <- data.frame(read.table("mycsv.csv", sep = ",", header = TRUE, na.strings = ".", as.is = TRUE)) > mydata P1 P2 P3 P1.P2 P1.P3 P2.P3 A B C 1 1 2 3 4 5 6 7 8 9 > names(mydata) [1] "P1" "P2" "P3" "P1.P2" "P1.P3" "P2.P3" "A" "B" "C" > grep("~",names(mydata)) integer(0) &gt...