[EMAIL PROTECTED] wrote:
I'm having trouble with the standard deviation distribution
as shown on http://mathworld.wolfram.com/StandardDeviationDistribution.html .
(Eric Weisstein references Kenney and Keeping 1951, which I can't check.)
I believe the graphs they show, but when I code the
Dear R-helpers:
thank your for your attention.
i am a newer to R and i am doing some protein category classification based on
the amino acid sequence.while i have some questions urgently.
1. any packages for analysis amino acid sequence
2. given two sequences AAA and BBB,how can i combine them
Hello ,
I am hoping for some advice. I want to run a repeated measures ANOVA. The
primary problem is that my attempt to remove missing values created a
dataset of missing values.
The data set consists of 92 rows (1 row per participant) x 186 variables.
The steps of the analysis undertaken are
Dear R-helpers:
thank your for your attention.
i am a newer to R and i am doing some protein category classification based on
the amino acid sequence.while i have some question urgent.
1. any packages for analysis amino acid sequence
2. given two sequences AAA and BBB,how can i combine them
I'm not aware of package for amino acid sequence analysis, but it'd
better to search yourself using, e.g., RSiteSearch() and help.search().
Regarding your question consider the following:
strg1 - c(AAA, AA, )
strg2 - c(BBB, , BBB)
strg - paste(strg1, strg2, sep = )
nA - nchar(gsub(B, ,
[EMAIL PROTECTED] wrote:
sddist
function(s,n) {
sig2 - n*s*s/(n-1)
2*(n/(2*sig2))^((n-1)/2) / gamma((n-1)/2) * exp(-n*s*s/(2*sig2)) *
s^(n-2)
}
There is another way using more statistical knowledge:
Building on functions already available in R
one may note that the easiest way of
Several weeks ago the idea of an 'R Community'
mailing list was discussed. I was impatiently awaiting
the announcement for it, so asked Martin Maechler
about it yesterday. I misinterpreted his message on
the subject, so perhaps others did as well.
Martin is seeking a couple of volunteers to
If I understand correctly Rick it trying to fit a model with random
effects on a binary response when there are either 1 or 2 observations
per group. I think that is very optimistic because there is so little
information available per random effect (exactly 1 or 2 bits of
information per random
On 6/16/06, Leaf Sun [EMAIL PROTECTED] wrote:
Sorry, I thought it was a straightforward question inside which I was stuck .
I used nls( ) to estimate a and b in this function.
nls(y~ a*x^b,start=list(a=a1,b=b1)
seems the start list I gave was not able to reach convergence and it gave
Hello all,
We are pleased to announce the official launch of the R Wiki at
http://wiki.r-project.org. Although there are already many sources of R
documentation, this R Wiki is a complementary tool, in the sense that
users become R documentation authors, a little bit like users become
I have some function to solve your problem.
#Plot stock quote
PSQ - function(INSTRUMENT,FROM,TO,START){
#Procedure to download Stock Quote Information
DATA -yahoo.get.hist.quote(
instrument = INSTRUMENT,
destfile = paste(INSTRUMENT, .csv, sep = ),
start=FROM,
Hi all,
My question concerns 2 error messages; one in the gam package and one in
the mgcv package (see below). I have read help files and Chambers and
Hastie book but am failing to understand how I can solve this problem.
Could you please tell me what I must adjust so that the command does not
I have to the treatment effect base on the observational data.And the
treatment variable is multinomial rather than binary.Because the
treatment assignment is not random,so the selection-bias exists.Under
this condition,what's the best way to estimate the treatment effect?
I know that if the
Dears,
I saw in the R-project site some R-books. However, I'm new in this
community and I didn't figure out what are the best books.
Can you suggest me some reference books? My intentions with R is
concerned to Artificial Intelligence simulations, Classification and
general Statistics (e.g.,
On Fri, 16 Jun 2006, Frede Aakmann Tøgersen wrote:
Well Jia, you use 'all' as a name for your dataframe, but this is also a
function, see ?all. If I try it with
mydata - data.frame(z1,z2,z3)
all goes well.
Exactly: 'step' is in the stats namespace and so object 'all' in the base
namespace
On Wed, 14 Jun 2006, toby marthews wrote:
I'm trying to map some time series data on to dates and
because I'm using R's strptime' facility I get an hour's
break in my time series: my readings for between 1am and
2am on 27th March 2005 won't map onto anything (I'm using a
data logger to
KURT: In case you have time to look briefly at the 'garch' code, you
might be interested in this example. This email includes a
self-contained example that generates a warning that seems to me
unnecessary: Warning in sqrt(pred$e) : NaNs produced. This occurs
because pred$e[1] = -0.121.
Dear Anderson,
I'd highly recommend Modern Applied Statistics with S (4th Ed., by WN
Venables BD Ripley). This has lots and lots of material on the topics of
your interest.
Best wishes
Paul
--
View this message in context: http://www.nabble.com/Books-t1807065.html#a4926486
Sent from the R
I don't know of R code per se, but
Check papers in the econometric literature by Heckman.
And another good recent reference: Stat. Med. 2005 Jul 30;24(14):2197-212.
On the use of discrete choice models for causal inference
Rusty Tchernis1;?; ?, Marcela Horvitz-Lennon2; 3 and Sharon-Lise T.
Dear WizaRds,
having met some of you in person in Vienna, I think even more fondly
of this community and hope to continue on this route. It was great
talking with you and learning from you. Thank you. I am trying to work
through an artificial example in post stratification. This is my
[EMAIL PROTECTED] wrote:
I'm trying to build forecasting equation from weights of 2-2-1 neural net.
Running the nnet function gives me a vecto of 9 weights, but I don't know how
to build the equation form these values.
Can anyone advice? Or at least tell me where the nnet output is described
Le Lundi 12 Juin 2006 06:51, Lorenzo Isella a écrit :
Dear All,
A simple question: packages like fitdistr should be ideal to analyze
samples of data taken from a univariate distribution, but what if
rather than the raw data of the observations you are given directly
and only a histogram?
this is not really an R question but I imagine I would
use R to do it.
I have a loop that has thousands of iterations and
after each iteration the results are a
box leung pvalue
a regression t-stat
and a profit number
so, i have 1000's of these sets of the 3 numbers above.
obviously, the
From: [EMAIL PROTECTED]
Date: Sun Jun 18 17:47:39 CDT 2006
To: r-help@stat.math.ethz.ch
Subject: [R] more of a statistics/display question
one more piece of info in regard to below :
it's okay if i have turn the profit number
into a 0/1 indicator variable with 0 indicating
negative and 1
On Sunday 18 June 2006 08:27, Anderson de Rezende Rocha wrote:
Dears,
I saw in the R-project site some R-books. However, I'm new in this
community and I didn't figure out what are the best books.
Can you suggest me some reference books? My intentions with R is
concerned to Artificial
Won't 'interval censoring using 'survreg' in the 'survival' package
handle this?
(http://finzi.psych.upenn.edu/R/library/survival/html/survreg.html)
Hope this helps.
Spencer Graves
Vincent Goulet wrote:
Le Lundi 12 Juin 2006 06:51, Lorenzo Isella a écrit :
Dear
Hi list,
Is any body has implemented the projection pursuit regression (ppr)
method in MASS library.Running the ppr method is not a problem , But problem is
how we are getting the model parameters using least square method ,becoz the
function is looking like a continous function
Hi:
I would like to create a custom command in R which generates SQL, which is
then processed via RODBC.
For example, the user would type:
retrieve firstName('JOHN') middlleName('WILLIAMS') lastName('FORD')
This would generate the following SQL which would then be processed by
RODBC:
Hello.
I am having some trouble using nlm in R for windows version 2.2.1. I have
to minimise a function which includes in its insides the evaluation of
integrals through the adapt function.
In order to make R work less, I calculated the gradient and coded it. I
made sure that the code and
This generates the sql statement so just pass that to your database:
retrieve - function(...) {
args - list(...)
sub(and, select * from people where,
paste(rbind(and, names(args), =, dQuote(args)), collapse = ))
}
# test
retrieve(firstname = JOHN, middlename = WILLIANS, lastname =
30 matches
Mail list logo