Dear all
is there any linear programming routine available for R?
if not, can you suggest any alternatives? not need to be very powerful, I
get only a samll problem to resolve.
many thanks
yong
__
[EMAIL PROTECTED] mailing list
The first menu item on the `File' menu is called `Source R code' and calls
source(). Did you look through the menus? If not, it would be worth
familiarizing yourself with them.
On Mon, 14 Jun 2004, Evgueni Parilov wrote:
Thanks!
That was exactly what I wanted.
Evgueni
Ko-Kang Kevin
Enrique Bengoechea wrote:
Hi,
Could someone please point to me which function clears the R console under Windows? (exactly what
the R Gui Windows menu Edit Clear Console does).
Seems simple but I haven't been succesful with the help system (help.search for console, clear
console, screen...) nor
Hi there;
I am a neophyte to R though I have been messing around with programming in other
languages and environments for some years (my dog's name is punchcard to give you some
idea of how many years). I have been trying to make a package and install it, to no
avail. The functions I have
On Tue, 15 Jun 2004, Yong Wang wrote:
is there any linear programming routine available for R?
Yes. help.search shows
print.simplex(boot) Print Solution to Linear Programming Problem
simplex(boot) Simplex Method for Linear Programming Problems
simplex.object(boot)Linear
On Mon, 14 Jun 2004, Ravi Varadhan wrote:
Is it possible to perform computations in quadruple precision (more
generally, with more digits in the floating-point arithmetic than that
allowed by double precision) in R?
You are making the assumption that R uses double precision. Whereas R
uses
Scott == Scott Norton [EMAIL PROTECTED]
on Tue, 15 Jun 2004 00:09:03 -0400 writes:
Scott This probably has a super easy answer...but I claim
Scott newbie status! (I did search help lists but this
Scott question is hard to isolate keyword-wise)
Scott Basically, I'm trying to
Dear all,
I have 122 vectors named from L1 to L122. Now I hope to
take log to each of the series, say
L1 - log(L1)
...
L122-log(L122)
Can anyone show me a iterative way to make the job simple.
I mean the way something like
for(i in 1:122){
...
}
or other similar methods.
Many thanks.
Jin
Hi!
If I have similar problem.
I first run R CMD check. Then I change to the directory mypackage.Rcheck where R CMD
check generated a file mypackage-manual.tex.
Then I try to compile the file using latex by myself. In addition I am using a latex
error aware editor (texniccenter on windows).
On Tue, 15 Jun 2004, Martin Maechler wrote:
If boot was programmed as we do other typical ``classes''.
The print method would print much less,
the summary method would compute all the interested quantities
and return an object of class summary.boot and there would be
a simple
I don't see why this needs to be iterative: the computations could be done
in parallel. But via a for loop you could use
for(i in 1:122) {
nm - paste(L, i, sep=)
assign(nm, log(get(nm))
}
You would do better, I think, to have 122 similar data items in a list.
datalist - lapply(1:122,
Hi folks,
Suppose I have a series of cases each with categorical
factors A, B.
What is the best way to glue A and B together into a single
factor? For example, given
A0 B1 ...
A1 B1 ...
A0 B2 ...
A1 B0 ...
A0 B0 ...
A1 B2 ...
then I'd like to end up with a single factor with levels
A0B0,
Dear Jin,
you could try something like,
L1 - runif(10)
L2 - runif(20)
L3 - runif(30)
L4 - runif(40)
res - vector(mode=list, length=4)
ss - paste(log(L, 1:4, ), sep=)
for(i in 1:4) res[[i]] - eval(parse(text=ss[i]))
res
I hope this helps.
Best,
Dimitris
Dimitris Rizopoulos
Doctoral
Prof Brian Ripley [EMAIL PROTECTED] writes:
On Tue, 15 Jun 2004, Martin Maechler wrote:
If boot was programmed as we do other typical ``classes''.
The print method would print much less,
the summary method would compute all the interested quantities
and return an object of class
(Ted Harding) [EMAIL PROTECTED] writes:
...
then I'd like to end up with a single factor with levels
A0B0, A0B1, A0B2, A1B0, A1B1, A1B2
according to all the combinations which actually occur in the
data (e.g. if (A1,B2) did not occur in the data, then A1B2 would
not be a level in the
Hi,
I am willing to load multiple C files to R, which
are inter-dependent (functions used in one may be
defined in other). What I was trying is to first
compile all of them separately (using R CMD SHLIB ...)
and then load them one by one (using dyn.load(...)
), but it doesnt work (which seems
On Mon, 14 Jun 2004, A Friend wrote:
Does anyone have any experience with SJava especially on Windows?
Yes, people do have experience, as searching the list archives will show
you. The experiences reported tend not to be good ones (selection bias,
of course).
The problem tends to be that
Dear Ted,
you could try something like,
dat - data.frame(V1=sample(LETTERS[1:3], 10, rep=T),
V2=sample(letters[1:3], 10, rep=T))
mat - apply(as.matrix(dat), 1, function(x) paste(x, collapse=))
dat$glue - factor(mat)
dat
I hope this helps.
Best,
Dimitris
Dimitris Rizopoulos
Doctoral
On Tue, 15 Jun 2004, Utsav Boobna wrote:
Hi,
I am willing to load multiple C files to R, which
are inter-dependent (functions used in one may be
defined in other). What I was trying is to first
compile all of them separately (using R CMD SHLIB ...)
and then load them one by one (using
On 15-Jun-04 Peter Dalgaard wrote:
(Ted Harding) [EMAIL PROTECTED] writes:
...
then I'd like to end up with a single factor with levels
A0B0, A0B1, A0B2, A1B0, A1B1, A1B2
according to all the combinations which actually occur in the
data (e.g. if (A1,B2) did not occur in the data,
On Tue, Jun 15, 2004 at 01:43:21AM -0700, Utsav Boobna wrote:
Hi,
I am willing to load multiple C files to R, which
are inter-dependent (functions used in one may be
defined in other). What I was trying is to first
compile all of them separately (using R CMD SHLIB ...)
and then load them
Hi, Ted:
Probably not the most elegant way, but not too tedious either.
a - as.factor(c(A0,A1,A2))
b - as.factor(c(B0,B1,B2))
as.vector(outer(a,b,FUN=paste,sep=))
[1] A0B0 A1B0 A2B0 A0B1 A1B1 A2B1 A0B2 A1B2 A2B2
Best, ft.
--
Fernando TUSELLe-mail:
Departamento de
Apologies, disregard my former message! I overlloked the requirement
that only combinations
present in the data must be taken. ft.
--
Fernando TUSELLe-mail:
Departamento de Econometría y Estadística [EMAIL PROTECTED]
Facultad de CC.EE. y Empresariales
Pamela Hall wrote:
Hi there;
I am a neophyte to R though I have been messing around with programming in other languages and environments for some years (my dog's name is punchcard to give you some idea of how many years). I have been trying to make a package and install it, to no avail. The
Does anyone know how to read a SYSTAT .syd file on Linux?
(Splus 6 does it, but it is easier to find a Windows box
with Systat than to download their demo. I'm wondering
if there is a better way than either of these options.)
Jon
--
Jonathan Baron, Professor of Psychology, University of
Dear all
Suppose I have teeth like data similar like
x - 1:200
y - 0.03*x[1:100]+rnorm(100, mean=.001, sd=.03)
z - 3-rep(seq(1,100,10),each=10)*.03+rnorm(100,mean=.001, sd=.03)
plot(x,c(y,z))
and I want to have a gradient estimations for some values from increasing part of
data
like
y.agg -
Jonathan Baron [EMAIL PROTECTED] writes:
Does anyone know how to read a SYSTAT .syd file on Linux?
(Splus 6 does it, but it is easier to find a Windows box
with Systat than to download their demo. I'm wondering
if there is a better way than either of these options.)
Jon
--
Jonathan
Dear users
I have a problem with the dr function: dimension reduction.
I give you my example, and i'll be pleased to read your comments.
#let be X a matrix 50*100:
library(dr);
X- matrix(rnorm(50*100,5,1),50,100);
#and let be Y a vector response:
Y- sample(0:1,50,replace=T);
#I choose (for
Dear Petr,
Probably I don't understand exactly what you are looking for.
However your plot(x,c(y,z)) suggests a broken-line model for the response
c(y,x) versus the variables x. Therefore you could estimate a segmented
model to obtain (different) slope (and breakpoint) estimates. See the
package
Can anyone explain to me why the AIC values are so different when
using glm.nb and glm with a negative.binomial family, from the MASS
library? I'm using R 1.8.1 with Mac 0S 10.3.4.
library(MASS)
dfr - data.frame(c=rnbinom(100,size=2,mu=rep(c(10,20,100,1000),rep(25,4))),
+
Having not tried this, it is dangerous to speculate, but it appears to
me that there
would be no problem passing rq arguments (crucially, only tau, the
specification
of the quantile of interest) to fit.mult.impute, since the call to the
fitter procedure
includes a ... argument. The real
Hello everyone, is there a package/packages for factor analysis,
particularly PCA?
thanks,
Katja
Katja Löytynoja
Taitoniekantie 9 A 218
40 740 Jyväskylä
Finland
tel.+35814 608058
cell.+35850 336 0174
[EMAIL PROTECTED]
__
[EMAIL PROTECTED]
Katja Loytynoja [EMAIL PROTECTED] writes:
Hello everyone, is there a package/packages for factor analysis,
particularly PCA?
help.search(factor analysis)
help.search(principal components)
(Whether PCA qualifies as factor analysis is debatable, though...)
--
O__ Peter Dalgaard
R comes with support for factor analysis and PCA (*not* the same thing)
in package stats which is normally loaded.
Try
help.search(factor analysis)
help.search(principal components)
On Tue, 15 Jun 2004, Katja Loytynoja wrote:
Hello everyone, is there a package/packages for factor analysis,
Have you checked really carefully for Unbalanced braces? Many editors Please do so
an look whether opening and closing of barces happens in the correct order. Also,
watch out for escaped braces and stuff like that. Some appropriate editor might help,
but don't rely on it (since it's not the
Maybe you can show us one of the Rd files that doesn't work, so some one on
this list can see what the problem might be?
Cheers,
Andy
From: Pamela Hall
Have you checked really carefully for Unbalanced braces?
Many editors Please do so an look whether opening and closing
of barces happens
On Tue, 15 Jun 2004, Benjamin Esterni wrote:
I have a problem with the dr function: dimension reduction.
It seems that you are using it inappropriately.
I give you my example, and i'll be pleased to read your comments.
#let be X a matrix 50*100:
library(dr);
You should not be
Pamela == Pamela Hall [EMAIL PROTECTED]
on Tue, 15 Jun 2004 08:41:27 -0400 writes:
...
Pamela I must be doing something else wrong. Its not the braces.
Pamela This is what I have done that may not be right?
Pamela 1. I put all of the Rd files into 1 file with a
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
On 15 Jun 2004 at 13:52, Vito Muggeo wrote:
Dear Petr,
Probably I don't understand exactly what you are looking for.
However your plot(x,c(y,z)) suggests a broken-line model for the
response c(y,x) versus the variables x. Therefore you could estimate
a segmented model to obtain
On Tue, 15 Jun 2004 13:11:45 +0200 Petr Pikal wrote:
Dear all
Suppose I have teeth like data similar like
x - 1:200
y - 0.03*x[1:100]+rnorm(100, mean=.001, sd=.03)
z - 3-rep(seq(1,100,10),each=10)*.03+rnorm(100,mean=.001, sd=.03)
plot(x,c(y,z))
and I want to have a gradient
On Mon, 14 Jun 2004 22:43:27 +0100 (BST), (Ted Harding)
[EMAIL PROTECTED] wrote :
On 14-Jun-04 Duncan Murdoch wrote:
Please try 1.9.1 beta. This should be fixed now...
...
Does this mean that barplot()
got un-fixed between 1.8.0 and 1.9.0?
Yes, and this is why we *really strongly* encourage
Hi list,
I'm attempting to re-create a Repeated Measures Compositional Analysis
as
described in the work by Aebischer et. al. (Ecology. 1993. 74(5):
1313-1325).
In this paper they describe transitions of data into a log ratio
difference matrix, from which they obtain two matrices using a
You have the code, so take a look for youself. There's an error in the
aic formula in negative.binomial (and neg.bin)
On Tue, 15 Jun 2004, Ken Knoblauch wrote:
Can anyone explain to me why the AIC values are so different when
using glm.nb and glm with a negative.binomial family, from the
Hi!
I have graph with 1 or more (up 1000k) nodes and only few edges (can be weighted
or not). It is a really sparse matrix.
I can generate this graph using R. Which structure I can, should, use to store the
edges? (mSparse?)
Having this datastructures which clustering algorithm I can use
lsoda doesn't pass along the names attribute of the state vector, y. If
you want to use the names of the state vector in your code, you need to
reassign it inside your ode function. I should either fix this in the
code for lsoda, or at least document it! For example, I can run the
following
Yep, shame on me. But on the other hand, this is a good example of a
vague menu
item. I think that instead of Source R code..., it should be Load
source R code...
unless it can do more than just loading a source from R file.
Evgueni
Prof Brian Ripley wrote:
The first menu item on the `File'
Dear
Thank you for contacting Texas Instruments Incorporated, an educational
technology company. Your message has been routed to a TI-Cares Agent and
you will receive a solution shortly. For 24/7 customer self-service,
please visit our knowledgebase at education.ti.com/support.
Have a
roger koenker wrote:
Having not tried this, it is dangerous to speculate, but it appears to
me that there
would be no problem passing rq arguments (crucially, only tau, the
specification
of the quantile of interest) to fit.mult.impute, since the call to the
fitter procedure
includes a ...
Evgueni Parilov wrote:
Yep, shame on me. But on the other hand, this is a good example of a
vague menu
item. I think that instead of Source R code..., it should be Load
source R code...
unless it can do more than just loading a source from R file.
No. source() does not only load, but also
On Tue, 15 Jun 2004 12:50:04 -0400, Evgueni Parilov
[EMAIL PROTECTED] wrote :
Yep, shame on me. But on the other hand, this is a good example of a
vague menu
item. I think that instead of Source R code..., it should be Load
source R code...
unless it can do more than just loading a source from
Thanks Uwe!
Now I understand how it works and why it was called Source R code
rather than
Load source R code. I did not really intend to confuse anybody. I
just wanted
to have full information on a subject.
I agree that most of the menu items are self-explained. But some of them
may not!
Jens Schumacher wrote:
knussear wrote:
Hi list,
I'm attempting to re-create a Repeated Measures Compositional Analysis as
described in the work by Aebischer et. al. (Ecology. 1993. 74(5):
1313-1325).
In this paper they describe transitions of data into a log ratio
difference matrix, from which
Hi R-users!
I'd like to know how to avoid automatic coercion of numeric variables
passed to array() to character, when there is one categorical variable
between others. I have tried I() function, but it is just functional
in data.frame.
Best Regards!
Eduardo Dutra de Armas
YEAH YEAH YEAH!! I got a bunch of functions to install with their help pages and even
help.search works!!! The INDEX file is just fine. YEAH YEAH YEAH!
This is what I have learned:
1. DO NOT concatenate the *.Rd files. This happens in the installation process. This
is probably completely
The package I am writing is for the Center for Tropical Forest Science, CTFS. This
center is a collaboration of 15+ institutions world wide that are investigating
properties of tropical forest dynamics, species diversity, species distributions. The
investigation is composed of the same
On Jun 15, 2004, at 10:11 AM, Douglas Bates wrote:
Jens Schumacher wrote:
knussear wrote:
Hi list,
I'm attempting to re-create a Repeated Measures Compositional
Analysis as
described in the work by Aebischer et. al. (Ecology. 1993. 74(5):
1313-1325).
In this paper they describe transitions of
From: Pamela Hall
YEAH YEAH YEAH!! I got a bunch of functions to install with
their help pages and even help.search works!!! The INDEX
file is just fine. YEAH YEAH YEAH!
This is what I have learned:
1. DO NOT concatenate the *.Rd files. This happens in the
installation process.
My $0.02...
From: Pamela Hall
Hi again;
I have placed many examples in the CTFS Rd pages because my
audience really really likes lots of examples in order to
understand the options available in a functions. At this
time I have set them all to \dontrun{}.
however, this isn't
Thank you! The next step in the conversion still fails and I can't seem to find any
examples in the archives. The result of the function 'as.POSIXct(strptime())' within
the 'sapply' comes back as numeric rather than POSIXct as expected:
ds - cbind(1:2, c(02/27/92 23:03:20, 02/27/92
I'm looking for an optimal approach to access Oracle databases via RWeb
applications. I'm new to R but familiar with programming functions and web
pages for the S+ Statserver. I'm now going through the motions of migrating
S+/Statserver applications to R/RWeb as a feasability exercise. I can
On Tue, 15 Jun 2004, Evgueni Parilov wrote:
Yep, shame on me. But on the other hand, this is a good example of a
vague menu
item. I think that instead of Source R code..., it should be Load
source R code...
unless it can do more than just loading a source from R file.
It is `source' not
Sounds like a permissions/ownership problem. odbcConnect returnign -1
just means that the ODBC device manager failed.
On Tue, 15 Jun 2004, Fowler, Mark wrote:
I'm looking for an optimal approach to access Oracle databases via RWeb
applications. I'm new to R but familiar with programming
I'm trying to install a package in windows. I have a package directory,
which contains all the essentials: .Rd in the Man directory, DESCRIPTION
file, etc. I copied it to my R\bin directory, where the Rcmd file is located.
Then, when I try to run Rcmd build mypackage from that directory in
On Tue, 15 Jun 2004, Yev Kirpichevsky wrote:
I'm trying to install a package in windows. I have a package directory,
which contains all the essentials: .Rd in the Man directory, DESCRIPTION
file, etc. I copied it to my R\bin directory, where the Rcmd file is located.
Then, when I try to run
On Tue, 15 Jun 2004, Pamela Hall wrote:
Hi again;
I have placed many examples in the CTFS Rd pages because my audience
really really likes lots of examples in order to understand the options
available in a functions. At this time I have set them all to
\dontrun{}.
however, this isn't
Thanks to everyone who responded to my 'Perl' problem! Indeed, it wasn't
in the Path.
I now have a different problem and would appreciate any insights:
after I run Rcmd build mypackage, I get the following error:
'sh' is not recognized as a command, program, or batch file
Thanks in advance,
Please follow the instructions in readme.packages very, very, very, very
carefully.
Andy
From: Yev Kirpichevsky
Thanks to everyone who responded to my 'Perl' problem!
Indeed, it wasn't
in the Path.
I now have a different problem and would appreciate any insights:
after I run Rcmd
Hi everyone,
I am trying to build a R library called nnNorm but I have some troubles
checking and installing it.
Here is the setup:
If I don't use a inst\doc directory(with the vignettes files) in my source,
the install step is working fine:
C:\testRCMD INSTALL nnNorm
-- Making
On Tue, 15 Jun 2004 18:17:36 -0400, [EMAIL PROTECTED] wrote:
hello, everybody,
I met a problem that I want to generate a random uniform number by using
function with c interface.
What I found is only unif_rand(), but its range is [0,1]. How can I create a
uniform random variable by using a
I am trying to perform a model 3 ANOVA for a 2 factor (say factor A and
factor B) anova in which factor A is fixed and factor B is random.
Therefore, the error term for the test of factor A should be the A:B
interaction term and the error terms for B and A:B should be the model
residual
(BNote: JP stands for Japanese.
(B
(BALERT!!
(BThis e-mail contained one or more virus-infected files and have been rejected.
(B(JP:$B%3%s%T%e!<%?%&%#%k%9$rH/8+$7$^$7$?$N$G!"%a!<%k$NAw?.$rCf;_$7$^$7$?!#(J)
(B
(BThe following attachments were infected:
(B(JP:[EMAIL
Hi R People:
I have a monthly time series x.ts which runs from 1/1995 through 12/2003.
x.ts - ts(x,start=1995,freq=12)
str(x.ts)
Time-Series [1:108] from 1995 to 2004: -1.638 -0.236 0.830 -0.548 0.363
...
My question: is there a way to print the observations from 1/1999 to 6/1999,
please?
?window
Date: Tue, 15 Jun 2004 18:43:04 -0500
From: Laura Holt [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Subject: [R] time series object
Hi R People:
I have a monthly time series x.ts which runs from 1/1995 through 12/2003.
x.ts - ts(x,start=1995,freq=12)
str(x.ts)
Time-Series
Hi Murray
I am trying to perform a model 3 ANOVA for a 2 factor (say factor A and
factor B) anova in which factor A is fixed and factor B is random.
...
In addition, I have tried using lme to perform this function,
but again without much success.
What did not work? And, did you read
Hi Chris,
I am trying to fit the following linear model to logged per capita
fecundity data (ie number of babies per female) for a mouse:
RsNRlS - glm(formula = ln.fecundity ~ summer.rainfall + N +
lagged.rainfall + season,
)
I am using this relationship in a simulation model, and
76 matches
Mail list logo