On Thu, 12 Oct 2006, Gavin Simpson wrote:
> On Wed, 2006-10-11 at 22:29 -0500, Frank E Harrell Jr wrote:
>> x <- c(0,7,8,14,15,120,242)
>> y <- c(122,128,130,158,110,110,92)
>>
>> lowess(x,y)
>>
>> $x
>> [1] 0 7 8 14 15 120 242
>>
>> $y
>> [1] 122. 128. 132.2857 158.
This issue is discussed in much more detail in the manuals for 2.4.0, in
particular in the new 'R Internals' manual. You will find the basic rules
in the Blue Book (Becker, Chambers, Wilks, 1988) but they have not always
been applied consistently.
On Wed, 11 Oct 2006, Michael Toews wrote:
> H
On Wed, 11 Oct 2006, Caroline TRUNTZER wrote:
> Dear all,
> I would like to estimate the parameters of a mixture of two beta
> binomial distributions. I saw the function fitdistr from the MASS
> library could be a solution but I have some difficulties. Does anyone
> know how to specify to fitdistr
On Wed, 11 Oct 2006, Stefan Grosse wrote:
> Joe Byers schrieb:
>> JGR, when I manually compiled it for my system by specifying my location
>> of java 1.5.0_06, worked fairly nicely. The autoinstall does not find
>> my installation of this java version since Redhat EL4 does not support
>> it. Red
On Wed, 11 Oct 2006, Dale Steele wrote:
> I'm trying to write a small function (below) to compute Box & Cox
> transformations of x for arbitrary values of lambda.
It seems to compute the log-likihoods of transformed models.
> I'd like to
> specify a range of values for lamba (min,max,step) and a
On Tue, 10 Oct 2006, Yury Yuryev wrote:
> Dear friends,
>
> I am sorry, I again rise that boring question about font embedding in
> EPS figure.
It gives us an opportunity to point you to the solution included in R
2.4.0.
> I found some discussions on this topic but there were no strait solution
On Tue, 10 Oct 2006, January Weiner wrote:
> Oh, thanks, that was hint enough :-) I see it now. I turns that R does
> not understand
>
> e-10
>
> ...which stands for 1e-10 and is produced by some of the bioinformatic
> applications that I use (notably BLAST).
And that is not standard C notation.
On FC5 Linux:
gannet% cat > foo.dat
a 1 2e-4
b 2 3e-8
gannet% R
...
> read.table("foo.dat")
V1 V2V3
1 a 1 2e-04
2 b 2 3e-08
> sapply(read.table("foo.dat"), class)
V1V2V3
"factor" "integer" "numeric"
so please tell us your environment and give a reproducible ex
Which Solaris is this? I am not seeing any problems with our Solaris 8
systems. Was the compiler built on that exact system? What gcc options
were used (-std=gnu99, as recommended, for instance)?
It looks like a system header problem: __builtin_isnan ought to be
built in, not an external ent
On Fri, 6 Oct 2006, Christian Bieli wrote:
Hi all
I have to generate some test data for import in an sql database. The
database is meant for web-based data entry in a study taking place in a
german speaking region, so factor levels of the variables include umlauts.
The variables in the datafram
On Fri, 29 Sep 2006, Sarosh Jamal wrote:
> I've been getting the following warning on each package update/install (Rv2.4
> on SunOS9):
Each? This only applies to a few packages, and is only a warning.
> Warning message:
> 'DESCRIPTION' file has 'Encoding' field and re-encoding is not possible.
On Mon, 9 Oct 2006, Mstislav Elagin wrote:
> Hi, All,
>
> I have come across a problem building a DLL from .f90 source (R 2.3.1,
> Windows XP). When using the R CMD SHLIB procedure, the DLL itself was
> being built, but its export table was empty.
>
> Among the output from R CMD SHLIB the followin
On Wed, 4 Oct 2006, paul sorenson wrote:
> What is the recommended way to convert/coerce and integer to a POSIXct
> please?
>
> d <- as.POSIXct(Sys.Date())
> i <- as.integer(d)
>
> as.POSIXct(i)
> Error in as.POSIXct.default(i) : do not know how to convert 'i' to class
> "POSIXlt"
>
> This appears
On Fri, 29 Sep 2006, Søren Højsgaard wrote:
I write a package foo which requires a package bar (from CRAN) to work.
So in the DESCRIPTION file I write Depends: bar. I would like it to be
so that when one installs foo, then it is automatically checked whether
bar is installed, and if not then b
On Thu, 28 Sep 2006, David Kaplan wrote:
> I can't seem to figure out why I'm getting this error. The output is
> copied right off the screen. Notice how in some cases the back slash is
> missing. In other cases, it can't read a file that I know is there.
This is rw-FAQ 2.16 and FAQ 7.8
Since
On Thu, 28 Sep 2006, Joe W. Byers wrote:
> Berton Gunter wrote:
>> Please folks -- use indexing.
>>
>> myframe<-myframe[,c(1,5,2,3,4)]
>>
>> Which begs the question: why bother rearranging the columns anyway, since
>> one can get them used, printed, etc. in any order you wish anytime you want
>> j
On Thu, 28 Sep 2006, hadley wickham wrote:
>>
>> And, to answer the specific question: Yes, R has lazy evaluation,
>> everywhere. Arguments are always evaluated if and when they are
>> needed.
>>
>
> But doesn't R has a rather limited force of lazy evaluation? - you
> have no control over it, apar
On Thu, 28 Sep 2006, Hans-Peter wrote:
> Hello,
>
> What's the best/recommended way to drop the dimension of a data.frame?
> Probably "test[,,]" but I am not sure. I would prefer an explicit
> function, but drop doesn't work with frames and I didn't find
> something in the helpfiles/mailarchive.
Jeff
>
> On Tue, 2006-03-28 at 17:00 +0100, Prof Brian Ripley wrote:
>> Rgui now supports \r in the same way as rterm.
>
>
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford,
matrix(some_list, nr, nc, byrow=TRUE) may be what you are looking for.
R arrays can be of any vector type, including list. I'd get used to R's
Fortran ordering rather than force transposes all the time.
On Wed, 27 Sep 2006, Evan Cooch wrote:
> In the process of moving a number of my scripts f
On Tue, 26 Sep 2006, Richard M. Heiberger wrote:
> I like this plan and have now played with the concept. I did the following
> on Windows in cygwin. It would also work in Unix, and I think could be
> tickled
> to work on the standard MS cmd line in Windows. It would certainly work
> on Window
On Tue, 26 Sep 2006, Giuseppe Antonaci wrote:
> I think this is not a R-devel question. Sorry to all if I'm wrong,
> please let me know.
In what sense is this not a programming question?
> I managed to build R successfully with the default BLAS but when I
> change the MKRULES to use ATLAS BLAS a
On Tue, 26 Sep 2006, [EMAIL PROTECTED] wrote:
>
> Original-Nachricht
> Datum: Tue, 26 Sep 2006 09:56:53 +0100 (BST)
> Von: Prof Brian Ripley <[EMAIL PROTECTED]>
> An: [EMAIL PROTECTED]
> Betreff: Re: [R] rpart
>
>> On Mon, 25 Sep 2006, [EMAIL PRO
On Mon, 25 Sep 2006, Justin Rhodes wrote:
> Dear R-help,
>
> I recently tried implementing glmmPQL in 2.3.1,
I thought *I* had implemented it: are you talking about my function in
package MASS or your own implementation?
> and I discovered a few differences as compared to 2.2.1.
You appear to
On Mon, 25 Sep 2006, [EMAIL PROTECTED] wrote:
> Dear r-help-list:
>
> If I use the rpart method like
>
> cfit<-rpart(y~.,data=data,...),
>
> what kind of tree is stored in cfit?
> Is it right that this tree is not pruned at all, that it is the full tree?
It is an rpart object. This contains both
On Mon, 25 Sep 2006, Bruce LaZerte wrote:
> # R version 2.3.1 (2006-06-01) Debian Linux "testing"
>
> # Is the following behaviour a bug, feature or just a lack of
> # understanding on my part? I see that this was discussed here
> # last March with no apparent resolution.
Reference? It is the th
I've not seen an actual answer to this, which is that this is a
misunderstanding as to how NextMethod works.
First,
> + x <- unclass(x)
looks wrong. NextMethod uses the next method at the call to the generic,
and subsequent changes to the object 'x' do not alter the class that would
be d
On Mon, 25 Sep 2006, Gad Abraham wrote:
> Hi,
>
> It seems to me that the Mersenne-Twister PRNG can be initialised using
> one integer instead of 624 integers, since inside RNG.c code there's a
> function defined as MT_sgenrand(Int32).
>
> How do I actually set this seed within R?
set.seed(), on
On Sat, 23 Sep 2006, Luca Telloli wrote:
> Hello R-Users,
> I'm new to R so I apologize in advance for any big mistake I might
> be doing. I'm trying to fit a set of samples with some probabilistic
> curve, and I have an important question to ask; in particular I have
> some data, from which
On Sun, 24 Sep 2006, Taka Matzmoto wrote:
> Hi R users
>
> Which version of S does the current R version use? Venables and Ripley's book
> (2000) says R is version 3 of S language. Is there any change since 2000 ?
Yes, the 'methods' package implements much of the changes in S version 4,
and the
On Fri, 22 Sep 2006, michael papenfus wrote:
> I am using Windows XP and R 2.3.1.
> During my lastest session I updated my packages and now when I try to start
> R 2.3.1
> I get the following error message:
>
> Fatal error: unable to restore saved data in .RData in an error window and
> Error in l
On Fri, 22 Sep 2006, [EMAIL PROTECTED] wrote:
> Hi Folks,
>
> I've just come across a kind of problem which leads
> me to wonder how to approach it in R.
>
> Basically, each a set of items is subjected to a series
> of "impacts" until it eventually "fails". The "force"
> of each impact would depen
On Fri, 22 Sep 2006, Duncan Murdoch wrote:
> On 9/22/2006 10:29 PM, X.H Chen wrote:
>> Dear R users,
>>
>> I have a question about R installation under Cygwin. When running
>> ./configure, I can't pass the checking phase due to an error:
>> --with-readline=yes (default) and headers/libs are not av
On Fri, 22 Sep 2006, Thomas Lumley wrote:
> On Fri, 22 Sep 2006, Thilo Kellermann wrote:
>
>> Hi,
>>
>> the data of the model fit is stored in lm$model and should work
>>
>
> Not reliably. In the first place, you should use the accessor function
> model.frame(model) rather than model$model, wh
On Fri, 22 Sep 2006, Tova Fuller wrote:
> Hello all,
>
> I have read as many merge issues as I possibly could tonight and
> although I presume this is a small error, I have not found the
> solution to my problem.
>
> I'm trying to merge two data sets: dat0 and TransTable. As you can
> see below,
The way nls looks for variables in its formula is very unusual. For
non-parameters it has
for(var in varNames[!varIndex])
mf[[var]] <- eval(as.name(var), data)
inside nls, and that means that the scoping rules are first 'data' then
the body of nls(). That's surely not what one wo
On Wed, 20 Sep 2006, Jeffrey Horner wrote:
> Gamal Azim wrote:
>> Is it possible to remotely save an R session then
>> terminate R? Of course the destructive task after
>> 'then' is rather straightforward by itself.
>
> Yes. Sending the process a SIGUSR1 signal:
>
> $ kill -USR1 pid
>
> will save
On Thu, 21 Sep 2006, Thorsten Muehge wrote:
>
> Hello R Experts,
> I want to aggregate parameters by week. But our production week ends Friday
> night instead of Sunday Night which is the default value in R.
The default in ISO8601, not just in R, but that is %W, not %U as used
below.
> In order
On Wed, 20 Sep 2006, Dylan Beaudette wrote:
> Getting a very strange error with a new install of R from source on x86;
>
> make[3]: Leaving directory `/tmp/R.INSTALL.r20887/cluster/src'
> ** R
> ** data
> ** moving datasets to lazyload DB
> Error in factor(c(1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,
On Wed, 20 Sep 2006, Dina Said wrote:
> Hi,
>
> I'm using fedora core 3. I downloaded R and installed it.
> After that, I invoked R by the root user and write the
> command install.packages("Rcmdr", dependencies=TRUE). This
> command installs Rcmdr automatically. However, whenever I
> type library
On Wed, 20 Sep 2006, Uwe Ligges wrote:
> Poizot Emmanuel wrote:
>> R under linux (Ubuntu or other distributions) is not exactly as it is
>> under windows.
>> Under linux, you have to use a shell (or a terminal) to launch R and use
>> it.
>> The main difference between windows and linux version of
On Wed, 20 Sep 2006, Murray Jorgensen wrote:
> I am in a discriminant analysis situation with a frame containing
> several variables and a grouping factor, if you like:
>
> set.seed(200906)
> exampledf <- as.data.frame(matrix(rnorm(50,5,2),nrow=10,ncol=5))
> exampledf$Group <- factor(rep(c(1,2,3),
On Tue, 19 Sep 2006, Gavin Simpson wrote:
> Dear list,
>
> Following advice posted to this list a while back by Prof Ripley [1], I
> have been trying to draw a per mille character [2] in an axis label.
>
> This should give the correct character:
>
> plot(1:10, ylab = "\u2030")
>
> but all I get is
On Tue, 19 Sep 2006, Sean O'Riordain wrote:
> Hi Duncan,
>
> Thanks for that. In the light of what you've suggested, I'm now using
> the following:
>
> # generate a random integer from 0 to t (inclusive)
> if (t < 1000) { # to avoid memory problems...
>M <- sample(t, 1)
> } else {
>
On Mon, 18 Sep 2006, Peter Dalgaard wrote:
> "Sarosh Jamal" <[EMAIL PROTECTED]> writes:
>
>> Hi there,
>>
>> I was updating the R-cmdr add-on (v.1.1-6 to the latest v.1.2) for R
>> (v.2.2.0) in a SunOS9 environment and came across some warnings during my
>> installation - it seems to download the
You need to print() explicitly inside sourced-ed script: auto-printing is
not active.
On Mon, 18 Sep 2006, Ulrich Leopold wrote:
> Dear list,
>
> First my information:
> platform i386-pc-linux-gnu
> arch i386
> os linux-gnu
> system i386, linux-gnu
> status
>
See
?try
?tryCatch
On Mon, 18 Sep 2006, Stefan Th. Gries wrote:
> Hi
>
> I wrote a script which retrieves links from websites and loads them with scan:
>
> ...
> website<-tolower(scan(current.pages[i], what="character", sep="\n",
> quiet=TRUE))
> ...
>
> However occasionally, the script finds b
On Sun, 17 Sep 2006, jim holtman wrote:
You can use 'names'
d = as.data.frame(matrix(c(1, 2, 3, 3), 2,2))
x <- table(unlist(d))
x
1 2 3
1 1 2
names(x)
[1] "1" "2" "3"
You can, but Søren was equally correct: names() on a 1D array actually
looks up the dimnames. OTOH, we do really want
On Fri, 15 Sep 2006, Rainer M Krug wrote:
> Hi
>
> Linux SuSE 10
>
> platform i686-pc-linux-gnu
> arch i686
> os linux-gnu
> system i686, linux-gnu
> status
> major 2
> minor 3.1
> year 2006
> month 06
> day01
> s
On Thu, 14 Sep 2006, Werner,Arelia [PYR] wrote:
> I am trying to run a cross-correlation using the "ccf()" function. When
> I select plot = TRUE in the ccf() I get a graph which has ACF on the
> y-axis, which would suggest that these y-values are the auto-correlation
> values.
But cross-correlati
On Fri, 15 Sep 2006, Kiermeier, Andreas (PIRSA - SARDI) wrote:
> The figure margins come from what is set in par("mar"), eg
>
>> layout(matrix(c(1:10),5,2),heights=c(1,rep(2,4)))
>> par("mar")
> [1] 5.1 4.1 4.1 2.1
>>
>
> There is not enough space left to plot anything with those margins. You
> w
On Thu, 14 Sep 2006, Deepayan Sarkar wrote:
[...]
> Ah, I'd forgotten about parse. A link from ?expression might be reasonable.
But it says
Details:
'Expression' here is not being used in its colloquial sense, that
of mathematical expressions. Those are calls (see 'call') in R,
It's a known problem. Run the X11 device with the arguments stated,
as it says
On Tue, 12 Sep 2006, Gael Even wrote:
> Hi,
>
> I use R in remote access on a Debian server.
> I need X11() to create graphical things so I use xvfb-run (Virtual X server
> environnement) R to allowing X11() ca
On Mon, 11 Sep 2006, Gregor Gorjanc wrote:
> Prof Brian Ripley stats.ox.ac.uk> writes:
> ...
> > Check out tests/internet.R. nsl() checks if you can resolve host names,
> > which has worked well enough there.
>
> Thank you prof. Ripley for this pointer. I am pos
On Tue, 12 Sep 2006, Paul Smith wrote:
> Is there some way of exporting R plots to epslatex, i.e., to a file
> with the eps file and another one with the LaTeX commands
> (representing the text in the plots), likewise Gnuplot does? If so,
> could you please indicate it to me?
R has an xfig driver
On Tue, 12 Sep 2006, stephenc wrote:
> Hi
>
>
>
> I hope this isn't off topics, but I have always found when I stepAIC() some
> glm I get an improvement in accuracy and kappa, but I have just done a case
> where I got a marginal deterioration. Is this possible, or should I be
> going through
On Tue, 12 Sep 2006, Alexander Nervedi wrote:
> Hi
>
> I am having trouble with dotcharts, and I keep getting the error message:
>
> Error in Summary.data.frame(..., na.rm = na.rm) :
> only defined on a data frame with all numeric or complex variables
>
> I am sure there is a really sim
On Mon, 11 Sep 2006, John Kornak wrote:
>
> Thanks to Professor Ripley and Sarah Goslee for their helpful responses
> which have progressed my installation of rimage a little further.
> Installing fftw2 instead of fftw3 of course solves the problem and I
> apologize for initially missing that r
On Mon, 11 Sep 2006, James W. MacDonald wrote:
> Jeff Bricker wrote:
> > There's probably a better way, but what I do in this situation is to
> > capture output to a .pdf that I can page through.
>
> There are two ways that might be considered better, if better implies
> doing things from within
On Mon, 11 Sep 2006, John Kornak wrote:
> Dear R people,
>
> I am trying to install rimage using install.packages("rimage") but am
> receiving the following errors despite having FFTW installed. I would
> appreciate any help to get this fixed.
>
> checking fftw.h usability... no
> checking fft
On Mon, 11 Sep 2006, Dong Wang wrote:
>
>
> I want to make a shared object from Fortran code to be loaded with
> dyn.load(), but have the following messages:
> # g77 -c kerimp1.f
Where did that come from? It's not the appropriate command line, as the
message below very clearly tells you.
> #
As nnet is done almost entirely in compiled C, you may well find that
already most of the computation is in a compiled language.
Please look at `Writing R Extensions' and profile your code to find the
bottlenecks.
While you are looking at that manual, please also consider the section on
tidyin
On Mon, 11 Sep 2006, [EMAIL PROTECTED] wrote:
> On 11-Sep-06 Gregor Gorjanc wrote:
> > It seems that 'internal' method was used (I use R 2.3.1 under
> > Linux) as indicated in help page of download.file. I could
> > use wget or lynx methods, but these two must be available,
> > so this is not real
On Mon, 11 Sep 2006, Duncan Murdoch wrote:
> On 9/11/2006 8:20 AM, Prof Brian Ripley wrote:
> > On Mon, 11 Sep 2006, Duncan Murdoch wrote:
> >
> > > On 9/11/2006 6:15 AM, Robin Hankin wrote:
> > > > Hi Sundar
> > > >
> > > >
&
On Mon, 11 Sep 2006, Duncan Murdoch wrote:
> On 9/11/2006 6:15 AM, Robin Hankin wrote:
> > Hi Sundar
> >
> >
> > thanks for this. But I didn't make it clear that I'm interested in
> > extreme numbers
> > such as 1e300 and 1e-300.
That's not relevant, unless you are interested in extremely sm
On Sun, 10 Sep 2006, Maciej Blizi?ski wrote:
> Hello all R-help list subscribers,
>
> I'd like to create a regression tree of a data set with binary response
> variable. Only 5% of observations are a success, so the regression tree
> will not find really any variable value combinations that will
On Sun, 10 Sep 2006, Rainer Volz wrote:
> Tomas Willebrand szooek.slu.se> writes:
>
> > I have R Version 2.3.1 (2006-06-01) installed on a AMD 64 machine with
> > SUSE 10.1.
Installed how is the usual question: from sources or an RPM?
> > I have Sun Java version 1.5.0-sun installed.
> >
> >
You haven't told us your OS.
SQLBase would appear to have a Windows ODBC driver, so RODBC should be
usable. I doubt it any other solution is (and RSQLite works only with
SQLite).
On Fri, 8 Sep 2006, SPIESSER Vincent wrote:
> I am trying to extract data from a database with R in order to produc
On Thu, 7 Sep 2006, alex lam (RI) wrote:
> Dear list,
>
> I have been trying to run the function "qvalue" under the package qvalue
> on a vector with about 20 million values.
>
> > asso_p.qvalue<-qvalue(asso_p.vector)
> Error: cannot allocate vector of size 156513 Kb
> > sessionInfo()
> Version
Which version of R?
Please try 2.4.0 alpha, as it has a different and more efficient
algorithm for the case of 1-1 matches.
On Wed, 6 Sep 2006, Adam D. I. Kramer wrote:
> Hello,
>
> I am trying to merge two very large data sets, via
>
> pubbounds.prof <-
> merge(x=pubbounds,y=prof,by.x="user"
n Wed, 6 Sep 2006, Rolf Turner wrote:
> Prof. Brian Ripley wrote:
>
> > On Wed, 6 Sep 2006, Christos Hatzis wrote:
> >
> > > See ?sweep
> > >
> > > sweep(a, 2, a[1,],"/")
> >
> > That is less efficient than
> >
> > a/r
On Wed, 6 Sep 2006, Christos Hatzis wrote:
> See ?sweep
>
> sweep(a, 2, a[1,],"/")
That is less efficient than
a/rep(a[1,], each=nrow(a))
>
> -Christos
>
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
> [EMAIL PROTECTED]
> Sent: Wednesday,
cex works for me.
On Tue, 5 Sep 2006, Nair, Murlidharan T wrote:
> Which is the parameter that is used to decrease the size of ylabs
> plotted in biplot? I tried playing with cex and cex.lab I am not getting
> it right
>
> pc <- princomp(USArrests)
>
> biplot(pc, xlabs = rep("", nrow(USArrests
On Tue, 5 Sep 2006, Patricia Bautista wrote:
> Hi!, I wonder to know if someone can explain me how I
> can access R object code. Briefly, what I need to do
> is: given a user defined function in R I want to
> access its object code with C because I need evaluate
> the function but using C.
Are y
On Tue, 5 Sep 2006, Tong Wang wrote:
> Hi everyone,
> I am doing some intensive computation: 5 regressions of the
> form Y~X with each y of size (1,1000) , even if I break invoke gc() for
> a few time in the loop, it still breaks down at some point with the
> error message:
gc() doe
On Mon, 4 Sep 2006, Georg Otto wrote:
>
> Hi,
>
> I have a csv file where the number of filled columns varies in the
> different rows:
>
> Sun 5-Feb-06,15,,,01:30:00,0:06:00,
> Mon 6-Feb-06,,
> Tue 7-Feb-06,7,,,00:41:00,0:05:51,
> Wed 8-Feb-06,,
>
> I would like to use read.table (or w
On Mon, 4 Sep 2006, Björn Thalheim wrote:
> Hi,
>
> I noticed that, when working with R on a command line, I cannot enter
> anything into the command line which is longer than 1023 characters.
Actually, 1000 or 1022 or 1023 bytes, and it depends on your 'command
line' (and you have not even tol
Please do read the rw-FAQ, Q2.9 (and the posting guide).
In particular, Windows never gives 4GB to a single 32-bit user process.
On Sun, 3 Sep 2006, Davendra Sohal wrote:
> Hi,
> I'm using R on Windows and upgraded the computer memory to 4GB, as R was
> telling me that it is out of memory (for
?par, specifically mar and mai, xaxt and yaxt. Note that there are my
default on 2.1 lines in the right margin and 4.1 in the left one.
You should find the description in 'An Introduction to R' helpful.
On Sun, 3 Sep 2006, Chris Grobler wrote:
> Dear All,
>
>Having a bit of trouble with p
On Sun, 3 Sep 2006, Paul Smith wrote:
> Dear All
>
> I am trying to fit Pareto distribution to some data. MASS package does
> not support Pareto distribution. Is there some alternative way?
Actually fitdistr{MASS} does if you supply the pdf for a Pareto.
That is not in base R, but easy to write
On Sat, 2 Sep 2006, Rick Bischoff wrote:
> Infinite values?
Or NA or NaN values
> On Sep 2, 2006, at 3:48 PM, Nair, Murlidharan T wrote:
>
> > I am getting the following error when I an trying to use princomp
> > princomp(unique.data)
> > Error in cov.wt(z) : 'x' must contain finite values only
On Sat, 2 Sep 2006, Alexander Nervedi wrote:
> Hi
>
> I have
>
> >x<-matrix(c(1,2,3,4), ncol = 2)
> >x
> [,1] [,2]
> [1,]13
> [2,]24
>
> I'd like these two be divided by their column totals. so 1/3 and 2/3 and 3/7
> and 4/7. The obvious
> >x/colSums(x)
> [,1]
I believe people (including me) did not reply because you appeared not to
have done your homework. The help page for ?nls _does_ have a reference
to the 'port' documentation, and RSiteSearch("trust region") is
informative and leads to an R package that does trust-region optimization.
(So woul
Your example is not actually a regression, and not reproducible. Here is
one that is both:
> fm1 <- lm(weight ~ ns(height, df = 5), data = women)
> predict(fm1, newdata=women[1,], se=TRUE)
Error: variable 'ns(height, df = 5)' was fitted with class "nmatrix.5" but
class "nmatrix.1" was supplied
R's svd does this by default (at least according to Wikipedia's
definition). Take a closer look at the help page, and in particular 'nu'
and 'nv'.
On Fri, 1 Sep 2006, Lord Tyranus wrote:
> Hi wizards, I have seen the function svd of R for singular value
> decomposition, but I need to computes
It is the same issue. Martin's reply was for a specific value of 'a': the
factor is gamma(a).
Please do study the help page for pgamma (as the posting guide did ask you
to), as it has all the details.
Note that the Mathematica version can easily lead to representation
difficulties for even mo
Please use Rcmd INSTALL --build and see if that works for you. It is how
Uwe Ligges and I build the public repositories.
On just Rcmd INSTALL if you don't want to distribute it.
(R CMD and Rcmd are the same thing, and I prefer the older form which is
more efficient.)
On Fri, 1 Sep 2006, R Heb
install.packages() in R 2.4.0 will make a full path out of a relative
path, to avoid any confusion. (normalizePath is a good way to do that,
BTW).
On Fri, 1 Sep 2006, Rolf Turner wrote:
> Prof. Brian Ripley solved the problem. He wrote:
>
> > I was not aware that this works w
On Fri, 1 Sep 2006, Hans-Peter wrote:
> Hi,
>
> I output multiple "grid-package-based" plots to the postscript device.
> Because the graphics are complicated and consists of a lot of datapoints
> (~200'000) the files become really big. To avoid this big files and to
> shorten the creation, I curr
This warning is indeed really an error, but do you want all your downloads
to fail just because one does? install.packages() behaves the same way on
Unices.
I was not aware that this works with relative paths for any version of R.
Try using a full path, which always works for me.
If indeed you
On Fri, 1 Sep 2006, [EMAIL PROTECTED] wrote:
> Prof Brian Ripley wrote:
> > I would not have expected glm to be more than say 5x slower than lm if
> CPU
> > cycles and not memory were the limiting factor. In that case more RAM
> > might be all you need.
>
> The
On Fri, 1 Sep 2006, [EMAIL PROTECTED] wrote:
> Prof Brian Ripley wrote
> > 3) As I recall, you were doing model selection via AIC on 20,000
> >observations. You might want to think hard about that, since AIC is
> >designed for good prediction. I would do model
On Fri, 1 Sep 2006, Duncan Murdoch wrote:
> On 8/31/2006 10:40 PM, Kimpel, Mark William wrote:
> > I installed GSView and the file opens correctly and the output is as is
> > should be. As you suggest, this must be a bug in Adobe CS2 Illustrator
> > and Photoshop.
>
> It could conceivably still b
On Fri, 1 Sep 2006, [EMAIL PROTECTED] wrote:
> Peter Dalgaard wrote
> > Is this floating point bound? (When you say 30 factors does that mean
> > 30 parameters or factors representing a much larger number of groups).
> > If it is integer bound, I don't think you can do much better than
> > increas
On Wed, 30 Aug 2006, Hans-Joerg Bibiko wrote:
> If you are using 'only' English then
>
> str <- "dog"
> strsplit(str,NULL)[[1]]
>
> works perfectly and it is fast.
It does also work 'perfectly' and fast in 'Unicode' in all major European
and CJK languages (and many others): extending the iconv
On Fri, 1 Sep 2006, Hrishikesh Rajpathak wrote:
> Hi,
>
> I am a newbie to R and trying to implement parallelism in R. I am
> currently using R-2.3.1, and Cygwin to run R on Windows xp.
Did you build R under Cygwin (which is not a supported platform), or are
you running a native Windows
?formals
?nargs
depending if you mean the number of formal arguments or the number of
actual arguments.
On Thu, 31 Aug 2006, Lord Tyranus wrote:
> I create a function with R, I want to know how many argument number.
> Thanks in advance.
>
>
--
Brian D. Ripley, [EMAIL PROTEC
On Thu, 31 Aug 2006, David Ziger wrote:
> I cannot generate png files running R2.2.1 on a SunOS 5.9 machine. See
> below. I did a search on this problem and people who have virtually the
> same error have solved the problem by installing or updating libpng .
> We downloaded libpng from Sun but
On Thu, 31 Aug 2006, Jeff D. Hamann wrote:
> I've been trying to get some code to run on a 64-bit FreeBSD machine running
>
> R : Copyright 2006, The R Foundation for Statistical Computing
> Version 2.3.0 (2006-04-24)
> ISBN 3-900051-07-0
>
> and keep getting a core dump with the following resul
On Thu, 31 Aug 2006, Monica Pisica wrote:
> Hi list,
>
> I have a question about 'plot'. I am trying to plot values registered every
> month - or every other month. If i build a data.frame called mydata like
> this (as an example)
>
> jan 3 1 7
> mar 2 4 2
> may 1 3 2
> jul3
901 - 1000 of 5059 matches
Mail list logo
