Displaying 10 results from an estimated 10 matches for "pr4".
Did you mean:
cr4
2009 Jan 28
3
initial value in 'vmmin' is not finite
...rbind(x[5],x[6])+ rho=exp(x[7])/(1+exp(x[7]))+ ey = exp((yop%*%con1)/rho)+ ed = exp((dan%*%con2)/rho)+ eh = exp((hil%*%con3)/rho)+ ew = exp((wt%*%con4)/rho)+ ev = ey+ed+eh+ew+ den=(ey+ed+eh+ew)+ iv = rho*log(den)+ pp=exp(x[4]+iv)/(1+exp(x[4]+iv))+ pr1 =pp* ey/den+ pr2 =pp* ed/den+ pr3 =pp* eh/den+ pr4 =pp* ew/den+ pnp=1/(1+exp(x[4]+iv))+ likelihood = (pnp*yogurt[,1])+(pr1*yogurt[,2])+(pr2*yogurt[,3])+(pr3*yogurt[,4])+(pr4*yogurt[,4])+ lsum = log(likelihood)+ return(-colSums(lsum))+ }> p = optim(c(0,0,0,0,0.1,-2,-0.2),fr, hessian = TRUE, method = "BFGS")Error in optim(c(0, 0, 0, 0, 0...
2009 Mar 02
1
initial gradient and vmmin not finite
...)
con4 = rbind(x[5],x[6])
rho = exp(x[7])/(1+exp(x[7]))
ey = exp((yop%*%con1)/rho)
ed = exp((dan%*%con2)/rho)
eh = exp((hil%*%con3)/rho)
ew = exp((wt%*%con4)/rho)
ev = ey+ed+eh+ew
den=(ey+ed+eh+ew)
iv = rho*log(den)
pp=exp(x[4]+iv)/(1+exp(x[4]+iv))
pr1 =pp*(ey/den)
pr2 =pp*(ed/den)
pr3 =pp*(eh/den)
pr4 =pp*(ew/den)
pnp=1/(1+exp(x[4]+iv))
likelihood = (pnp*yogurt[,1])+(pr1*yogurt[,2])+(pr2*yogurt[,3])+(pr3*yogurt[,4])+(pr4*yogurt[,4])
lsum = log(likelihood)
return(colSums(lsum))
}
p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2)))
p
____________________________________________________________...
2012 Jul 17
2
Problem creation tensor
...em[,j],levels=1:5)
}
occ=factor(c(rep(1,dim(M1)[1]),rep(2,dim(M2)[1]),
...,rep(21,dim(M21)[1])),levels=1:21)
Z <- array(rep(0,m*21*5),c(m,21,5),
list(paste("item",1:m,sep=""),paste("Occ",1:21,sep=""),c("pr1","pr2","pr3","pr4","pr5")))
for ( i in 1:m){
as.matrix(table(occ, UserItem[,2]))
Z[i,,]=table(occ, UserItem[,i])
}
Z.CAND <- CANPARA(Z,dim=7)
I have implemented this code but I have one error in correspondance of:
for ( i in 1:m){
Z[i,,]=table(occ,UserItem[,i])
}
and error is:...
2005 Jun 27
1
A. Mani : colours in Silhouette
Hello,
In cluster analysis with cluster, how does one colour the silhouette
plots ? For example in using pam. There seems to be some problem there.
Everything else can be coloured.
Thanks,
A. Mani
Member, Cal. Math. Soc
2005 Jun 28
0
R-help Digest, Vol 28, Issue 28
...= NULL, sub = NULL,
>
> xlab = expression("Silhouette width " * s[i]), col = "gray",
> do.col.sort = length(col) > 1, border = 0, cex.names = par("cex.axis"),
> do.n.k = TRUE, do.clus.stat = TRUE, ...)
>
>
> data(ruspini)
> pr4 <- pam(ruspini, 4)
> si <- silhouette(pr4)
> plot(si,col = "red")
>
I tried that before with many more options and got a blank image. It must have
been due to the options.
> The issue is that whenever code is written there is always a choice as to
> what...
2011 Apr 09
1
loop and sapply problem, help need
...this question
M1 <- 1:10
lcd1 <- c(11, 22, 33, 44, 11, 22, 33, 33, 22, 11)
lcd2 <- c(22, 11, 44, 11, 33, 11, 22, 22, 11, 22)
lcd3 <- c(12, 12, 34, 14, 13, 12, 23, 23, 12, 12)
#generating variables through sampling
pvec <- c("PR1", "PR2", "PR3", "PR4", "PR5", "PR6", "PR7", "PR8", "PR9",
"PR10")
fun11 <- function(x){
smpool <- c(1,2,3,3)
x <- sample(smpool, 10, replace= TRUE)
}
newd <- sapply (pvec, fun11)
# function to recode this new generated data
fun3 &l...
2007 Jul 11
1
"Exclude" test fails on OS X (i386 Darwin 8.10.1)
Version 3.0.0cvs protocol version 30.PR4
Capabilities:
64-bit files, 32-bit inums, 32-bit timestamps, 64-bit long ints,
socketpairs, hardlinks, symlinks, IPv6, batchfiles, inplace,
append, no ACLs, xattrs, iconv
The exclude test fails on OS X Darwin Kernel Version 8.10.1
RELEASE_I386. It just hangs there until CTL-C. The same...
2003 Jun 17
1
lme() vs aov(y ~ A*B + Error(aa %in% A + bb %in% B)) [repost]
I've posted the following to R-help on May 15.
It has reproducible R code for real data -- and a real
(academic, i.e unpaid) consultion background.
I'd be glad for some insight here, mainly not for myself.
In the mean time, we've learned that it is to be expected for
anova(*, "marginal") to be contrast dependent, but still are
glad for advice if you have experience.
Thank
2009 Jan 28
2
t.test in a loop
Hi All,
I've been having a little trouble with creating a loop that will run a a
series of t.tests for inspection,
Below is the code i've tried, and some checks i've looked at.
I've used the get(paste()) idea as i was told previously that the use of the
eval should try and be avoided.
I've run a single syntax to check that my systax is correct and works
without any problems
2008 Jun 30
4
Rebuild of kernel 2.6.9-67.0.20.EL failure
Hello list.
I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without
modifications.
How did I try it?
Created a (non-root) build environment (not a mock )
Installed the kernel.scr.rpm and did a
rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee
prep-out.log
The build failed at the end:
Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL
Checking