Hi,
I have this string, in which I want to extract some of it's element:
x - Best-K Gene 11340 211952_at RANBP5 Noc= 3 - 2 LL= -963.669 -965.35
yielding this array
[1] 211952_at RANBP5 2
In Perl we would do it this way:
__BEGIN__
my @needed =();
my $str = Best-K Gene 11340 211952_at
Dear all,
Currently I have the following output
mydensity - density(x)
print(mydensity)
x y
Min. : -92.14 Min. :0.000e+00
1st Qu.: 356.66 1st Qu.:5.530e-09
Median : 805.45 Median :4.681e-05
Mean : 805.45 Mean :5.564e-04
3rd Qu.:1254.24 3rd
Hi,
Is there an R function or package that computes
p-value?
-GV
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented,
Hi,
I have two histograms created separately using
the following code. It creates two separate figures.
dat - read.table(file=GDS1096.modelout, head = FALSE )
__BEGIN__
dat - read.table(file=GDS1096.modelout, head = FALSE )
hist(dat$V2, main=AIC Freq, xlab = \# Component, breaks = 36, xlim =
Dear all,
I am trying to model number of samples from
a given series. The series are modelled according
Gamma function.
In order to estimate the # samples, I use BIC/AIC
with MLE (computed from dgamma function).
Here is the code I have.
__BEGIN__
mlogl - function( x_func, theta_func, samp) {
Hi,
Below I have a function mlogl_k,
later it's called with nlm .
__BEGIN__
vsamples- c(14.7, 18.8, 14, 15.9, 9.7, 12.8)
mlogl_k - function( k_func, x_func, theta_func, samp) {
tot_mll - 0
for (comp in 1:k_func) {
curr_mll - (- sum(dgamma(samp, shape = x_func,
scale=theta_func,
, i would just
tell you what it was.
On Thu, May 22, 2008 at 9:36 PM, Charilaos Skiadas wrote:
On May 22, 2008, at 9:14 PM, Edward Wijaya wrote:
Hi,
Below I have a function mlogl_k,
later it's called with nlm .
__BEGIN__
vsamples- c(14.7, 18.8, 14, 15.9, 9.7, 12.8)
mlogl_k
Hi all,
What's the construct to import a built-in package
(e.g. ARIMA) into the R script.
I can't seem to find them in the documentation.
Regards,
Edward
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
You might try cron job under Windows.
http://drupal.org/node/31506
HTH.
- Edward
On Thu, May 22, 2008 at 8:51 AM, Thomas Pujol [EMAIL PROTECTED] wrote:
I am using R in a Windows environment.
I store my data in a Microsoft SQL database that gets updated automatically
nightly. Once my SQL
Hi,
I tried to compute maximum likelihood under gamma distribution,
using nlm function. The code is this:
__BEGIN__
vsamples- c(103.9, 88.5, 242.9, 206.6, 175.7, 164.4)
mlogl - function(alpha, x) {
if (length(alpha) 1) stop(alpha must be scalar)
if (alpha = 0) stop(alpha must be
and Computer Science
Hanover College
On May 21, 2008, at 10:05 PM, Edward Wijaya wrote:
Hi,
I tried to compute maximum likelihood under gamma distribution,
using nlm function. The code is this:
__BEGIN__
vsamples- c(103.9, 88.5, 242.9, 206.6, 175.7, 164.4)
mlogl - function(alpha, x
Hi,
Currently the R script I have is executed with this command:
$ R CMD BATCH mycode.R
And the output is stored in mycode.Rout.
Is there a way I can issue command from shell (like above)
so that the output is printed to STDOUT?
It's troublesome to open the Rout file every time to debug.
Dear all,
How can I compute the log likelihood of a gamma
distributions of a vector.
I tried the following. But it doesn't seem to work:
samples-c(6.1, 2.2, 14.9, 9.9, 24.6, 13.2)
llgm - dgamma(samples, scale=1, shape=2, log = TRUE)
It gives
[1] -4.291711 -1.411543 -12.198639 -7.607465
.
Regards,
Edward
Edward Wijaya 写道:
Dear all,
How can I compute the log likelihood of a gamma
distributions of a vector.
I tried the following. But it doesn't seem to work:
samples-c(6.1, 2.2, 14.9, 9.9, 24.6, 13.2)
llgm - dgamma(samples, scale=1, shape=2, log = TRUE)
It gives
[1
Hi,
How can I convert the matrices to list.
For example I have this snippet:
samples-mymatrix[1,]
print(samples)
which prints:
V1 V2V3V4V5V6
1 103.9 88.5 242.9 206.6 175.7 164.4
How can I convert the object samples such that it prints:
[1] 103.9 88.5 242.9 206.6 175.7
, [EMAIL PROTECTED] wrote:
unclass or should work , I think. so
nonames-unclass(samples)
On Tue, May 20, 2008 at 10:16 PM, Edward Wijaya wrote:
Hi,
How can I convert the matrices to list.
For example I have this snippet:
samples-mymatrix[1,]
print(samples)
which prints:
V1
Peter Alspach
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Edward Wijaya
Sent: Wednesday, 21 May 2008 2:17 p.m.
To: r-help@r-project.org
Subject: [R] Converting Data Types
Hi,
How can I convert the matrices to list.
For example I have
Hi all,
Is there any EM (Expectation Maximization) package available for R?
CRAN Package archive (http://cran.r-project.org/)
doesn't seem to have it.
Regards,
Edward
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
18 matches
Mail list logo