Displaying 12 results from an estimated 12 matches for "miu".
Did you mean:
mi
2013 Jan 13
1
How to combine two loops?
...tion is:
Stage I
I have 2 groups and for each group I need to compute the following steps;
i) Generate 3 random numbers from normal distribution and square them.
ii) Repeat step 1 for 15 times and at the end I will get 15 random numbers.
I already done stage I using for loop.
n1<-3
n2<-3
miu<-0
sd1<-1
sd2<-1
asim<-15
w<-rep(NA,asim)
x<-rep(NA,asim)
for (i in 1:asim)
{
print(i)
set.seed(i)
data1<-rnorm(n1,miu,sd1)
data2<-rnorm(n2,miu,sd2)
w[i]<-sum(data1^2)
x[i]<-sum(data2^2)
}
w
x
Second stage is;
Stage II
For each group, I need to:
i) Sort the group;
i...
2012 Jun 14
0
fixed trimmed mean for j-group
...<-(exp(g*b)-1)/g*(exp((h*b^2)/2))
y<-(exp(g*c)-1)/g*(exp((h*c^2)/2))
z<-(exp(g*d)-1)/g*(exp((h*d^2)/2))
######################## FIXED SYMMETRIC TRIMMED MEAN
###############################
asim<-5000
pv<-rep(NA, asim)
for(j in 1:asim)
{
print(j)
set.seed(j)
n1=15
n2=15
n3=15
n4=15
miu=0
sd1=1
sd2=1
sd3=1
sd4=1
a=rnorm(n1,miu,sd1)
b=rnorm(n2,miu,sd2)
c=rnorm(n3,miu,sd3)
d=rnorm(n4,miu,sd4)
## data transformation
g=0
h=0
w<-a*exp(h*a^2/2)
x<-b*exp(h*b^2/2)
y<-c*exp(h*c^2/2)
z<-d*exp(h*d^2/2)
mat1<-sort(w)
mat2<-sort(x)
mat3<-sort(y)
mat4<-so...
2012 Jul 07
0
fixed trimmed mean for group
...##############FIXED SYMMETRIC TRIMMED MEAN#############
> asim<-5000
> pv<-rep(NA, asim)
> for(j in 1:asim)
> {
> print(j)
> set.seed(j)
> n1=15
> n2=15
> n3=15
> n4=15
> miu=0
> sd1=1
> sd2=1
> sd3=1
> sd4=1
>
> a=rnorm(n1,miu,sd1)
> b=rnorm(n2,miu,sd2)
> c=rnorm(n3,miu,sd3)
> d=rnorm(n4,miu,sd4)
>
> ## data transformation
> g=0
> h=0
>
>...
2012 Dec 17
1
seeking a help on if function
...(which mean equal trimming), the result will come out as "a" not the value of ssdw when equal trimming is done. I f anyone could with this, I really appreciate it. I need the value of ssdw so that I can proceed with computing the F test. Thanks.
Regards,
Hyo Min
UPM Malaysia
n1<-15
miu<-0
sd1<-1
data1<-rnorm(n1,miu,sd1)
## data transformation as G=0.5 and h=0.5 (data is normal)
G<-0.5
h<-0.5
w<-((exp(G*data1)-1)/G)*(exp((h*data1^2)/2))
group1<-sort(w)
g<-0.15 ## proportion of trimming
## group 1
uw<-floor(0.2*n1)
low1<-mean(group1[1:uw])
up...
2010 Feb 07
2
conditioned xyplot, many y variables
...)
JUNK$Species<-factor(JUNK$Species)
JUNK$ID<-factor(JUNK$ID)
summary(JUNK)
# Traditional graphics solution:
par(mfrow=c(3,2),mar=c(0,0,0,0)+0.0,oma=c(4,4,4,1),xpd=FALSE, las=0)
varNamesAndLabels<-data.frame(
name=c("z", "QQQ", "ppp")
, label=c("z (mIU/mL)", "QQQ (pg/L)", "ppp (mg/L)")
)
rownames( varNamesAndLabels)<- varNamesAndLabels$name
count_y_variables<-0
for(this_y_name in rownames( varNamesAndLabels) ) {
count_y_variables <- count_y_variables + 1
countSpecies<-0
for(thisSpecies in levels(JUNK$Spe...
2010 Dec 08
1
the output of function lars
..." "actions" "entry" "Gamrat"
[11] "arc.length" "Gram" "beta" "mu" "normx"
[16] "meanx"
I do not know what does RSS, arc.length, mu mean???
If mu means \hat{miu}, it does not seems to be right. I did not simulate my
data with mean.
besides I caculated R2 by myself but it cannot match the R2here
I think R2=1-sum((y-x%*%lobj$beta)^2)/var(y)*19 #I have 20 observations...my
R2 is lower than the one here
Thanks a lot
--
View this message in context: http:...
2012 Jul 03
2
EM algorithm to find MLE of coeff in mixed effects model
...sum(sapply(data.list.wob,loglike))
ECM=F
for (m in 1:300)
{
Sig.hat=Zi%*%psi.old%*%t(Zi)+var.old*diag(nrow(Zi))
W.hat=psi.old-psi.old%*%t(Zi)%*%solve(Sig.hat)%*%Zi%*%psi.old
Sigb=psi.old-psi.old%*%t(Zi)%*%solve(Sig.hat)%*%Zi%*%psi.old
det(Sigb)^(-0.5)
Y.minus.X.beta=t(t(y.m)-c(Xi%*%beta.old))
miu.m=t(apply(Y.minus.X.beta,MARGIN=1,function(s,B=psi.old%*%t(Zi)%*%solve(Sig.hat))
{
B%*%s
}
)) ### each row is the miu_i
tmp1=permutations(length(gausspar$nodes),2,repeats.allowed=T)
tmp2=c(tmp1)
a.mat=matrix(gausspar$nodes[tmp2],nrow=nrow(tmp1)) #a1,a1
#a1,a2...
2012 Jul 03
0
need help EM algorithm to find MLE of coeff in mixed effects model
...sum(sapply(data.list.wob,loglike))
ECM=F
for (m in 1:300)
{
Sig.hat=Zi%*%psi.old%*%t(Zi)+var.old*diag(nrow(Zi))
W.hat=psi.old-psi.old%*%t(Zi)%*%solve(Sig.hat)%*%Zi%*%psi.old
Sigb=psi.old-psi.old%*%t(Zi)%*%solve(Sig.hat)%*%Zi%*%psi.old
det(Sigb)^(-0.5)
Y.minus.X.beta=t(t(y.m)-c(Xi%*%beta.old))
miu.m=t(apply(Y.minus.X.beta,MARGIN=1,function(s,B=psi.old%*%t(Zi)%*%solve(Sig.hat))
{
B%*%s
}
)) ### each row is the miu_i
tmp1=permutatio...
2013 Apr 02
0
Question about removing missing data rows in the qplot
Hi,
I've a question about the qplot. I use facets=Antibiotics~Gram.Stainng to draw 6 plots in a panel. Now the problem is within each 'Gram.Staining' categorical level there are some antibiotics does not have value (in my case is log(MIU)). Because some antibiotics are positive gram staining and some are negative.
So my question is how to remove these missing value rows in the plot?
My code is
qplot(Bacteria,log(MIC),data=project2,facets=Gram.Staining~Antibiotics,color=Bacteria,geom='bar',stat='identity',fill...
2020 Feb 27
2
[PATCH] Update the 5 year logo to 10 year logo
...IJ43B-515
zSTg%&L{n=>tp=7Hno_!W&^L;bcrX*$bc1qDX%B9XwA})o6I0R<6lcOs!pvMd>)YKR
zxg#@Sn~I8!Jq$=`Lx5QQ{NG5g*^USpwLZVS{YMl6!zuA%1$822?X@JC^tZRxc$rc?
zxwsiqMa!<$6GYPd>C`XXFX}pa;&9mA)h7uRy<GQBLegOt18$t>VhWVFv|z~qpX^A9
zA9+8Ah5uVn&1D5|mm_?=5oUQG#XF_qt!KwIt@$Ejoy)DhE(LJf@miU{z@NXX@L8Qr
z7^`y{59c_{*!49T^T)D(PU||tm<xstR|PS<VyoCLsM5ShOH$Xtl6y5!-@W3*9!h!O
z{JE`D!ecZwSTOjcMkn<xV*OQm@sAA!^6ez~gK35gX_bI(REpGP2z8ZT=^z<MO#z!2
zTDy{x^U&VU<_`}Fix(XqfDH}GtSB;$aH^xFPh?N^AD#P@m#HKyZ`DEU`%SOcd`*yh
zar0{@{P0QJ4*C+r6)YAkW}xNNUxpVZa=;OIG7|{`Wca_y;0-e}8CHgM`Mk|GSMM?8
z...
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
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...H#4NAvx>Eecbi%*ei*j*=i2UGkd}V)VR>lyM^q=#;(Meh3hE)%
z#nl27(jB$?yTthyk5n9}@4LNO7mPu(Low)M43cmlC}l^2hePo~?|+Pmo%+jQ&5Hro
zZ{B`zS~8C*<<<1sR}*Xd9(aBfMIYTF))x(sgq~_n=yq(0pWyjh`B4d;sf)k2Yu}EJ
z7C#DJUu0MHiMsgklkLUDvg3+^ArF%k!zZ7AJ$NVZUy-+Q at wq9%9}3t09$yN&m3ZYw
ztr5kKaZUfl`WYBA|J04)DMIUE#jIcNO$Ue9EIzz_s^`7h_Wrj{dlI!t<LgGlZl9rU
zRf6XE0pFl*Z=!CK#AhQ1o#z^kiM95g%T={FF4^s>x1$;#exfLL$QJe3@&Uis)58aI
zy&lmIjxcdmZ8J+!c~y(!`I)7{oIWWM>*4LQp51$^@jz<K<bx2SdkCZ6kPt_Nh!Cd`
zA&j!(zCqcrX7x%~(-$_I)P$PdikjUhyPQTeAS#t0D&_b$SiLpp9N&Hzwf2t#x*ah&
zCl;Wpb)#Wb_sIrgfd...