Greetings everyone;
I am not an R user; however, I am attempting to install R 2.7.0 on a
multi-user server for use by students and faculty in the Statistics
department here. I have an old version of R installed on the same
server, and it works fine; however, I am running into a problem
in
I am looking for a way to automatically detect outliers greater then X
standard devations then (here is the hard part) APPEND those dummy
variables to an existing array of regression variables (called cReg).
Any suggestions on how to go about this?
(source) sample data set that I am modelling:
Date
Could someone please point me to where information about
the two entries %\VignetteDepends{} and %\VignettePackage{}
are documented?
I'm getting an error with %\VignetteDepends{} with my
package name in braces. Looking at the vignettes on my
system, some include the package name and some don't.
Sh
On May 27, 2008, at 12:45 PM, ctu wrote:
Hi,
I don't have any problem to install the R 2.7.0 in XP and windows
Vista.
I download R 2.7.0 from http://cran.cnr.berkeley.edu/
You could try it again.
If you still have a problem to install it. You might want to consult
with your department PC
On Tue, 27 May 2008, Stan Horwitz wrote:
On May 27, 2008, at 12:45 PM, ctu wrote:
Hi,
I don't have any problem to install the R 2.7.0 in XP and windows Vista.
I download R 2.7.0 from http://cran.cnr.berkeley.edu/
You could try it again.
If you still have a problem to install it. You might wan
My Linux is pretty rusty, but I think the error message means what is says:
you do not have readline header (*.h) and library files. They usually
reside in readline-devel package if I remember correctly.
Hope this helps,
Andy
__
Andy Jaworski
518-1-01
Process Lab
On 5/27/08, Charilaos Skiadas <[EMAIL PROTECTED]> wrote:
>
> On May 27, 2008, at 3:10 AM, Hofert Marius wrote:
>
>
> > Hi,
> >
> > Consider the following cloud()-plot consisting of two lines:
> >
> > library(lattice)
> >
> mydata=data.frame(x=rep(1:4,2),y=rep(1:2,each=4),z=c(1,2,3,4,4,3,2,1))
> >
Hello there!!!
Sorry to bother you all with such question and difficulties that I have been
facing on.
Recently I have been searching for packages to run GMM estimatives with R.
I have been searching for such packages for a while, but since I am a new
user of R system,
my quest so far was unsuce
I am trying to use R lm() with quantitative and qualitative predictors, but am
getting different results than those that I get in SAS.
In the R ANOVA table documentation I see that "Type-II tests corresponds to the
tests produced by SAS for analysis-of-variance models, where all of the
predictors
This definitely does the trick.
I knew there was an easier way!
Petr Pikal wrote:
>
> [EMAIL PROTECTED] napsal dne 27.05.2008 09:31:16:
>
>>
>> Hi Petr
>>
>> My mistake in omitting "replace=T" in the first part.
>> Unfortunately I oversimplified my problem as I'm not actually dealing
> with
were you meaning Gaussian mixture model? If so, try the mclust package.
But you might want to try to learn how to search using R website
resourse, try
RSiteSearch('GMM')
X
Correia 写道:
> Hello there!!!
>
> Sorry to bother you all with such question and difficulties that I have been
> facing on.
K. Elo wrote:
Hi,
Monica Pisica wrote:
> - There is no perfect “beginner” book.
How about
- Crawley, Michael (2007). The R book, Wiley & Sons.
- Maindonald, John & John Braun (2007): Data Analysis and Graphics Using
R (2nd edition), Cambridge University Press.
As a political scientist (with
I have a matrix of frequency counts from 0-160.
x<-as.matrix(c(0,1,0,0,1,0,0,0,1,0,0,0,0,1))
I would like to apply a function creating a new column (x[,2])containing
values equal to:
a) log(x[m,1]) if x[m,1] > 0; and
b) for all x[m,1]= 0, log(next x[m,1] > 0 / count of preceding zero values
+1)
Thanks to everyone who responded to my plea for help, using ...
./configure --with-readline=no
worked fine, and finally reading doc/html/R-admin.htm helped.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do re
On Tue, 27-May-2008 at 06:17AM -0700, DavidM.UK wrote:
|>
|> I don't think it's an issue with Linux functionality at all, it's a
|> case of having your linux machine configured correctly -- which
|> isn't a straight forward issue.
Since configuration has a direct bearing on functionality, we're
Leo, do you mean GMM as in Generalized Method of Moment? If so i've posted a
similar question here before and the answer appears to be none at this time.
Horace
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Leo Correia
Sent: Tuesday, May 27, 2008 11:53
On Tue, 27 May 2008, T.D.Rudolph wrote:
I have a matrix of frequency counts from 0-160.
x<-as.matrix(c(0,1,0,0,1,0,0,0,1,0,0,0,0,1))
I would like to apply a function creating a new column (x[,2])containing
values equal to:
a) log(x[m,1]) if x[m,1] > 0; and
b) for all x[m,1]= 0, log(next x[m,1
Dear all,
I stumbled over a problem recently when trying to use srt with text() on a
windows device.
What I intended to do was to plot a simple regression line, and to rotate
a piece of text such that the text has the same angle as the regression
line.
However, the text is always plotted in a sl
I think this comment for ?par, meant for both crt and srt, applies:
crt
A numerical value specifying (in degrees) how single characters
should be rotated. It is unwise to expect values other than multiples
of 90 to work. Compare with srt which does string rotation.
So I would say that even
Hi,
I'm trying to get the beginning of my zoo object's index through start
(d) but sometimes get an error:
Error in start(d):
Unused argument (s) (c(NA, 8.792, 8.725, ...)
The data(d) is retrieved from Bloomberg through RBloomberg. I am
pulling 6 tickers (which means I'm getting 6 c
Note that the scale of x-axis and y-axis is different in your plot. One
simple way to avoid this is to keep the data unit in the x direction is
equal that in the y direction, by setting asp=1 in calling plot function.
X
Dr. Christoph Scherber 写道:
Dear all,
I stumbled over a problem recently
Thanks to all for the postings so far!
I found that setting asp=0.5 and then dividing the slope by 2 seems to do
the trick:
##
x=1:10
y=x*2-rnorm(1:10)
plot(x,y,pch=16,asp=0.5)
abline(lm(y~x))
yval=predict(lm(y~x),list(x=rep(2,length(x[1]
slope=as.numeric(lm(y~x)[[1]][2])
text(2,yval,"Regre
On Tue, 27 May 2008, Xiaohui Chen wrote:
Note that the scale of x-axis and y-axis is different in your plot. One
simple way to avoid this is to keep the data unit in the x direction is equal
that in the y direction, by setting asp=1 in calling plot function.
Or when that is not possible or n
My best guess at your data based on what you posted works fine:
> d <- zoo(c(NA, 8.792, 8.725))
> start(d)
[1] 1
> packageDescription("zoo")$Version
[1] "1.5-3"
> R.version.string # Vista
[1] "R version 2.7.0 Patched (2008-05-11 r45672)"
Try sending me the output of dput(d) along with the versio
Hi,
A friend just introduced me to R today and I think it is really nice after
browsing its web site. I'm eager to cook up some web-based interface to use
R at work for more platform-independent access. I plan to implement it in
a Python-based framework such as TurboGears or Django on either Fre
Check out http://www.rpad.org
On Tue, May 27, 2008 at 6:26 PM, Hseu-Ming Chen <[EMAIL PROTECTED]> wrote:
> Hi,
> A friend just introduced me to R today and I think it is really nice after
> browsing its web site. I'm eager to cook up some web-based interface to use
> R at work for more platform-i
Hi Werner,
One possibility is to use the xlsReadWritePro package
which allows to do this quite easily.
Regards,
Moshe.
--- Werner Wernersen <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I have an Excel (2000 / 2003) workbook with multiple
> sheets on which I would like to perform several
> statistical
Hello,
I am looking for a function that computes a consensus DNA sequence by
plurality. I found "consensus" in bio3d which almost does what I need.
However, it would be important for me to include ambiguities and not to
omit every position that is less than the threshold set in "consensus".
I
In fact x[4,2] should = log(x[5,1]/2]
whereas x[3,2] = log(x[5,1/3])
i.e. The denominator in the log function equals the number of rows between
m==0 and m>0 (inclusive, hence the "+1")
Hope this helps!...
Charles C. Berry wrote:
>
> On Tue, 27 May 2008, T.D.Rudolph wrote:
>
>>
>>
>> I ha
Dear R-Gurus,
I have ended up with a calculation problem where I need to use brobs.
I have to work my way through a vector with a for loop to act on each
element in a calculation (refering to the previous
value in the new vector of results -- so as far as I know I can't use
"apply") -- this produ
Hi all,
Sorry about this, my brain finally ticked over.
I just
need to convert the list of brobs to a numeric vector (which I can do,as
the answer to
the calculation need not be a brob -- just the terms along the way).
so
nums<-sapply(t,as.numeric)
Does the trick in my case, I think.
Ma
Hi,
I have devised some applications of R as a scripting shell language and
started using it flawlessly a few months ago.
Never minded publishing it because it looked so obvious. However, a search
in the archives revealed there is no recent easy-to-use method proposed.
The latest reference I f
Does this do what you want:
> x<-c(0,1,0,0,1,0,0,0,1,0,0,0,0,1)
> y <- rle(x)
> result <- lapply(seq_along(y$lengths), function(.indx){
+ if (y$values[.indx] == 0)
log(y$values[.indx+1]/seq(y$lengths[.indx]+1, by=-1,
length=y$lengths[.indx]))
+ else rep(log(y$values[.indx]), y$lengths[.ind
On Tue, 2008-05-27 at 13:38 -0700, Horace Tso wrote:
> Leo, do you mean GMM as in Generalized Method of Moment? If so i've posted a
> similar
> question here before and the answer appears to be none at this time.
>
Generalized Method of Moments (GMM) is an estimation method that can be
appli
On Tue, May 27, 2008 at 09:28:09PM -0300, Alexandre Aguiar wrote:
> I have devised some applications of R as a scripting shell language and
> started using it flawlessly a few months ago.
>
> Never minded publishing it because it looked so obvious. However, a search
> in the archives revealed th
Well done Jim.
It always ends up a little more complicated than I originally hoped, but I
certainly couldn't improve on it!
jholtman wrote:
>
> Does this do what you want:
>
>> x<-c(0,1,0,0,1,0,0,0,1,0,0,0,0,1)
>> y <- rle(x)
>> result <- lapply(seq_along(y$lengths), function(.indx){
> + i
In fact x[4,2] should = log(x[5,1]/2]
whereas x[3,2] = log(x[5,1/3])
i.e. The denominator in the log function equals the number of rows between
m==0 and m>0 (inclusive, hence the "+1")
Hope this helps!
Charles C. Berry wrote:
>
> On Tue, 27 May 2008, T.D.Rudolph wrote:
>
>>
>>
>> I have a ma
Dear R News,
I am inquiring whether the articles that are published in the 'R News'
have undergone a peer-review/refereeing process, for the purpose of
reporting the academics publications here at the University of Sydney to
the Australian Government.
Thank you for your assistance.
Rega
Dear Glenn,
As I responded to your direct inquiry to me, articles in R News are
peer-reviewed.
Regards,
John
On Wed, 28 May 2008 11:48:18 +1000
"Glenn Howell" <[EMAIL PROTECTED]> wrote:
> Dear R News,
>
>
>
> I am inquiring whether the articles that are published in the 'R
> News'
> have u
On 28/05/2008, at 1:48 PM, Glenn Howell wrote:
Dear R News,
I am inquiring whether the articles that are published in the 'R News'
have undergone a peer-review/refereeing process, for the purpose of
reporting the academics publications here at the University of
Sydney to
the Australian Go
Dear all,
I could not draw a forest plot for meta-analysis under ramdon models using the
rmeta package. The rmeta has a default function for MH (fixed-effect) model.
Has the rmeta package been updated for such a function? Or someone revised it
and kept a private code?
I would appreciate it
Dear all,
I want to create a list that contains 0,1,2,3, ..., 1 as its elements. I
used the following code, which apparently doesn't work very well.
a <- 0
for(i in 1:1) {
a <- list(a, i)
}
The result is not what I wanted. So how to create the bind lists recursively
so that the la
On May 27, 2008, at 10:43 PM, Daniel Yang wrote:
Dear all,
I want to create a list that contains 0,1,2,3, ..., 1 as its
elements. I used the following code, which apparently doesn't work
very well.
a <- 0
for(i in 1:1) {
a <- list(a, i)
}
The result is not what I wanted. So ho
On 28/05/2008, at 2:43 PM, Daniel Yang wrote:
Dear all,
I want to create a list that contains 0,1,2,3, ..., 1 as its
elements. I used the following code, which apparently doesn't work
very well.
a <- 0
for(i in 1:1) {
a <- list(a, i)
}
The result is not what I wanted. So how
This is somewhat subtle.
Rolf's solution (here corrected to...)
a <- list()
for(i in 0:1) a[[i+1]] <- i
is the best of the loop solutions (or at least the best I know of). The
apparently similar
a <- list(0)
for(i in 1:1) a <- c(a, list(i))
will take a lot longer, although the result
Unfortunately Reduce is just as slow as the slowest of
the solutions below:
Reduce(c, 0:1, list())
On Tue, May 27, 2008 at 11:18 PM, <[EMAIL PROTECTED]> wrote:
> This is somewhat subtle.
>
> Rolf's solution (here corrected to...)
>
> a <- list()
> for(i in 0:1) a[[i+1]] <- i
>
> is the
Dear all,
May anyone explain to me how I run a linear programming or Data Envelopment
Analysis (DEA models) into R?
Thanks a lot.
Best Regards.
Marcus Vinicius
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://s
In the following example:
x <- rnorm(1:100)
y <- seq(from=-2.5, to=3.5, by=0.5)
z <- as.matrix(table(cut(x,c(-Inf, y, +Inf
## I wish to transform the values in z
j <- log(z)
## Yet retain the row names
row.names(j)<-row.names(z)
Now, how can I go about creating a scatterplot with row.names(
Stan Horwitz wrote:
Thanks to everyone who responded to my plea for help, using ...
./configure --with-readline=no
worked fine, and finally reading doc/html/R-admin.htm helped.
Unless of course you actually wanted the arrow keys to work on the
command line. In which case, read R-admin again...
Does anyone have a program for comparing interaction means (similar to HSD.test
of agricolae) for a split or strip plot design? For those who are familiar with
SAS I am looking for something similar to the DIFF option of the LSMEANS
statement of PROC MIXED and the PDMIX800 macro. The diff option
But pre-allocation still helps
a <- vector("list", 10001)
for(i in 0:1) a[[i+1]] <- i
gives
user system elapsed
0.042 0.001 0.057
on my system, against
user system elapsed
0.743 0.000 1.907
for Bill's 'best' solution.
On Wed, 28 May 2008, [EMAIL PROTECTED] wrote:
Thi
On my system
help.search("linear program")
gives (slightly edited)
linp(limSolve) Linear Programming
solveLP(linprog)solve Linear Programming / Optimization
problems
lp.object(lpSolve) LP (linear programming) object
lpcdd(rcdd) linear
On Tue, 27 May 2008, T.D.Rudolph wrote:
In the following example:
x <- rnorm(1:100)
y <- seq(from=-2.5, to=3.5, by=0.5)
z <- as.matrix(table(cut(x,c(-Inf, y, +Inf
## I wish to transform the values in z
j <- log(z)
## Yet retain the row names
row.names(j)<-row.names(z)
Hmm. The rownames
Hi Marcus!
On Wednesday 28 May 2008 05:56, Marcus Vinicius wrote:
> Dear all,
> May anyone explain to me how I run a linear programming or Data
> Envelopment Analysis (DEA models) into R?
Package "DEA" (on CRAN):
http://cran.r-project.org/web/packages/DEA/index.html
Package "FEAR" (NOT on CRAN
54 matches
Mail list logo