Hi, are the methods for multiple testing p value adjustment (Shaffer, Westfall,
free) implemented in the function adjusted() in multcomp package so called
closed testing procedure? what about those methods (holm, hochberg, hommel, BH,
BY) implemented in the p.adjust() in the stats package?
Hi,
I have a question regarding the modeling methodology of the following problem:
* I have two data sets {X_i,y_i} {X_i,z_i}, i=1..N,
where y_i = f(X_i) + i.i.d. Gaussian noise
and z_i = g(X_i) + i.i.d. Gaussian noise
* I apply bayesian linear regression to each of them and obtain
I am not a Bayesian. In the non-Bayesian case you would use SUR to model both
equations simultaneously. If both use the exact same matrix of data, X
(i.e., the value are numerically absolutely identical), then SUR will
collapse to OLS. In that sense you get a combined estimate using SUR that
Hi all..
I'm very new R, and i'm analyzing microarray data using Bioconductor..
Recently i was given microarray data to analyze. The problem is whenever i
run MAS5 presence calls algorithm,
it throws an error saying NA/NaN/Inf in foreign function. How do i remove
such NA/NaN/Inf's ??
I tried
Dear Igors
On 9 September 2011 22:00, Igors igors.lahanc...@gmail.com wrote:
Thank you for your answer. However it doesn't solve my problem fully.
The problem is that I have much bigger data set than I sent to you (it was
only a small part : 874 obs.). My full data set is 546718 obs.
If I
On Sep 11, 2011, at 5:38 AM, anand m t wrote:
Hi all..
I'm very new R, and i'm analyzing microarray data using Bioconductor..
Recently i was given microarray data to analyze. The problem is
whenever i
run MAS5 presence calls algorithm,
it throws an error saying NA/NaN/Inf in foreign
What R packages are available for performing classification tasks?
That is, when the predictor has done its job on the dataset (based on
the training set and a range of variables), feedback about the true
label will be available and this information should be integrated for
the next classification
Hello all:
As everyone knows the stock quotes in Yahoo can be downloaded directly
on EXCEL spreadsheet.
we can also use the simple sentence of package tseries of R to download
Yahoo quote:
x - get.hist.quote(instrument = ibm, start = 2010-01-01,
quote = close)
On Sep 11, 2011, at 11:42 AM, Jay wrote:
What R packages are available for performing classification tasks?
That is, when the predictor has done its job on the dataset (based on
the training set and a range of variables), feedback about the true
label will be available and this information
require(quantmod)
getSymbols('IBM', from='2010-01-01')
write.csv(Cl(IBM), file='ibm-2010-present.csv')
I can open the resulting file in Excel without problems. Depending on
your preference, you may want to use write.csv2 instead.
On 11 September 2011 06:59, Yumin zpx...@gmail.com wrote:
The
Hi,
I used the rseek search engine to look for suitable solutions, however
as I was unable to find anything useful, I'm asking for help.
Anybody have experience with these kinds of problems? I looked into
dynaTree, but as information is a bit scares and as I understand it,
it might not be what
Try this:
http://cran.r-project.org/web/views/MachineLearning.html
On 09/11/2011 12:43 PM, Jay wrote:
Hi,
I used the rseek search engine to look for suitable solutions, however
as I was unable to find anything useful, I'm asking for help.
Anybody have experience with these kinds of problems?
Try:
library(Hmisc)
?na.delete
Ken Hutchison
On Sep 11, 2554 BE, at 5:38 AM, anand m t anandro...@gmail.com wrote:
Hi all..
I'm very new R, and i'm analyzing microarray data using Bioconductor..
Recently i was given microarray data to analyze. The problem is whenever i
run MAS5 presence
Dear R-help,
Apologies if this isn't exactly the right place for this question.
I am trying to run R in emacs using ESS. I have done this
successfully many times - it works right now on my home computer in
windows 7 and in ubuntu.
The problem is that when I start emacs the regular splash
Dear all,
Can anyone take a look at my program below?
There are two functions: f1 (lambda,z,p1) and f2(p1,cl, cu).
I fixed p1=0.15 for both functions. For any fixed value of lambda (between
0.01 and 0.99),
I solve f1(p1=0.15, lambda=lambda, z)=0 for the corresponding cl and cu
values.
Then I
Alex,
This query belongs on the ESS list
ess-h...@r-project.org
Please place any followup questions there.
You seem to have accepted the default location of R, so ESS should just
work.
The only line you need in ~/.emacs
is
(require 'ess-site)
The other lines you mention cause interference and
Bill,
Thanks for the improvment. I had not been concerned in looking at the
type; was primarily interested in the sizes. Will go back and update
the function.
On Thu, Sep 8, 2011 at 11:53 AM, William Dunlap wdun...@tibco.com wrote:
In my.ls() you ought to convert the pos argument
to an
### code ###
x=sample(LETTERS[1:26],100,T)
prob=function(y){
count=sum(x==y)
total=length(x)
count/total
}
### end ###
How do I quote the letters A,B,C,D ect, in order to make the prob function
return the weights of the desired Letter?
Thanks!
--
Please include a meaningful 'Subject' because these conversations are
archived, and serve as a valuable help resource for the community. I
don't believe I understand what you want.
Stephen
On 09/11/2011 03:32 PM, li li wrote:
Dear all,
Can anyone take a look at my program below?
There
You quote them in quotes. A
Also: bad scoping.
Michael Weylandt
On Sep 11, 2011, at 6:11 PM, casperyc caspe...@hotmail.co.uk wrote:
### code ###
x=sample(LETTERS[1:26],100,T)
prob=function(y){
count=sum(x==y)
total=length(x)
count/total
}
###
Hi casperyc,
Try
# option 1
x - sample(LETTERS[1:26], 100, TRUE)
prob - function(y){
count=sum(x==y)
total=length(x)
count/total
}
prob('A')
prob('B')
prob('Z')
# option 2
tx - prop.table(table(x))
tx
tx['A']
# option 2.1
foo - function(x, l = NULL){
Hello,
you would get more answers if you code had proper indentation and comments.
Also, please provide a meaningful topic. You should also explain how
this is an R question and not just a debug my code request. What are
are you trying to achieve? Which of the numerous variables you
declared
Dear List,
After including cluster() option the coxreg (from eha package)
produces results slightly different than that of coxph (from survival)
in the following time-dependent treatment effect calculation (example
is used just to make the point). Will appreciate any explaination /
comment.
Sorry: there was an error in the weight calculation, fixed version is
the following, but still the final estimates differ as explained in
the original email:
#
require(survival)
require(eha)
data(heart)
head(heart)
follow - heart$stop - heart$start
fit - glm(transplant
On Fri, Sep 9, 2011 at 9:38 PM, Bigelow, Seth sbige...@fs.fed.us wrote:
I wish to display a single-panel Lattice figure with grouped data and fitted
regression lines. I don't seem to be able to get the
individual regression lines to display, e.g.;
d - data.frame(q = rep(1:6, each=10), cc =
Hi,
I want to store the output of my program to a file. However, With subsequent
runs of my code, I'd like to append to the same log file.
Currently, I'm using:
outfile - file(log.txt, open=w)
cat(results, file=outfile)
This works, but will overwrite the log file each time.
Is there a
I have date data as a numeric and hourly data in 0 to 2300 hours in a dataframe.
d - rep(20110101,24)
h - seq(from = 0, to = 2300, by = 100)
df - data.frame(LST_DATE = d, LST_TIME = h, data = rnorm(24, 0, 1))
S - chron(dates. = as.character(df$LST_DATE), times. =
27 matches
Mail list logo
