Displaying 20 results from an estimated 1327 matches for "m2".
Did you mean:
m2p
2006 Mar 13
3
hfsc and dropped packets
Hi,
I''m trying to get a handle on hfsc. Here is my configuration:
root@jmnrouter:/jmn# tc class show dev vlan1
class hfsc 1: root
class hfsc 1:1 parent 1: ls m1 0bit d 0us m2 225000bit ul m1 0bit d 0us m2
225000bit
class hfsc 1:10 parent 1:1 rt m1 191000bit d 25.0ms m2 135000bit ls m1 0bit
d 0us m2 135000bit ul m1 0bit d 0us m2 225000bit
class hfsc 1:20 parent 1:1 rt m1 22008bit d 25.0ms m2 56008bit ls m1 0bit d
0us m2 56008bit ul m1 0bit d 0us m2 225000bit
class hfsc 1...
2006 Apr 23
2
distribution of the product of two correlated normal
Hi,
Does anyone know what the distribution for the product of two correlated
normal? Say I have X~N(a, \sigma1^2) and Y~N(b, \sigma2^2), and the
\rou(X,Y) is not equal to 0, I want to know the pdf or cdf of XY. Thanks
a lot in advance.
yu
[[alternative HTML version deleted]]
2010 Jun 16
2
data frame
...CLUSTER YEAR variable Delta R_pivot
M1 2005 EC01 NA NA
M1 2006 EC01 2 NA
M1 2007 EC01 4 5
M2 2005 EC01 NA NA
M2 2006 EC01 5 NA
M2 2007 EC01 8 7
M1 2005 EC02 NA NA
M1...
2004 Jul 13
1
MLE, precision
Hi, everyone
I am trying to estimate 3 parameters for my survival
function. It's very complicated. The negative
loglikelihood function is:
l<- function(m1,m2,b) -sum( d*( log(m1) + log(m2)
+ log(1- exp(-(b + m2)*t)) ) + (m1/b - d)*log(m2 +
b*exp(-(b + m2)*t) ) + m1*t - m1/b*log(b+m2) )
here d and t are given, "sum" means sum over these
two vairables.
the parameters are assumed small, m1, m2 in
thousandth, m2 in millionth.
I used t...
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
Dear R users,
I'm trying to fit a gamma-frailty model on a simulated dataset, with 6 covariates, and I'm running into some results I do not understand. I constructed an example from my simulation code, where I fit a coxph model without frailty (M1) and with frailty (M2) on a number of data samples with a varying degree of heterogeneity (I'm running R 2.3.1, running takes ~1 min).
library(survival); set.seed(10000)
lambda <- 0.01 # Exp. hazard rate
# Beta coefficients for Age,TC,HDLC,SBP,Diab,Smok
beta <- c(0.0483,0.0064,-0.0270,0.0037,0.4284,0.5234)...
2010 Jun 17
1
big big problem
...TER YEAR variable value1 value2
M1 2005 EC01 NA NA
M1 2006 EC01 2 5
M1 2007 EC01 4 5
M2 2005 EC01 NA NA
M2 2006 EC01 5 7
M2 2007 EC01 8 7
M1 2005 EC02 NA NA...
2009 Aug 18
1
create a table in the console!!
...quot;, sep=",")
> Data2 <- read.csv2("c:/GWPMax1.csv",sep=",")[1:20,1:2]
> M <- merge(Data1,Data2, by.x = "Policy.Number",by.y = "Policy.Number",all.x = TRUE,all.y = TRUE )
> (No_GWPMax<-nrow(M[M[,25]=="NA",]))
[1] 8
> M2<- merge(Data1,Data2, by.x = "Policy.Number",by.y = "Policy.Number")
> M2$GWP_Max <- as.numeric(as.character(M2$GWP_Max ))
> class1 <- M2[M2[,25]>0 & M2[,25]<1000,]
> (NbpolicyClass1 <- nrow(class1))
[1] 5
> class2 <- M2[M2[,25]>1000 &...
2009 Aug 18
1
Tr : create a table in the console!!
...quot;, sep=",")
> Data2 <- read.csv2("c:/GWPMax1.csv",sep=",")[1:20,1:2]
> M <- merge(Data1,Data2, by.x = "Policy.Number",by.y = "Policy.Number",all.x = TRUE,all.y = TRUE )
> (No_GWPMax<-nrow(M[M[,25]=="NA",]))
[1] 8
> M2<- merge(Data1,Data2, by.x = "Policy.Number",by.y = "Policy.Number")
> M2$GWP_Max <- as.numeric(as.character(M2$GWP_Max ))
> class1 <- M2[M2[,25]>0 & M2[,25]<1000,]
> (NbpolicyClass1 <- nrow(class1))
[1] 5
> class2 <- M2[M2[,25]>1000 &...
2009 Jul 20
2
I need to obtain all the rows in m1 in which m2 the elements of m2 are present
...quot;,"6","7"),
ncol=2)
m1
[,1] [,2]
[1,] "a" "1"
[2,] "7" "2"
[3,] "a" "3"
[4,] "i" "4"
[5,] "o" "5"
[6,] "u" "6"
[7,] "i" "7"
m2<-c("a","7")
I need to obtain all the rows in m1 in which m2 the elements of m2 are
present
[,1] [,2]
[1,] "a" "1"
[2,] "7" "2"
[3,] "a" "3"
[7,] "i" "7"
Thx again
--
Alberto Lora Mich...
2006 May 09
3
tc del class not working
When I start my script:
* - Creating classes on br1 for upload control ...
* - tc class add dev br1 parent 2:0 classid 2:46 hfsc ls m1
576.0Kbit d 2000ms m2 192.0Kbit ul m2 384Kbit ...
[ ok ]
* - tc class add dev br1 parent 2:46 classid 2:47 hfsc sc umax
1500b dmax 30ms rate 80Kbit . [ ok ]
* - tc class add dev br1 parent 2:46 classid 2:48 hfsc ls m2
152.0Kbit ul m2 152.0Kbit ... [ ok ]
* - tc class add dev br1 parent 2:46 classid 2:...
2005 Oct 13
12
HFSC Advanced Limiting
I am investigating changing our broadband bandwidth control scripts over to HFSC. What I would like to accomplish is this:
{ * DOWNLOAD *
I need a download class for each user (we''ll call it Class 1) that offers a minimum of lets say 384Kbps and a maximum of 1.5Mbps.
The rates between 384Kbps and 1.5Mbps should be shared among all users.
Each user should always be able to get at
2011 Feb 23
0
parallel bootstrap linear model on multicore mac
People of R(th),
I have been ramming my head against this problem, and I wondered if
anyone could lend a hand. I want to parallelize a bootstrap of a linear
model on my 8-core mac. Below is the process that I want to parallelize
(namely, the m2.ph.rlm.boot<-boot(m2.ph,m2.ph.fun, R = nboot) command).
This is an extension of the bootstrapping linear models example in
Venables and Ripley to lmer. Please excuse my rather terrible
programming skills. I am always open to suggestions. Below the example I
describe what methods I have tried...
2004 May 27
2
block diagonal matrix function
...matrices and rowname() and colname() inherited appropriately.
Two questions:
(1) Is there a better way to do this? (kronecker() isn't applicable here)
(2) How do I generalize it to take an arbitrary number of matrices as
inputs?
TIV
Robin
"blockdiag" <-
function (m1, m2, p.tr = 0, p.ll = 0)
{
## p.tr and p.ll are padding values
topleft <- m1
topright <- matrix(p.tr, nrow(m1), ncol(m2))
colnames(topright) <- colnames(m2)
lowleft <- matrix(p.ll, nrow(m2), ncol(m1))
lowright <- m2
rbind(cbind(topleft, topright), cbind...
2008 Jan 29
2
add/subtract matrices, ignoring NA or missing values
Hi,
For example, given two 2x2 matrices m1 and m2. I would like to add/subtract
element by element
> m1
[,1] [,2]
[1,] NA NA
[2,] 1 2
> m2
[,1] [,2]
[1,] 1 NA
[2,] NA 2
> m1 + m2
[,1] [,2]
[1,] NA NA
[2,] NA 4
How can I ignore the NA, and get this ? Hope the solution can be extende...
2013 Jun 27
2
[LLVMdev] Problem with linking modules which use a shared type
Hi,
I stumbled upon a strange thing regarding types and module linking, which I don't quite get, so maybe someone around here can enlighten me.
Consider the following code, which will create 2 modules m1 and m2, and a named structured type %T = { i32 }; m1 contains only a function definition f(%T), m2 contains a function declaration f(%T) and a function definition h(%T), where h will call f in its body. Note that all functions in m1 and m2 are based upon the same structured type %T:
LLVMContext cont...
2006 May 11
9
HFSC and prioritization
...the user is not doing anything else. No matter
how much bandwidth P2P wants, if something else needs bandwidth, the
other traffic should win out and receive the bandwidth.
These are my actual rules:
# Base user class
tc class add dev wivl4 parent 5:0 classid 5:130 hfsc ls m1 1536.0Kbit d
2000ms m2 512.00Kbit ul m2 1024Kbit
# Priority queue
tc qdisc add dev wivl4 parent 5:130 handle 134: prio bands 3
tc qdisc add dev wivl4 parent 134:1 handle 135: hfsc default 1
tc qdisc add dev wivl4 parent 134:2 handle 136: hfsc default 1
tc qdisc add dev wivl4 parent 134:3 handle 137: hfsc default 1...
2010 Apr 14
1
envelope in spatstat
...= wind)
However the marks are alignments (lines), and so have to be treated
differently to normal scalar marks. Hence to create a mcf object with the
appropriate test function for alignment marks, I input 'func' (below)
suggested by Stoyan & Penttinen (1989):
func <- function(m1,m2) { sin(abs(m1-m2))^2}
mcf <- markcorr(points, func, normalise = TRUE, method="density")
So far, so good. However, usinf 'envelope' and 'rlabel' I would like to
check if the pattern in the data is lost when randomly relabeling the mark
for each point. If the t...
2010 May 14
1
Replacing double loop by apply
...he following question.
is it possible to replace to following double loop by some apply calls?
#######################################################################
m1 <- data.frame(v1=factor(letters[1:5]),
v2=factor(letters[2:6]),
v3=factor(letters[3:7]))
m2 <- data.frame(v1=factor(letters[3:7]),
v2=factor(letters[1:5]),
v3=factor(letters[4:8]))
ind <- matrix(logical(nrow(m2) * ncol(m2)),
nrow = nrow(m2), byrow = TRUE)
for (j in 1:ncol(m2))
{
for (i in 1:nrow(m2))
{
ind[i, j] <- any...
2007 Feb 27
2
Creating a contended section of bandwidth with HTB and IMQ
Hi All,
I''m trying to create a contended section of bandwidth using IMQ. I have the
imq0 device up and running, with traffic passing through it.
Firstly, I need to throttle the entire device imq0 to 2mbit/s.
I would then like to add throttle rules for individual IP addresses, allowing
them to pass up to 512kbit/s each, as long as imq0 has not reached its
2mbit/s.
The configuration
2004 Dec 06
2
Re : LOOPS
...odes.
Basically, the following codes were written to establish the sire gametes
or phases for daughter design for one markers two alleles.
Here are the SAS code:
do i=1 to 744;
do j=745 to 1540;
m[j]=0;
if sire[j]=anml[i] then do;
if m1[j]=m1[i] and m2[j]=m2[i] then m[j]=0;
else if m1[j]=m1[i] then m[j] =1;
else if m1[j]=m2[i] then m[j==2;
else if m2[j]=m1[i] then m[j]=1;
else if m2[j]=m2[i] then m[j]=2;
else m[j]=0;
end do;
end;
end;
Thanks Stella...