methods.
It does, in detail, in sections 2.4 and 2.5.
Thanks you,
Shengzhe
On 8/26/05, Prof Brian Ripley [EMAIL PROTECTED] wrote:
On Fri, 26 Aug 2005, Shengzhe Wu wrote:
Thanks for your reply. Actually I called function as below.
p1 = predict(object, newdata, dimen=1)
p2 = predict(object
I suspect you have not print()-ed your graphics, see FAQ Q7.22.
It is then possible to include dev.off() within the function. E.g.
testit - function(fn = test.eps)
{
trellis.device(postscript, file=fn, onefile = FALSE, horizontal=FALSE)
print(stripplot(voice.part ~ jitter(height), data =
On Fri, 26 Aug 2005, Tarca, Adi wrote:
Hi everyone,
According to R reference manual, the nnet function uses the BFGS method
of optim to optimize the neural network parameters.
What the help page says is
...: arguments passed to or from other methods.
That means methods of nnet().
Try profiling. Doing this many times to get an overview, e.g. for sapply
with df=1:1000:
% self% total
self secondstotalsecondsname
98.26 6.78 98.26 6.78 FUN
0.58 0.04 0.58 0.04 unlist
0.29 0.02
Look at ?[.factor:
finaldataset$Species - finaldataset$Species[,drop=TRUE]
solves this.
On Fri, 26 Aug 2005, Martin Lam wrote:
Hi,
Since I've had no replies on my previous post about my
problem I am posting it again in the hope someone
notice it. The problem is that the
Are there NAs in the variable?
SYNTAX==Ditrans and SYNTAX!=Ditrans are not mutually exclusive.
On Fri, 26 Aug 2005, Stefan Th. Gries wrote:
Dear all
I have a problem with splitting up a data frame called ReVerb:
» str(ReVerb)
`data.frame': 92713 obs. of 16 variables:
$ CHILD: Factor
That is one of the S4 vs R differences. See the complements.
On Fri, 26 Aug 2005, Phineas Campbell wrote:
I'm working through the examples in Venables and Ripley in the 'New-style
Classes' chapter.
On a call to representation, in the lda example, it is unable to find the
class named.
Is
On Fri, 26 Aug 2005, Luwis Tapiwa Diya wrote:
I have been trying to write my own user defined function in Rpart.I
imitated the anova splitting rule which is given as an example.In the
work I am doing ,I am calculating the concentration index(ci) ,which
is in between -1 and +1.So my deviance
?cut
This is in `An Introduction to R', the manual which ships with R and basic
reading.
On Fri, 26 Aug 2005, David James wrote:
What is the quickest way to create many categorical variables
(factors) from continuous variables?
This is the approach that I have used:
# create sample data
On Fri, 26 Aug 2005, David James wrote:
I create a zooreg object that runs from Jan-1-2002 0:00 to Jun-1-2005
0:00...
regts.start = ISOdatetime(2002, 1, 1, hour=0, min=0, sec=0, tz=)
regts.end = ISOdatetime(2005, 6, 1, hour=0, min=0, sec=0, tz=)
regts.zoo - zooreg( NA, regts.start,
The survival package is a recommended package in R and contains survreg()
which uses the AFT definitions for Weibull survival. This is well
documented, and MASS (the book) has comparisons of PH and AFT
parametrization for a Weibull example.
I think you mean Frank Harrell's `Design' package.
On Sat, 27 Aug 2005, Benn Fine wrote:
I want to fit a logistic regression model under a specified
null hypothesis, say Ho:Beta_k=1
Is there a way to constrain a parameter in glm.fit ??
You should be calling glm(). Then you can use offset() in your formula.
(You can also use it as an
Looks like the package is not written for current R. native.enc was used
prior to R 2.1.0. Please take the advice of the posting guide and ask the
maintainer.
BTW, there is no package ttda on CRAN, and you have not told us where you
found it or what version you used. If you mean
What is an `mle object'? If it is something produced by mle() in stats4,
in what sense is `prediction' relevant?
On Mon, 29 Aug 2005 [EMAIL PROTECTED] wrote:
Hello.
I am looking for a reference to compute prediction intervals from a mle
object.
Thank you.
Tony.
[[alternative
On Mon, 29 Aug 2005, Ferran Carrascosa wrote:
Hi,
I have a matrix with 700.000 x 10.000 cells with floating point data.
I would like to work with the entire table but I have a lot of memory
problems. I have read the ?memory
I work with Win 2000 with R2.1.0
The only solution that I have
I was puzzled as to what the question actually was. If you set
options(contrasts=c(contr.sum, contr.poly))
the interaction contrasts are precisely those created manually. Nothing
fancier is required. But I am not sure what you want to do with them once
you have them.
On Tue, 30 Aug 2005,
Are you talking about the pdf() device? You never say so, but you do say
The font encoding error does not appear in the postscript
device driver, and it works fine.
This report is far too vague: what did you do and what reported the file
was corrupt? (Note that Acrobat reader is famous for
On Wed, 31 Aug 2005, Berwin A Turlach wrote:
G'day Nam-Ky,
NKN == Nam-Ky Nguyen [EMAIL PROTECTED] writes:
NKN I intend to burn some R CDs to colleagues in Vietnam. I want
NKN to put all binary files for base as well as contributed
NKN packages (for both Windows and Linux).
I
On Tue, 30 Aug 2005, b wrote:
Solved the problem with the 'R' installation!
One of the processes that loads up on startup was
conflicting with it. It was a program called
Gameutil.exe which seems to enhance games by
altering display modes. For some reason a conflict
existed. Anyhow I dont
On Wed, 31 Aug 2005, Philippe Grosjean wrote:
deggle wrote:
Hello,
I'm very new in working with tcl/tk in R and have a problem which will
probably
sound silly to most of you.
Here is the code I have problems with:
readcelfiles - function()
{
require(tcltk)
tt -
On Wed, 31 Aug 2005, Duncan Murdoch wrote:
Prof Brian Ripley wrote:
On Wed, 31 Aug 2005, Philippe Grosjean wrote:
deggle wrote:
Hello,
I'm very new in working with tcl/tk in R and have a problem which will
probably
sound silly to most of you.
Here is the code I have problems
On Thu, 15 Sep 2005, Barry Rowlingson wrote:
Prof Brian Ripley wrote:
substring(str, c(1, 26), c(25,length(str)))
nchar(str) surely?
Yes, or anything larger: I actually tested 1.
regexps can be rather slow though. Here's two functions:
But that's not the way to do this repeatedly
The problem is a known one and solved in R 2.1.1-patched. On your OS the
result is not as accurate as most, but the tolerance set was too tight so
the test failure is not something to worry about.
Please install R-patched instead, as it has many bug fixes in place.
On Thu, 15 Sep 2005,
On Fri, 16 Sep 2005, Dieter Menne wrote:
Charles Annis, P.E. Charles.Annis at StatisticalEngineering.com writes:
Here's an example with a glm; lm() works the same way but has fewer internal
objects.
mod3 - glm(tree ~ altitude, family = binomial)
You can use names() to find out what's
You send HTML mail (please see the posting guide and do as we ask) and it
seems to have not produced a working text version.
My guess is that you have a changed version of USArrests around. In any
case, there is no problem in R 2.1.1 (sic), so this is local to you.
On Fri, 16 Sep 2005,
Try to give a vector result with one element for each of element of x,
e.g.
integrandtotest - function(x) colSums(outer(t, x, ^))
works, although in fact this integration can be done analytically (it is a
sum of exponentials).
On Fri, 16 Sep 2005, A.Brennan wrote:
Hi
i am having a
See the rw-FAQ ..., Q2.8 in the latest version.
On Sat, 17 Sep 2005, Srinivas Iyyer wrote:
Dear group,
apologies if this is a stupid question. I searched
CRAN sites. I am afraid I missed it.
Can any one help me if I can update my windows version
of 1.9.1 to 2 or higer.
--
Brian D.
On Mon, 19 Sep 2005, Chris Buddenhagen wrote:
What I would like to know is if there is a way to know which specific
samples fall within a terminal branch created by tree? The classification
summarizes data, but I want to know which specific samples fall within each
classification (branch in
Your message *was* received, and you can check the archives to see it at
https://stat.ethz.ch/pipermail/r-help/2005-September/077889.html
You need to read the code to answer the question for yourself. There is
lots of code interpreting csplit in the rpart package. These lines might
be a
On Wed, 21 Sep 2005, Antonio, Fabio Di Narzo claimed:
This code gives an error:
a - ts(1:10, start=0, freq=10)
b - ts(1:10, start=1, freq=10)
ts.intersect(a,b)
In reality, it gives a helpful warning:
ts.intersect(a,b)
NULL
Warning message:
non-intersecting series in: .cbind.ts(list(...),
On Wed, 21 Sep 2005, Robin Hankin wrote:
If I do this:
par(mfrow=c(2,2))
persp(matrix(1:4,6,6),box=F,phi=33,theta=33)
persp(matrix(1:4,6,6),box=F,phi=33,theta=33)
persp(matrix(1:4,6,6),box=F,phi=33,theta=33)
persp(matrix(1:4,6,6),box=F,phi=33,theta=33)
(
or indeed
This really is the wrong list. Please do read the posting guide as we do
ask (before posting). In so far as this is relevant to R at all, it is
within the description of the R-devel list.
On Wed, 21 Sep 2005, Bogner, Konrad (LfU) wrote:
I'm trying to load a dynamic link library and it seems
The problem was that a tolerance in window.default was abs(start)*ts.eps,
which is no tolerance at all if start = 0.
So the special feature here was that the intersection of the series
started at zero. I had already tested and committed a fix.
window.default does use a tolerance for en, just
This was not really accurate.
Being `read-only' is a Windows phenomenon on open files (even files open
for reading). In 2.0.1 and 2.1.1 the file was open whilst it was parsed
and so `read-only' during that time. In 2.0.1 it was then closed, but in
2.1.1 it was (potentially) opened several
Concise answers to verbose questions:
1) Use options(scipen) (and probably change the margin sizes).
Or something like
options(scipen=10)
par(mar=c(5,8,4,2)+0.1)
plot(x, y, axes=FALSE)
axis(2, las=2)
axis(1, labels=FALSE)
axis(1, at = c(-2e8, 2e8), labels = expression(-2 %*% 10^8, 2 %*% 10^8))
The original reply was deliberately (I guess) vague. (I've removed the
history, as attributions had already been removed, in violation of
copyright law. If you cite someone, you MUST credit the author.)
Sometimes a little knowledge is a dangerous thing, and we have had a
number of partially
How to generate an import library for R.dll for Visual Studio is discussed
in README.packages: please do read it.
Do note the recommendation not to install R in a path with spaces - your
tool is one of many that does not work with spaces. But you don't want a
Borland import library (any more
This is just a question of rounding error. You are computing 11/999,
which is not exactly representable. Most computers use IEC60559
arithmetic, which gives a precision of .Machine$double.eps ~ 2e-16. So
you can expect rounding error around 11/999 * .Machine$double.eps ~ 2e-18.
However,
datai[1,1] appears to be a 'POSIXct' object. Which date that is depends
on the locale, and as.Date uses UTC (see ?as.Date). For me:
d - as.POSIXct(2005-07-01)
d
[1] 2005-07-01 BST
format(d, tz=GMT)
[1] 2005-06-30 23:00:00
as.Date(d)
[1] 2005-06-30
Use
as.Date(d + 23.99*3600)
to avoid
On Fri, 23 Sep 2005, Thomas Lumley wrote:
On Wed, 21 Sep 2005, Martin Maechler wrote:
Dirk == Dirk Enzmann [EMAIL PROTECTED]
on Wed, 21 Sep 2005 13:18:32 +0200 writes:
Dirk The selection problem can be solved by
Dirk dr2000=read.spss('myfile')
Dirk
On Sun, 25 Sep 2005, Chun-Ying Lee wrote:
Dear R user:
I bulid a package, and in the package I use the function nls
to solve some questions. If I have two sets of data, and I want to
summary these two data's nls output, I write the command in the
package source code like:
{
..
On Sun, 25 Sep 2005, Andrew Robinson wrote:
Hi R-helpers,
I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R
Version 2.1.1.
I have constructed a package using package.skeleton(), when I try
$ R CMD build foo
* checking for file 'foo/DESCRIPTION' ... OK
* preparing
This works if you omit the deriv() step.
Use R's options(error=dump.frames) and debugger(). This gives
Browse[1] rhs
[1] 0.433 0.4272571 0.3994105 0.3594037 0.3270730 0.3104752 0.3000927
[8] 0.2928445 0.2874249 0.2831787
attr(,gradient)
VrVm alphalamda
First a warning: loess in R is only loosely related to loess in S, being
derived from a C implementation (by the same authors).
In R I don't think you can do this. Those details are never exposed, and
are hidden in an undocumented C/Fortran workspace.
On Mon, 26 Sep 2005, Carlisle Thacker
On Mon, 26 Sep 2005, Martin Henry H. Stevens wrote:
Hello all,
1. Does Matrix 0.98-7 fix any of this?
2. Assuming no, how does one acquire Matrix 0.95-13?
It is in the Archive on CRAN, e.g.
http://cran.r-project.org/src/contrib/Archive/M/Matrix_0.95-13.tar.gz
Cheers, and thank you kindly
On Mon, 26 Sep 2005, Dean Sonneborn wrote:
I am attempting to read in a SAS 9.1 data file. After starting R I
change to the directory containing the sas data file and use the dir
command to confirm that it is there. Then I run the following R-code:
library(foreign)
sashome - /Program
You will have to show us the error! It will be shown in
d-p-q-r-tests.Rout.fail (unless this was a segfault or similar).
It is not OK to skip the test, but note that this test is random and does
fail about 1 in 50 times, so you could just try rerunning it.
On Mon, 26 Sep 2005, Jeff Ross
This is what happens if you don't have a usable X11 display. Did you
perhaps use a root account on a console owned by a normal user?
On Tue, 27 Sep 2005, Fernando Mayer wrote:
Dear R-users,
i'm a very newbie in linux, but decided to build R from source.
Following the R Installation and
On Tue, 27 Sep 2005, Richard Nixon wrote:
Hi
Just one of those niggles...
I've just been trying to plot a filled circle.
I thought that this would do it
plot(1,1,type=n)
points(1,1,pch=1,bg=blue,cex=5)
#bg: background (fill) color for open plot symbols
But I need to do this instead
you
need the console running the user that is going to be running R.
(I see this when running tests remotely, as then I have no X11 display set
to save bandwidth. I just ignore it.)
Thanks,
Fernando Mayer.
Prof Brian Ripley escreveu:
This is what happens if you don't have a usable X11
Prof Brian Ripley [EMAIL PROTECTED] wrote:
You will have to show us the error! It will be shown in
d-p-q-r-tests.Rout.fail (unless this was a segfault or similar).
It is not OK to skip the test, but note that this test is random and
does fail about 1 in 50 times, so you could just try
I've not seen a reply to this, nor ever seen it.
Please make a reproducible example available (do see the posting guide).
On Sun, 25 Sep 2005, Ajay Narottam Shah wrote:
I have a situation with a large dataset (3000+ observations), where
I'm doing lags as regressors, where I get:
Call:
On Tue, 27 Sep 2005, Eric van Gyzen wrote:
On Sun, Sep 25, 2005 at 12:07:02PM +0100, Prof Brian Ripley wrote:
On Sun, 25 Sep 2005, Andrew Robinson wrote:
Hi R-helpers,
I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R
Version 2.1.1.
I have constructed a package using
On Wed, 28 Sep 2005, Robert Bagchi wrote:
Hi Patrick
thanks for your advice. I have now tried glmmPQL, and it worked fine -
I'm getting consistent results between plots and models fitted by
glmmPQL. Plus it allows predict() and resid() which is another advantage
over lmer at present.
quick
This is entirely a function of the database format the ODBC Excel driver
expects. It is *not* an error in odbcConnectExcel as your subject
accuses.
To use a spreadsheet as a database you need to have column names.
On Thu, 29 Sep 2005, [ISO-2022-JP] ÅÄÃæ Áï wrote:
Dear R-help
I would like
On Wed, 28 Sep 2005, Sérgio Nunes wrote:
Hi,
I'm pretty new to R, I've just installed version 2.1.1 on Windows.
I have a simple (it seems) doubt - how do I change the Console default
Language ?
It's set to Portuguese but I would like to view it in English.
Do read the administration manual
On Thu, 29 Sep 2005, Mike Cheung wrote:
Dear all,
Edmond Ng (http://multilevel.ioe.ac.uk/softrev/reviewsplus.pdf) provides
an example to fit the mixed effects meta-analysis in Splus 6.2. The
syntax is:
lme(fixed=d~wks, data=meta, random=~1|study, weights=varFixed(~Vofd),
On Thu, 29 Sep 2005 [EMAIL PROTECTED] wrote:
booop booop a écrit :
1...Kindly tell me is it possible to form
a matrix which contains a no of matrices.. for eg..
if a,b,c,d are matrices
and e is a matrix which contains a,b,c,d as rows and columns..
I don't think you can use
On Thu, 29 Sep 2005, Mulholland, Tom wrote:
Well I downloaded the data using the link in your message which suggests
that the code is right. I don't have its loaded (I assume it's from the
irregular time series package)
The package is called its, as the message correctly said.
so I can't
Note your first error message is
connections.c, line 2580.18: 1506-052 (S) Duplicate case label for value 4.
Labels must be unique.
and the relevant line is
#if SIZEOF_LONG == 8
case sizeof(long):
so you have somehow got options to give you 4-byete longs although
configure
On Sat, 1 Oct 2005, Ramón Casero Cañas wrote:
I'd like to update MASS from version 7.2-11 to version 7-2.19. I am
running R 2.0.1 on ubuntu, that installs MASS with the package r-cran-vr.
First, MASS is part of VR, so it is VR you update.
Second, please do *READ THE HELP PAGE*. The first
Use mtext() to put labels in the margin line you want. title() writes
xlab and ylab in line 3 AFAIR. So you could use something like
par(mai=c(0.6,0.6,0.2,0.2))
plot(D$dates, D$inflation,
type=l, lwd=2, col=cadetblue4, cex.axis=0.6)
mtext(X label string, 1, line=2, cex = 0.6)
mtext(Y
On Thu, 29 Sep 2005, Christian Hennig wrote:
?confint
Thank you to all of you.
As far as I see this is not mentioned on the lm help page (though I
presumably don't have the recent version), which I would
suggest...
and I would suggest that you study a good book on the subject.
(confint
On Fri, 30 Sep 2005, sloan jones wrote:
I have used the RODBC package to read in data I have stored in an Access
file. When I am using data from files other than ACCESS I have no
problem using the survival package to work with dates; however, with the
ACCESS data the dates are reading-in
On Sat, 1 Oct 2005, Haiyong Xu wrote:
I just want to ask if there is any function that can convert decimal
number to binary vector.
What do you mean by a binary vector? Function intToBits might be what you
are looking for, e.g.
x - as.integer(1234)
y - intToBits(x)
options(width=50)
y
On Sat, 1 Oct 2005, Paul Roebuck wrote:
On Sat, 1 Oct 2005, Marc Schwartz wrote:
On Sat, 2005-10-01 at 06:07 -0500, Paul Roebuck wrote:
Can someone tell me how to fix the left margin of plot region
such that the tick values don't overwrite the x-axis label?
I haven't been able to set the
On Sun, 2 Oct 2005, Kemp S E (Comp) wrote:
Hi,
I am using the arima.sim function to generate some AR time series.
However, the function does not seem to produce exactly the same time
series when I specify the innov parameter. For example
r - rnorm(300)
x - arima.sim(300,
On Sun, 2 Oct 2005, Paolo Cavatore wrote:
Does anyone know whether R has its own modeling language for optimization
problems (like SIMPLE in NuOPT for S-plus)?
No. Note that SIMPLE is the language of NUOPT, not of S-PLUS. There is
an (extra-cost) interface module S+NUOPT, but it is an
I think this an issue about the amount of graphics memory. You are asking
for an image of about 17*2*3 = 102Mb, and you need more than that.
From the help page:
Windows imposes limits on the size of bitmaps: these are not
documented in the SDK and may depend on the version of
On Sun, 2 Oct 2005, Jonathan Baron wrote:
Try
R CMD INSTALL -l lib pkgs
The help file is in the utils package. I'm sure this is
documented in the manual too.
But you still need to ensure that library() can find the packages.
I think the simplest way is that we have set up by default for
Either use \dontrun in a \examples{} section of a help page (see Writing
R Extensions), or something like
num_reps - if(interactive()) readline(How many reps do you have... ) else 500
On Mon, 3 Oct 2005, Ken Termiso wrote:
Hi all,
I've got a package i've written that i am trying to check,
] On Behalf Of Prof Brian Ripley
Sent: Sunday, October 02, 2005 10:38 AM
To: Paolo Cavatore
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] modeling language for optimization problems
On Sun, 2 Oct 2005, Paolo Cavatore wrote:
Does anyone know whether R has its own modeling language for optimization
On Mon, 3 Oct 2005, Tuszynski, Jaroslaw W. wrote:
It is unclear what you are trying to do, but check bin2raw in caTools
package:
print(x - (1:5)*pi)
[1] 3.141593 6.283185 9.424778 12.566371
[5] 15.707963
print(y - bin2raw(x))
[1] 18 2d 44 54 fb 21 09 40 18 2d 44 54 fb 21 19
On Mon, 3 Oct 2005, Daniel Pick wrote:
I have successfully downloaded the sources and built R as a shared
library on a Red Hat Enterprise Level 3 box. I am now trying to build the
GNOME GUI, but configure is barfing on glade. According to the system
logs, the RPM for
help.search(convex programming)
No help files found with alias or concept or title matching 'convex
programming' using fuzzy matching.
...
Thanks,
Roger
On 10/3/05, Prof Brian Ripley [EMAIL PROTECTED] wrote:
On Mon, 3 Oct 2005, Huntsinger, Reid wrote:
Have you looked at the R interface
On Mon, 3 Oct 2005, Marc Schwartz (via MN) wrote:
On Mon, 2005-10-03 at 22:00 +0100, Ted Harding wrote:
Rolf ( Marc)
On 03-Oct-05 Marc Schwartz (via MN) wrote:
On Mon, 2005-10-03 at 16:31 -0300, Rolf Turner wrote:
A student in one of my courses has asked me about getting R graphics
output
On Tue, 4 Oct 2005, Weiwei Shi wrote:
Hi, Jim:
I tried your code and get the following error:
trn1-read.table('trn1.svm', header=F, na.string='.', sep='|')
Med-apply(trn1, 2, median, na.rm=T)
Ind-which(is.na(trn1), arr.ind=T)
trn1[Ind]-Med[Ind[,'col']]
Error in [-.data.frame(`*tmp*`, Ind,
On Wed, 5 Oct 2005 [EMAIL PROTECTED] wrote:
Hi,
Good Day Everybody.
Can anybody help me
I am working with R ...and i am using nnet() function for regression, but
after trained model, i m getting all predictions as 0.99 or 1.,
why is so, i dont understand?
As a guess, you
On Wed, 5 Oct 2005, Leonardo Lami wrote:
Hi,
I have the same problem.
I found a solution but I think there is something of more simple and correct.
I take a group and I put 0 the values of the other group, after I use the
function glm to abtain the Fisher's discriminant function for this
On Wed, 5 Oct 2005, Robin Hankin wrote:
I have a square matrix Ainv of size N-by-N where N ~ 1000
I have a rectangular matrix H of size N by n where n ~ 4.
I have a vector d of length N.
I need X = solve(t(H) %*% Ainv %*% H) %*% t(H) %*% Ainv %*% d
and
H %*% X.
It is possible to
On Wed, 5 Oct 2005, Duncan Murdoch wrote:
On 10/5/2005 9:31 AM, Prof. Paul R. Fisher wrote:
Hi all
I am a complete newbie to this list (just subscribed) and a newcomer to
R (an S user from olden times). I have been using scatter3d to create a
3d scatter plot with surface. The graphic is
On Thu, 6 Oct 2005, Luc Vereecken wrote:
Version 2.x no longer has a -gui=none option. Running BATCH jobs
supposedly disables a gui though.
There is no --gui=none option in 2.1.x (but there is in 2.0.x).
The gui is not disabled in 2.1.x either according to the help page nor in
practice
# Not sure why this did not work
#sapply(tw, setlen, len)
It probably did, but you discarded the result. Try
tw - sapply(tw, setlen, len)
On Thu, 6 Oct 2005, sosman wrote:
Marc Schwartz wrote:
On Wed, 2005-10-05 at 22:19 +1000, sosman wrote:
I have some data in a CSV file:
On Thu, 6 Oct 2005, Omar Lakkis wrote:
Can someone, please, explain the difference is results below (notice
the isdst value)
Timezones. (That is carefully documented on the help page for
as.POSIXlt.)
unlist(as.POSIXlt('2005-7-1'))
sec min hour mday mon year wday yday isdst
0
Your R is an obselete version. There are some changes in 2.2.0 which
affect the hinting given to X11 WMs, so please try the current version of
R. (I know it had only recently been released, but the posting guide
refers you to the development versions.)
I don't know if this will solve it (I
On Fri, 7 Oct 2005, Robin Hankin wrote:
Hi
I tried to compile R-2.2.0 just now.
On MacOS X.
configure worked fine, but
compilation stopped with
Please ask MacOS-specific questions on the r-sig-mac list.
(Missing RestFP, SaveFP.)
I believe this is well-known and indicates an incompatible
Please tell us which version of R.
If this is 2.2.0 (I think it is) you need to review the tools you are
using, as mingw-runtime-3.8 is required. (It says so explicitly in that
manual, and that is the current version of that package.)
For R-patched you will also need a recent enough version
On Fri, 7 Oct 2005, Andreas Cordes wrote:
I have the problem that for the step procedure stops due to missing
values. There are no options in Step or stepAIC to handle missing
values. Is there any way to run stepwise modelselection in R in an
automated way in this case?
Try the hint it gives
On Fri, 7 Oct 2005, Dimitri Giunchi wrote:
Hi all,
I was unable to obtained a smoothed line using the loess function.I used the
following code reported in the examples of R documentation:
cars.lo - loess(dist ~ speed, cars)
Then I tried to plot both the data and the smoothed line
Why not use a text connection?
On Fri, 7 Oct 2005, Paul Roebuck wrote:
I have a data file from which I need to read portions of
data but data location/quantity can change from file to file.
I wrote some code and have a working solution but it seems
wasteful to have to do it this way. Here's
If A is the matrix the answer is 2*(1 - 2^(-A)), which took about 10secs
for an example of your size.
From \sum_{i=1}^n x^{1-i} = (1-x^{-n})/(1-x), E OE.
On Fri, 7 Oct 2005, Tim Smith wrote:
I have a 7000x7000 matrix, and each element is an integer. For each
element, I want to apply the
On Fri, 7 Oct 2005, Prof Brian Ripley wrote:
If A is the matrix the answer is 2*(1 - 2^(-A)), which took about 10secs
for an example of your size.
From \sum_{i=1}^n x^{1-i} = (1-x^{-n})/(1-x), E OE.
NB: I have assumed n = 1 here, as people nornally do when using 1:count.
On Fri, 7 Oct
On Fri, 7 Oct 2005, Patrick Burns wrote:
If there weren't an analytic solution to your problem,
then you could build a vector of the answers from 1
to the maximum in the matrix. Call that 'wtvec'. Then:
ans - array(NA, dim(A), dimnames(A))
ans[] - wtvec[as.vector(A)]
should get you what
On Sat, 8 Oct 2005, Marc Schwartz wrote:
On Sat, 2005-10-08 at 09:28 -0400, Denis Chabot wrote:
Hi,
A few times I tried to control the number and position of tick marks
in plots with the yasp or xasp parameters. For example, a y axis was
drawn by default with tick marks at 0, 20, 40, 80 and
The example on the png help page, Rplot%03d.png, _is_ what you want.
(More details are on ?postscript.)
formatC() is an S/R peculiarity: sprintf() is the cross-language way to do
this sort of thing.
On Sat, 8 Oct 2005, bogdan romocea wrote:
Dear useRs,
Is there a way to 'properly' format %d
On Sat, 8 Oct 2005, Denis Chabot wrote:
Hi, sorry about the bad syntax, though the right syntax would not have worked
either, according to your tests (Mark, Brian, Peter).
It DOES work according to my tests! (Do give us the credit for testing
our advice: we would appreciate your showing
On Sat, 8 Oct 2005, Ajay Narottam Shah wrote:
I am doing
library(rpart)
m - rpart(y ~ x, D[insample,])
D[outsample,]
y x
8 0.78391922 0.579025591
9 0.06629211 NA
10 NA 0.001593063
p - predict(m, newdata=D[9,])
Error in model.frame(formula,
On Sat, 8 Oct 2005, Erin Hodgess wrote:
Dear R People:
I have a question about the lm.ridge function, please.
In the example, there is one set of output values in the select
function but another in the comment section.
Am I missing something please?
The values in the examples were
On Mon, 10 Oct 2005, Bernd Weiss wrote:
I successfully set up a local MySQL-database. Connecting via RODBC is
not problem, the same in fetching 3 of 4 tables. But trying to
connect to table 4 fails.
author-sqlFetch(test,author)
Error in fromchar(unclass(x)) : character string is not in a
On Mon, 10 Oct 2005, Philippe Grosjean wrote:
Hello,
You must explicitly use print(), show() on an object -here, use
print(i)- in a loop or alternatively, use cat() to display string like:
cat(loop, i, \n)
With RGui under Windows, there is another subtility: if you have turn on
'Misc -
201 - 300 of 4745 matches
Mail list logo