Displaying 6 results from an estimated 6 matches for "gexp".
Did you mean:
exp
2008 Jun 23
2
Pairwise Partitioning of a Vector
...0.1 70.0 73.0 75.0 83.9 93.1 97.6 98.8 113.9
PAIR2
part1 = 30.9 60.1
part2 = 70.0 73.0 75.0 83.9 93.1 97.6 98.8 113.9
....
PAIR9
part1 = 30.9 60.1 70.0 73.0 75.0 83.9 93.1 97.6 98.8
part2 = 113.9
I'm stuck with this kind of loop:
__BEGIN__
# gexp is a Vector
process_two_partition <- function(gexp) {
sort.gexp <- sort(as.matrix(gexp))
print(sort.gexp)
for (posb in 1:ncol(gexp)) {
for (pose in 1:ncol(gexp)) {
sp_b <- pose+1
sp_e <- ncol(gexp)
# This two doesn't do what I wa...
2008 Jun 12
1
Data.matrix fail to convert data.frame into matrix
...atrix function doesn't
seem to work.
__ BEGIN__
dat <- read.table("mydata", comment.char = "!" , na.strings = "null");
# Select n-genes by random sample
# n = 1
nosamp <- 1
geneid <- sequence(nrow(dat))
geneid.samp <- sample(geneid,nosamp)
geneid.samp
gexp<- dat[geneid.samp,]
gexp.arr <- data.matrix(gexp, rownames.force = NA)
print(is.matrix(gexp.arr))
print(gexp.arr)
__END__
Yielding this output:
__BEGIN__
> print(is.matrix(gexp.arr))
[1] TRUE
> print(gexp.arr)
V1 V2 V3 V4 V5 V6 V7 V8
10354 803.1 11...
2008 Jul 01
1
Plotting Bi-Gamma Distribution
...ve.
Notice that red line can correctly fit the histogram that has two peaks
(i.e. red curve also has two peaks).
But the gamma curve there only has one curve.
Is there a way I can fit the gamma function such that it also
yields two peaks?
Below is my code that generate the graph:
__BEGIN__
dgexp <- density(gexp.arr)
gexp.ymax <- max(dgexp$y)
# Histograms
gehist <- hist(gexp.arr,
col="blue",
main=paste(genm,"-",desc, ", Mean=",ge_mean, ", SD=",ge_sd, ",
Krt=",ge_kurt, &q...
2012 Feb 22
3
gamlss results for EXP and LNO seem to have reversed AIC scores
Hi,
I'm a bit puzzled by the gamlss fitting of exponential and lognormal data.
Gamlss seems to think that exponentially distributed data fits better with a
lognormal distribution, and vice versa.
For example,
X <- rexp(1000)
Gexp <- gamlss(X~1,family=EXP) # X~1 is X tilde 1
GAMLSS-RS iteration 1: Global Deviance = 2037.825
GAMLSS-RS iteration 2: Global Deviance = 2037.825
Glno <- gamlss(X~1,family=LNO)
GAMLSS-RS iteration 1: Global Deviance = 2185.763
GAMLSS-RS iteration 2: Global Deviance = 2185.763
Gexp$ai...
2008 Jun 26
0
Unsqueezing Stacked Figures
...es are compressed)
Is there a way I can resize the figure?
So that it can show proportional and
complete plots?
The code I have is the following:
(I can also provide the figure if needed).
__BEGIN__
trellis.device("png", color=TRUE)
png(out_fname)
par(mfrow = c(2,1))
scat <- plot(gexp.arr,
type = "l",
col="red",
xlab="Sample",
ylab="Exp Lvl",
main="My Plot",
split=c(1,1,1,2), more = TRUE
)
gehist <- hist(gexp.arr,
col="blue&qu...
2012 Dec 07
1
Error using constrOptim in constraint definition
Hello,
I'm trying to run constrOptim. It returns to me an error about the fact that constraints arguments (ui and ci) are non compatibles:
> optout= constrOptim(startparams, f=ImpulseSS, grad=grImpulse, ui=UI, ci=CI, data=gexp[k,], t=t)
Error in ui %*% theta : non-conformable arguments
I would like to point out that I can calculate that product in the command line:
> UI %*% startparams
[,1]
[1,] 2.297217e+05
[2,] 1.699765e+05
[3,] 2.037806e+05
[4,] 6.535699e+00
[5,] 1.014253e+01
[6,] 4.21...