Dear R-users,
In the function constrOptim there is an option to get an approximation
to the hessian of the surrogate function R at MLE by declaring
hessian=TRUE in the calls to the function optim. I would like to ask
if it is advisable to get an approximate hessian for the funcrion f as
follows:
Sorting is a wonderful topic! Especially because you can discuss different
fundamental ideas like brute force, divide and conquer, and questions of
efficiency, tradeoffs of space and time, etc. In this way sort.6 is
one of a
sequence of sorting algorithms. Each of them demonstrates a specific
Hello,
I'm a french student working on cDNA microarrays.
I decided to use the bioconductor packages and I have a problem with
normalization functions.
As everybody who manipulates microarray data know, there are some spots
on the slide which are to be removed. In the marrayLayout there is a slot
Good morning !
Today I found a strange, for my poor knowledge of R, behaviour of
'which' on a matrix:
HAL9000 str(cluster.matrix)
num [1:227, 1:6300] 2 2 2 2 2 2 2 2 2 2 ...
HAL9000 class(cluster.matrix)
[1] matrix
HAL9000 ase - cluster.matrix[1:5,1:5]
HAL9000 ase
[,1] [,2] [,3] [,4] [,5]
[EMAIL PROTECTED] wrote:
Good morning !
Today I found a strange, for my poor knowledge of R, behaviour of
'which' on a matrix:
HAL9000 str(cluster.matrix)
num [1:227, 1:6300] 2 2 2 2 2 2 2 2 2 2 ...
HAL9000 class(cluster.matrix)
[1] matrix
HAL9000 ase - cluster.matrix[1:5,1:5]
HAL9000 ase
Hi there fellow R-Users,
Can anyone recommend a good book on the theory and practice of applying
GARCH models.
Also, does any one know of any R related subject material in addition to
library(tseries).
Regards
Wayne
Dr Wayne R. Jones
Senior Statistician / Research Analyst
KSS Limited
St
Hi,
I have a dynamic programming problem with many state variables, let's
call them l, n, m, etc. The transition probabilities are originally
given in an array form, eg
transtition[l,m,n,ll,mm,nn]
give the probability of going from l,m,n to ll,mm,nn.
However, the numerical solution is best
Hi there fellow R-Users,
Can anyone recommend a good book on the theory and practice
of applying
GARCH models.
Hello Wayne,
* Campbell, John, Lo, Andrew W., MacKinlay, A. Craig, The
Econometrics of Financial Markets, 1996, Princeton, NJ:
Princeton University Press.
On Thu, Apr 01, 2004 at 05:11:31PM +0200, Rau, Roland wrote:
Hello,
-Original Message-
From: Jingky P. Lozano [SMTP:[EMAIL PROTECTED]
Sent: Thursday, April 01, 2004 2:15 PM
To: [EMAIL PROTECTED]
Subject:[R] Question on Data Simulation
Dear mailing list,
Hi Folks,
I am interested in extending the repertoire of link functions
in glm(Y~X, family=binomial(link=...)) to include a tan link:
eta = (4/pi)*tan(mu)
i.e. this link bears the same relation to the Cauchy distribution
as the probit link bears to the Gaussian. I'm interested in sage
advice
Dear Marc,
I may misunderstand your question, so perhaps this answer isn't what you're
looking for:
A short time ago, Henrik Bengtsson posted a function to r-help (search for
plotSymbols in the list archive) that displays all available symbols. I
think that it's fair to say that fewer than 25
On Fri, 2 Apr 2004 [EMAIL PROTECTED] wrote:
I am interested in extending the repertoire of link functions
in glm(Y~X, family=binomial(link=...)) to include a tan link:
eta = (4/pi)*tan(mu)
i.e. this link bears the same relation to the Cauchy distribution
as the probit link bears to
After making
E - eigen( something )
I would like to extract those eigenvectors which have an eigenvalue of
1. If I had an isone() function, I would simply say
E$vectors[,which(isone(E))]
but the problem is that I have no such thing. I found all.equal, so I
could test for all.equal(x, 1), but
This is a question that I should have sent to 'ess-help', but I
take my chances...
In 'NEWS' (1.9.0) I read the good news
o Underscore '_' is now allowed in syntactically valid names, and
make.names() no longer changes underscores. Very old code
...
but when I try it in emacs
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hello all -
As I progress in R I am trying to automate functions I would have
normally farmed out to Excel, SPSS or Statistica. Single factor anova
is one of them. For example, a dataset from NIST StRD
Has anyone achieved a current successful build of R
1.8 or 1.9 or earlier for HPUX Tru 64 OSF 5.1?
There is no binary version of R for HPUX Tru 64 OSF
5.1. The R admin manual mentions that the native make
fails on Alpha/OSF (aka Tru 64) and gnu make must be
used instead. There are problems
GB == Göran Broström [EMAIL PROTECTED]
on Fri, 2 Apr 2004 16:29:00 +0200 writes:
GB This is a question that I should have sent to 'ess-help', but I
GB take my chances...
(yes, you should have..)
GB In 'NEWS' (1.9.0) I read the good news
GB o Underscore '_' is now
Hello,
This message is to announce that the R Reference Manual is now
available in book form.
There are two volumes, which cover all the commands in the base
package. They are available for order from all major bookstores.
The R Reference Manual - Base Package (2 volumes)
Hardback
On Fri, 2 Apr 2004, John Fox wrote:
I may misunderstand your question, so perhaps this answer isn't what you're
looking for:
A short time ago, Henrik Bengtsson posted a function to r-help (search for
plotSymbols in the list archive) that displays all available symbols. I
think that it's
At 05:16 AM 4/2/2004, Carlisle Thacker watched in amazement as electrons
turned into magical things called words:
Marc,
It is very difficult for the eye to distinguish even 25 symbols or 25
colors on the same plot. I find that my brain tends to saturate at 5 of
each, and using 5 symbols each
So what is your question? With R-1.8.1 on WinXPPro, I get:
y - scan(clipboard)
Read 48 items
f - factor(rep(1:2, 24))
fit - aov(y ~ f)
print(unclass(sum.fit)[[1]], dig=15)
Df Sum Sq Mean Sq F value Pr(F)
f1 3.63834187585414e-09
Marc R. Feldesman feldesmanm at pdx.edu writes:
Is there any effective way to get distinct geometric plotting symbols and
colors for plots involving more than 25 groups?
Not sure if geometric shapes is essential but if not you could plot them
with different letters (or LETTERS). The
Sorry, I apparently left out the line
sum.fit - summary(fit)
Andy
From: Liaw, Andy
So what is your question? With R-1.8.1 on WinXPPro, I get:
y - scan(clipboard)
Read 48 items
f - factor(rep(1:2, 24))
fit - aov(y ~ f)
print(unclass(sum.fit)[[1]], dig=15)
Df
A rainbow has a continuous distribution of wavelengths. These are
not at all the same thing as colors!
Quoting Prof Brian Ripley [EMAIL PROTECTED]:
You would do well to get people to recognize more than a dozen spot
colours, too -- we are tuned to recognizing quite large blobs of colour.
Liaw, Andy wrote:
So what is your question? With R-1.8.1 on WinXPPro, I get:
Obviously I was asking about how one would go about running one. I
appreciate your answer. Thanks.
--
chris ciotti ([EMAIL PROTECTED])
PGP ID: 0xE94BB3B7
__
[EMAIL
R stores its arrays in reverse odometer order where the leftmost
array index varies fastest. The result of:
matrix(1:4,2)
shows this.
Your pack and unpack are in odometer
order with the rightmost index varying fastest. If you were
to move to R's convention then you could move back and
forth
Hello!!
I need some help! I tried everything, but nothing worked!
I have a vector c with dates in it, in the format 2004-04-01 and i need to
convert it in the form
04/01/2004 or 01/04/2004 !
How can i do that??
Thanks
Bye Martina
--
__
[EMAIL
Hi, i can convert excell to list in R with package RODBC ()but i don't
understand 2 mistake
1) Don't read the last row of the table excell
2) Don' t take the hours
My excell file call prueba4.xls and have the following rows:
where prueba4.xls was make in excell (office xp) and have one
On Fri, 2 Apr 2004, Martina Azaroglu wrote:
I need some help! I tried everything, but nothing worked!
I am sure you did not try _everything_, but it would have been helpful to
give some idea of what you were trying.
I have a vector c with dates in it, in the format 2004-04-01 and i need to
On Fri, 2 Apr 2004 [EMAIL PROTECTED] wrote:
And the second question is the command odbcFechRows() don't work,
That is not a question, that is an assertion! There is no such function,
and I presume you mean odbcFetchRows().
for example if write
tbl-odbcFetchRows(canal) in place of sqlFetch
Convert your character dates to one of R's date classes
(POSIXlt, chron or the new Date class in 1.9.0) and then
format the date. The m/d/y format you want is actually the
default format in chron:
z - c( 2004-01-20, 2004-02-22 ) # test vector
# load chron and set default year to 4 digits
I have shapefiles for the state climatic divisions for the United States
and read.shape brings them in wonderfully. Now I wish to run through a
list of several thousand observation sites to find out in which division
each is located. I figure that I can compute the winding number for each
site
On Fri, 2 Apr 2004, Clint Bowman wrote:
I have shapefiles for the state climatic divisions for the United States
and read.shape brings them in wonderfully. Now I wish to run through a
list of several thousand observation sites to find out in which division
each is located. I figure that I
Roger,
Thanks for your reference. Since I can get the polygon coordinates (and
have the coordinates of my sites, I can cobble together a function that
will do the trick.
Again, thanks,
Clint
On Fri, 2 Apr 2004, Roger Bivand wrote:
On Fri, 2 Apr 2004, Clint Bowman wrote:
I have
Members of this list may be interested in
http://www.vividsolutions.com/jts/jtshome.htm. There is a C++ port asw
well by the postgis folks.
Tim
On Fri, 2004-04-02 at 12:29, Clint Bowman wrote:
Roger,
Thanks for your reference. Since I can get the polygon coordinates (and
have the
Where Can I Have More information about Formulae in R (Nesti8ng,
Crossing Factors, Etc). I´ve tried to obtain this information in R Help
and that source were a kind of inconclusive for my Doubts.
[]s
Leonard Assis
Estatístico
CONFE 7439
__
[EMAIL
Then you need to look elsewhere, like the reference in ?formula and the
books in the R FAQ.
I am his co-author, but I still think Bill Venables' account in MASS
(Venables Ripley, 1994, 7, 9, 2002) is unequalled.
On Fri, 2 Apr 2004, Leonard Assis wrote:
Where Can I Have More information about
I believe the white book has a chapter on the topic.
HTH,
Andy
From: Leonard Assis
Where Can I Have More information about Formulae in R (Nesti8ng,
Crossing Factors, Etc). I´ve tried to obtain this information
in R Help
and that source were a kind of inconclusive for my Doubts.
[]s
Hi,
Could you give me your opinions on this idea, please ?
Predicting Y from X is similar* to finding the Expectation{Y|X}.
The joint probability Pr(X,Y) 'contains' this expectation.
If this joint probability follows a 'simple' law, it is easier to
estimate it directly than using regression
I want a list of the number of times some factor levels appear together,
similar to the following SQL statement:
SELECT A, B, COUNT(C) FROM TBL GROUP BY A, B
How do I do that with a data.frame in R?
Thanks,
Jesper Frickmann
Statistician, Quality Control
Novozymes North America Inc.
77 Perry
I believe this should be an easy thing to do...
I have a function I repeatably call which takes input parameters and
outputs columns to various data frames. I also wish to keep a summary of
certain values as I call the function. I though keeping the values in a
vector then appending the vector by
I need to do a simple cumulative sum by group and add the result to the
data. I found an earlier thread in the help files with a few suggestions.
Somewhat, one of the suggestions does not work with my data, and I don't
really understand why ?
The error am getting using the my data below is...
You might want to look at ?table and ?xtab. See also ?tapply for use of
general functions (other than just COUNT) with GROUP BY.
HTH,
Vadim
-Original Message-
From: JFRI (Jesper Frickmann) [mailto:[EMAIL PROTECTED]
Sent: Friday, April 02, 2004 1:06 PM
To: [EMAIL PROTECTED]
Subject:
If you first look at the output of the inner core of your expression:
split(d.1, d.1$v1)
you might get a hint. Then make this change.
cbind(d.1, v4=c(lapply(split(d.1, d.1$v1), function(x) cumsum(x$v3)),
recursive=T))
There is no variable named y in your data frame, so x$y is the
wrong
Marko,
Looks fine to me. Why do you think the syntax is incorrect? Works for me in
1.8 on Windows.
a - 1
a
[1] 1
a - c(a,2)
a
[1] 1 2
a - c(a,2,4,5)
a
[1] 1 2 2 4 5
steve
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Mark O. Kimball
Sent:
Assuming you have a data.frame, dat, with coluns A and B, I think what you
want would be:
table(paste(dat$A,dat$B))
JFRI (Jesper Frickmann) wrote:
I want a list of the number of times some factor levels appear together,
similar to the following SQL statement:
SELECT A, B, COUNT(C) FROM TBL
Marko,
Looks fine to me. Why do you think the syntax is incorrect? Works for me in
1.8 on Windows.
I have a function I repeatably call which takes input parameters and
outputs columns to various data frames.
Looks fine to me although it is not completely clear that your request for
column and
On Friday 02 April 2004 04:52 pm, you wrote:
I believe this should be an easy thing to do...
I have a function I repeatably call which takes input parameters and
outputs columns to various data frames. I also wish to keep a summary
of
certain values as I call the function. I though keeping
From: =?iso-8859-1?q?J=20Swinton?= [EMAIL PROTECTED]
Has anyone achieved a current successful build of R
1.8 or 1.9 or earlier for HPUX Tru 64 OSF 5.1?
There is no binary version of R for HPUX Tru 64 OSF
5.1. The R admin manual mentions that the native make
fails on Alpha/OSF (aka Tru
Hi
this is my first time and i have no clue how this work
do you replay to my e-mail
or do u post the asnswer some other place
please let me know
the other question do u answer questions or give help
for questions written in SAS
for example i need the help to simulate from Bivariate
normal??
so
Hello,
I have spent considerable time trying to figure out that which I am about to describe.
This included searching Help, consulting my various R books, and trail and (always)
error. I have been assuming I would need to use a loop (looping over columns) but
perhaps and apply function
Amal -
Do:
library(MASS)
help(mvrnorm)
for documentation on the function you should use.
I found this, even without knowing the function name exactly,
by doing help.search(mvnorm). Please read the posting
guide and the FAQ.
- tom blackwell - u michigan medical school - ann arbor -
52 matches
Mail list logo