Displaying 20 results from an estimated 129 matches for "n2".

2012 Dec 04

2

Solve system of equations (nleqslv) only returns origin

...require(nleqslv)
###### Global Parameters ############
beeta=0.8
pq=10000
L=12600
theta=0.6
psale=0.6
mu=psale*(1-theta)
alphah=0.15
Cg=6240
Cs=2820
A= 100
D=0.0001
greekp=0.43
K=100000
##### Species Parameters ##########
b1=0.38
p1=16654
v1 = 0.28
N1=6000
g1=1
delta1=1
b2=0.4
p2=2797
v2 = 0.31
N2=10000
g2=1
delta2=1
### Define functions with vector x = c(Lg, Ls, gamma1, gamma2, lamda)
firstordercond <- function (x) {
y=numeric(4)
y[1]=(alphah/x[3])-(x[5]*((p1-(((theta+mu)*(((N1/A)*g1^greekp*x[1]^b1)+K))+((theta+mu)*(((1-exp(-2*D*v1*N1))*x[2])+K))))*(((N1/A)*g1^(greek...

2010 Jul 05

3

Function to compute the multinomial beta function?

Dear R-users,
Is there an R function to compute the multinomial beta function? That is, the normalizing constant that arises in a Dirichlet distribution. For example, with three parameters the beta function is Beta(n1,n2,n2) = Gamma(n1)*Gamma(n2)*Gamma(n3)/Gamma(n1+n2+n3)
Thanks in advance for any assisstance.
Regards,
Greg
[[alternative HTML version deleted]]

2013 Jan 03

0

help with NLOPTR

...aic)
###### Global Parameters ############
beeta=0.8
pq=10000
pf=10000
F=20
L=12600
theta=0.6
psale=0.6
mu=psale*(1-theta)
alphah=0.15
Cg=6240
Cs=2820
A= 100
D=0.0001
greekp=0.43
K=100000
##### Species Parameters##########
b1=0.38
p1=16654
v1 = 0.28
N1=6000
g1=1
delta1=1
b2=0.4
p2=2797
v2 = 0.31
N2=10000
g2=1
delta2=1
####################################### Objective function
############################################
eval_f = function (x) {
return(-1 * ( ( (1-alphah) *log(F) ) + ( alphah* (
(x[1]*(((N1/A)*(g1^greekp)*(x[3]^b1))+((1-exp(-2*D*v1*N1))*x[4])))
+ (x[2]*(((N2/...

2012 May 20

2

Remus network buffering problem

Hi all,
I have a following problem:
- Remus network buffering doesn''t work. It seems to be because of no vif is
reported by the function server.xend.domain on line 29 of
/usr/local/lib/python2.7/dist-packages/xen/remus/vm.py (seen using pdb) :
27 if domid:
28 try:
29 self.dominfo = server.xend.domain(domid, ''all'')
> /usr/local/lib/python2.7/dist-packages/xen/remus/vm.py(36)__init_...

2017 Jun 06

3

integrating 2 lists and a data frame in R

Hi Bogdan,
Kinda messy, but:
N <- data.frame(N=c("n1","n2","n3","n4"))
M <- data.frame(M=c("m1","m2","m3","m4","m5"))
C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
MN<-as.data.frame(m...

2000 Aug 30

3

family question

...ed such that each family stops having
children when more boys than girls are present ? Equal numbers of boys and
girls are got in the population, according to the simulation, is that obvious ?
ND <- NP <- NULL #ND - number boys, P: girls
for (j in 1:1000000) # number of families
{
n1 <- n2 <- NULL
for (i in 1:20000) {
n1[i] <- rbinom(1,1,0.5) # each equally likely - here number of boys
n2[i] <- 1- n1[i] # and girls
if(sum(n1)>sum(n2)) break
}
ND[j] <- sum(n1)
NP[j] <- sum(n2)
}
sum(ND)/sum(NP)
j
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

2017 Jun 06

1

integrating 2 lists and a data frame in R

> On Jun 6, 2017, at 4:01 AM, Jim Lemon <drjimlemon at gmail.com> wrote:
>
> Hi Bogdan,
> Kinda messy, but:
>
> N <- data.frame(N=c("n1","n2","n3","n4"))
> M <- data.frame(M=c("m1","m2","m3","m4","m5"))
> C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
> MN<-...

2017 Jun 06

1

integrating 2 lists and a data frame in R

...e, Jun 6, 2017 at 7:44 AM, David Winsemius <dwinsemius at comcast.net>
wrote:
>
> > On Jun 6, 2017, at 4:01 AM, Jim Lemon <drjimlemon at gmail.com> wrote:
> >
> > Hi Bogdan,
> > Kinda messy, but:
> >
> > N <- data.frame(N=c("n1","n2","n3","n4"))
> > M <- data.frame(M=c("m1","m2","m3","m4","m5"))
> > C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"),
> I=c(100,300,400...

2008 Oct 24

1

pwr.2p2n.test when the ratio of n1/n2 is known

Hi,
I am trying to do a power calculation for a difference in proportions
test where I want to estimate the sample size required. I know (well
estimate) that group one (n1) is 10% of the population and group 2 (n2)
is 90% of the population. I know the effect size (h). pwr.2p2n.test
only allows one variable to be left null whereas I would like both n1
and n2 to be determined where I know there relative proportions. Any
ideas how to do this? Is there a different function?
Thanks
Dan
--
****************...

2017 Jun 06

1

integrating 2 lists and a data frame in R

...emius <dwinsemius at comcast.net>
> wrote:
>
>>
>> > On Jun 6, 2017, at 4:01 AM, Jim Lemon <drjimlemon at gmail.com> wrote:
>> >
>> > Hi Bogdan,
>> > Kinda messy, but:
>> >
>> > N <- data.frame(N=c("n1","n2","n3","n4"))
>> > M <- data.frame(M=c("m1","m2","m3","m4","m5"))
>> > C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"),
>> I=c...

2008 Dec 10

10

AFR healing problem after returning one node.

I''ve got configuration which in simple includes combination of afrs and
unify - servers exports n[1-3]-brick[12] and n[1-3]-ns and client got
cluster configuration:
volume afr-ns
type cluster/afr
subvolumes n1-ns n2-ns n3-ns
option data-self-heal on
option metadata-self-heal on
option entry-self-heal on
end-volume
volume afr1
type cluster/afr
subvolumes n1-brick2 n2-brick1
option data-self-heal on
option metadata-self-heal on
option entry-self-heal on
end-volume
volume afr...

2017 Jun 06

1

integrating 2 lists and a data frame in R

Here's another approach:
N <- data.frame(N=c("n1","n2","n3","n4"))
M <- data.frame(M=c("m1","m2","m3","m4","m5"))
C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
# Rebuild the factors...

2011 Sep 14

1

ruby to solve a physics question

I am trying to solve one of my graduate level physics problems with
ruby...
Here is what I have so far...
a6=0.0
for n1 in -10..10
for n2 in -10..10
for n3 in -10..10
if n1!=0 and n2!=0 and n3!=0
p=Math.sqrt(n1**2+n2**2+n3**2+n1*n2/1.414+n1*n3/1.41+n2*n3/1.414)
a6+=(1/p)**6
end
end
end
end
puts a6
What I''ve got here is a 10x10x10 face-centered cubic lattice of atoms....

2011 Oct 06

2

Wilcox Test / Mann Whitney U Test

...rep(1,21),rep(4,3))
tab2 = rbind(table(dropouts),table(remain))
ordTest = wilcox.test(x=dropouts,y=remain,correct=FALSE,exact=FALSE)
cumsum(colSums(tab2))
W = max(c(sum(rank(cbind(dropouts,remain))[1:length(dropouts)]),sum(rank(cbind(dropouts,remain))[-(1:length(dropouts))])))
n1 = length(dropouts)
n2 = length(remain)
testStat = (S-n1*(n1+n2+1)/2)/(sqrt(n1*n2*(n1+n2+1)/12))
2*(1-pnorm(testStat))

1997 Jul 09

1

R-beta: Problem with `rpois'

There is a problem with `rpois'. It does seem to take care about the
order of the arguments. This is an example:
> rpois(n=1,lambda=2)
[1] 3
> rpois(lambda=2,n=1)
[1] 2 0
It obviously uses the first argument as the number of samples to be
drawn, which is wrong.
I used Version 0.49 Beta (April 23, 1997).
Fredrik

1997 Jul 09

1

R-beta: Problem with `rpois'

There is a problem with `rpois'. It does seem to take care about the
order of the arguments. This is an example:
> rpois(n=1,lambda=2)
[1] 3
> rpois(lambda=2,n=1)
[1] 2 0
It obviously uses the first argument as the number of samples to be
drawn, which is wrong.
I used Version 0.49 Beta (April 23, 1997).
Fredrik

1998 Nov 16

0

Re: Hotelling corrected

By accident, I left out the lines defining n1 and n2. Here it is as a
function.
Peter B.
hotelling <- function(d1,d2){
k <- ncol(d1)
n1 <- nrow(d1)
n2 <- nrow(d2)
xbar1 <- apply(d1,2,mean)
xbar2 <- apply(d2,2,mean)
dbar <- xbar2-xbar1
v <- ((n1-1)*var(d1)+(n2-1)*var(d2))/(n1+n2-2)
t2 <-...

2017 Jun 06

1

integrating 2 lists and a data frame in R

Dear Bert,
thank you for your response. here it is the piece of R code : given 3 data
frames below ---
N <- data.frame(N=c("n1","n2","n3","n4"))
M <- data.frame(M=c("m1","m2","m3","m4","m5"))
C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
how shall I integra...

1998 Jul 01

4

R-beta: R-0.62.1 under Digital Unix

I am grateful for the advice of Douglas Bates on my earlier problem in making
R-0.62.1, but I'm afraid I'm still having problems....
I have been installing the various updates to R for quite some time on my
alpha, and it is only now that I have been having really severe problems.
The three or 4 versions before 0.61.1 installed without error. For 0.61.1 I
needed to install GNU make.
For

2002 Mar 01

5

Power of t-test in R vs. S-PLUS

...ded")
Two-sample t test power calculation
n = 5
delta = 0.5488882
sd = 0.2
sig.level = 0.01
power = 0.8
alternative = two.sided
NOTE: n is number in *each* group
S-PLUS 2000 Professional Release 2:
> normal.sample.size(n1=5, n2=5, mean=0, sd1=0.2, sd2=0.2, power=0.8,
alpha=0.01, alternative="two.sided")
mean1 sd1 mean2 sd2 delta alpha power n1 n2 prop.n2
1 0 0.2 0.4322771 0.2 0.4322771 0.01 0.8 5 5 1
------------------------------------------
Dr. J?rn Quedenau
Coordinator Data Managemen...