Paul Lemmens [EMAIL PROTECTED] writes:
Hello!
I'm not grasping why cbind (in the code below) warns that
Warning message:
number of rows of result
is not a multiple of vector length (arg 2) in: cbind(z, p)
when I do
sections - function(length, parts)
{
p - 1:parts
Ronaldo Reis Jr. [EMAIL PROTECTED] writes:
Hi,
I make an analysis and depending of the order of the variables, the
significance change, look.
m1 - glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot)
anova(m1,test=F)
Analysis of Deviance Table
Model: binomial, link: logit
Robin Hankin [EMAIL PROTECTED] writes:
Hello everybody.
I've been experimenting with mapply(). Does anyone else have problems with:
R mapply(rep,times=1:4, MoreArgs=42)
(I get a seg fault).
So do I. It's a bug, but I suppose list(42) is what you need.
--
O__ Peter Dalgaard
Philippe Glaziou [EMAIL PROTECTED] writes:
Rau, Roland [EMAIL PROTECTED] wrote:
yesterday I took the R-1.8.0-source file and compiled it on my
own. As I am using Linux just for a couple of weeks, it was my
first compiling session with ./configure, make,
Everything went fine, except
Tord Snall [EMAIL PROTECTED] writes:
Null deviance: 13.1931 on 269 degrees of freedom
Residual deviance: 9.9168 on 268 degrees of freedom
AIC: 13.917
...
BUT, note the under dispersion. I GUESS it is because I have surveyed a
moss on marked trees at three occations (with two years
Nathan Cooper [EMAIL PROTECTED] writes:
I have an indexing question. I have a data set that looks like this:
b
[[1]]
[1] 22 23 24 25 26
[[2]]
[1] 6 7 8 9 NA
etc. from [[1]] to [[1000]]
Then I need to use the sample function to take two samples from b[[1]]
to b[[1000]] each
Thomas Bock [EMAIL PROTECTED] writes:
Dear list,
I can not understand why the expression in
the subject does not work correct:
dcrn[which(fn == inve[2])]
numeric(0)
inve[2]
[1] 406.7
dcrn[which(fn == 406.7)]
[1] 1.3994e-07 1.3988e-07 1.3953e-07 1.3966e-07 1.3953e-07 1.3968e-07
I've rolled up R-1.8.0.tgz a short while ago. This is a new version
with major changes (see below). Notably, the Macintosh version for OS
X has been substantially improved; the old Carbon interface is no
longer being supported.
Also notice that the underscore will no longer work as an assignment
Spencer Graves [EMAIL PROTECTED] writes:
This seems to me to be a special case of the general problem of
a parameter on a boundary.
Umm, no...
I have this problem with my data. In a GLM, I have 269 zeroes and
only 1 one:
I don't think that necessarily gets you a parameter estimate
Uwe Ligges [EMAIL PROTECTED] writes:
Yair Snir wrote:
Hello,
I have a problem, I can't install the package 'mgu74av2cdf'. I
downloaded the zip file, yet when asked the R console to install it
from a zip file, I got the answer:
Error in file(file, r) : unable to open
Laura Quinn [EMAIL PROTECTED] writes:
I am dealing with a huge matrix in R (20 columns, 54000 rows) and have
lots of missing values within the dataset which are currently displayed as
the value -999.00 I am trying to create a new matrix (or change the
existing one) to display these values as
Rick Bilonick [EMAIL PROTECTED] writes:
I'm still having problems installing rimage - the installation can't
find the fftw headers. As suggested, I installed the fftw rpm (for RH
9 from freshrpms). It installed without any errors or warnings. Yet I
get exactly the same error message - it
Javier Arsuaga [EMAIL PROTECTED] writes:
I am trying to install R in Linux 8.0 and I downloaded
R-1.7.1-1.i386.rpm and did rpm -hiv R-1.7.1-1.i386.rpm
and I am getting the following message:
warning: R-1.7.1-1.i386.rpm: V3 DSA signature: NOKEY, key ID
97d3544e
error: Failed dependencies:
Federico Calboli [EMAIL PROTECTED] writes:
Dear All,
I suspect this is kind of dumb, but when I was under the thrall of the
dark lord (read, using a W2K box), all my work in R files came out as
foo.RData. I moved on to GNU/Linux, and all the old .RData files keep on
working as they used.
Allan Tingey [EMAIL PROTECTED] writes:
Hi,
We are considering using R on a large problem that will require four
64 bit cpu's and 64Gb of ram and some flavor of unix. We are just
wondering if R is going to be happy in this environment. If have
experience with this please let me know.
Karin Gerard Preusting [EMAIL PROTECTED] writes:
Dear all,
I just started using R today. What I need to do is find a text
string in a (large) text file and than copy the some position
related lines (or text) to an other text file. Is this possible with
R?
It's hardly a typical R
Sam McClatchie [EMAIL PROTECTED] writes:
System info:
Red Hat 9.0
R Version 1.7.0
ESS 5.1.21
Emacs 21.2.1
---
Colleagues
At work I've had to migrate from Mandrake 9.1 to Red Hat 9.0 and I'm
reinstalling R. I am having a problem downloading from
R-1.7.1-1.i386.rpm
Pingzhao Hu [EMAIL PROTECTED] writes:
Hi All,
I a C function code, which has been successed in integrating
with Splus on Unix. Now I want to move it on windows.
So I want to integrate this C function with R on windows XP.
I download tools.zip and saved it in c:\bin directory,
download
Thomas W Blackwell [EMAIL PROTECTED] writes:
Ed -
You seem to have encountered a bug. I can reproduce Ed's difficulty
in a completely artificial example in which there are unused levels :
tmp - factor(rep(seq(10), seq(10))) # length(tmp) # [1] 55
ave(seq(50), tmp[-seq(5)])
Laurent Faisnel [EMAIL PROTECTED] writes:
Hi,
I'm trying to compile R-1.7.1 from source (on a RedHat 8.0) instead of
using the binary version, as it has often been advised.
However I don't manage to find a solution to the following error which
occurs during the make procedure :
gcc
Morgan Hough [EMAIL PROTECTED] writes:
I received an error trying to build R-1.7.1-src.rpm on RH9.0.93 (Severn)
and I was wondering if anybody could tell me what is going on with the
manual error.
I'd look further up the log for an indication of whether the PDF
manuals got built at all. You
[EMAIL PROTECTED] writes:
Dear Peter,
I don't know if this is the proper way to ask for help installing R, but if
not, I presume you can pass this on to the appropriate place.
I'm trying to install the latest binary on Redhat 9, and keep getting the same
error message no
[EMAIL PROTECTED] writes:
This is a relatively recent article that is somewhat accessible.
Jensen, D. R., and Solomon, Herbert (1994), Approximations to joint
distributions of definite quadratic forms, Journal of the American
Statistical Association, 89 , 480-486
It has references to
Jason Turner [EMAIL PROTECTED] writes:
On Thu, 2003-09-25 at 05:34, Christoph Lehmann wrote:
Hi
I have calculated lots (1000) of linear models and would like to store
each single result as an element of a 3D matrix or a similar storage:
something like
glm[i][j][k] = lm(...)
Remko Duursma [EMAIL PROTECTED] writes:
Dear R-helpers,
i get some strange results using a linear mixed-effects model (lme), of the type:
lme1 - lme(y ~ x, random=~x|group, ...)
For some datasets, i obtain very small standard deviations of the random effects. I
compared these to
Thomas W Blackwell [EMAIL PROTECTED] writes:
Ruben -
Why not simply save(x, file=new.file.name) ?
See help(save), help(files). The file name must be
quoted, and it must be passed as a named argument to save().
That's not the issue. tclvalue(a) instead of as.character(a) should be
[EMAIL PROTECTED] writes:
Dear R users,
I want to retrieve ... argument values within a function. Here is a small
exmaple:
myfunc - function(x, ...)
{
if (hasArg(ylim)) a - ylim
plot(x, ...)
}
x - rnorm(100)
myfunc(x, ylim=c(-0.5, 0.5))
Error in myfunc(x, ylim =
David Brahm [EMAIL PROTECTED] writes:
Our Sun Ultra-80 running Solaris 2.8 (platform = sparc-sun-solaris2.8) failed
to make R-1.8.0alpha_2003-09-15.tar.gz. The tail end of the output follows.
Note that we have gcc version 2.95.3, as a recent attempt to upgrade failed
miserably. Any ideas?
Rafael A. Irizarry [EMAIL PROTECTED] writes:
for word documents submitted to picky journals i usually use seomthing
like this:
bitmap(plot_1.png,width=6,height=6,res=600,pointsize=12,family=Times)
on my computer this resutls in quality just as good (to mu eye) as with
postscript. you
Martin Maechler [EMAIL PROTECTED] writes:
as.numeric (levels (y)[as.numeric (y)])
NaLi [1] 1 3 1
NaLi which is a bit awkward.
as.numeric(as.character(y)) !
{in some cases you might consider using as.integer() instead of as.numeric()}
Better to convert first and then index:
array chip [EMAIL PROTECTED] writes:
Hi all,
Thanks for all the suggestions. I was able to get the
data frame out by first deleting some other large
objects in the directory, and then changing the data
frame into matrix by as.matrix(), splitting the matrix
into 4 blocks and finally using
Ricardo Pietrobon [EMAIL PROTECTED] writes:
IDdatecost
1 2001-01 200.00
1 2001-01 123.94
1 2001-03 100.23
1 2001-04 150.34
2 2001-03 296.34
2 2002-05 156.36
I would like to obtain the median costs and boxplots
Douglas Bates [EMAIL PROTECTED] writes:
x - c(1:10) #data to be broken up into dummy variables
v - c(3,5,7) #breakpoints
p = 1#drop this column to avoid dummy variable trap
How can I get a matrix y that has the associated dummy variables for
columns?
William Noble [EMAIL PROTECTED] writes:
Hello,
I am trying to do an ANOVA on a microarray data set consisting of
22690 elements. The ANOVA is fine, but when I try to put the data in
a frame in order to exporting it, I get a stack overflow. I have
found documentation on dynamic memory in
Ulrich Leopold [EMAIL PROTECTED] writes:
Dear list,
could someone point me to the right command to subtract 2 columns in a
data.frame. Might be a bit embarrassing question. But I cannot figure
out how to do this simple command in R.
E.g.,
mydata$difference - mydata$x - mydata$y
--
Ulrich Leopold [EMAIL PROTECTED] writes:
E.g.,
mydata$difference - mydata$x - mydata$y
That's what I thought, but I get the following message:
propLSK.STONE.Pox0t30$Pox0t30STONE-propLSK.STONE.Pox0t30$Pox0t30
numeric(0)
Does it mean the resulting vector is empty? If yes, what
Thomas Lumley [EMAIL PROTECTED] writes:
I think it's clear that it will parse as either
(x==y) != z
or
x == (y!=z)
but not which.
The rule is that everything is left associative except assignment and
exponentiation (and IF, for some reason). If in doubt, just remember
that it is the
Richard A. O'Keefe [EMAIL PROTECTED] writes:
Thomas Lumley [EMAIL PROTECTED] wrote:
I would have said that the behavior of
(if (cond) names else dim)(x) - 10
is undefined in the S language, along with things like the order of
evaluation of the apply functions.
Paul Meagher [EMAIL PROTECTED] writes:
I am wanting to construct a probability distribution for height and then,
hopefully, visually and analytically demonstrate that it is normally
distributed.
These are the commands I have developed so far:
fat - read.table(fat.dat, header=TRUE)
mu
Michael Fugate [EMAIL PROTECTED] writes:
## BEGINNING OF CODE ###
# a fake dataset to make the bumps with
nn - 30 # of data points
mm - 7 # number of support sites for x(s)
# create sites s
ss - seq(1,10,length=nn)
# create the data y
e1 -
Thomas Lumley [EMAIL PROTECTED] writes:
I would have said that the behavior of
(if (cond) names else dim)(x) - 10
is undefined in the S language, along with things like the order of
evaluation of the apply functions.
Actually, following Luke's analysis, I think it is pretty squarely an
Barry Rowlingson [EMAIL PROTECTED] writes:
Peter Dalgaard BSA wrote:
BTW: If you ever actually need to do something like that, try
eval(substitute(foo(x)-10,list(foo=as.name(if (cond) names else
dim
please no!!! What's wrong with:
if(cond){
names(x) - 10
} else
rab [EMAIL PROTECTED] writes:
I installed R 1.7.1 on my intel desktop (RH 9) without any problems.
I installed the same binary (R-1.7.1-1.i386.rpm) on my Toshiba
Satellite 2535cds (RH 9). R starts as usual all the way to how to quit
then prints:
Illegal instruction
and exits to the
Spencer Graves [EMAIL PROTECTED] writes:
A common, punishing error for me, with DF being a data
frame, is the following:
if(DF$a = 1) ...
where I intended to write if(DF$a == 1) This error
first replaces column a of DF with the trivial vector 1 (of length
[EMAIL PROTECTED] (Bernhard Bals) writes:
Hi,
I'm just about to install the Precompiled Binary Distribution of R
from CRAN for SuSE Linux 8.2: R-base-1.7.1-1.i386.rpm
My system is SuSE Linux 8.2 Professional Edition on a Pentium 4:
# uname -a
Linux linux 2.4.20-4GB #1 Mon Mar 17
Thomas Petzoldt [EMAIL PROTECTED] writes:
Marc Schwartz wrote:
In general, if you want to leave the existing device open and have a new
device open for a new plot, you simply call the device name that you
want to open (ie. under Linux you would use X11() ) to open a new
plotting device
Monica Palaseanu-Lovejoy [EMAIL PROTECTED] writes:
Hi,
I am a beginner user of R. I have a trivial question I am almost
ashamed I cannot figure it out does not matter how many times I
am reading the help.
I have a table in .txt format, tab delimited. I can read it with
giovanni caggiano [EMAIL PROTECTED] writes:
Dear All,
A couple of questions about the nls package.
1. I'm trying to run a nonlinear least squares
regression but the routine gives me the following
error message:
step factor 0.000488281 reduced below `minFactor' of
0.000976563
even
Dirk Eddelbuettel [EMAIL PROTECTED] writes:
On Tue, Aug 26, 2003 at 03:17:19PM -0400, Liaw, Andy wrote:
Has anyone tried using R on the the AMD Opteron in either 64- or 32-bit
mode? If so, any good/bad experiences, comments, etc? We are considering
getting this hardware, and would like
Simon Fear [EMAIL PROTECTED] writes:
I tried to define a function like:
fnx - function(x, by.vars=Month)
print(by(x, by.vars, summary))
But this doesn't work (does not find x$Month; unlike other functions,
such as
subset(), the INDICES argument to by does not look for variables in
Marc Schwartz [EMAIL PROTECTED] writes:
On Wed, 2003-08-27 at 09:21, [EMAIL PROTECTED] wrote:
Hi guys,
Yesterday I posted my first couple of questions (see bottom of this
message) to this forum and I would like to thank you guys for all the
useful feedback I got. I just would like
Sven Garbade [EMAIL PROTECTED] writes:
Hi all,
suppose I've got a vector y with some data (from a repeated measure
design) observed given the conditions in f1 and f2. I've got a model
with two unknown fix constants a and b which tries to predict y with
respect to the values in f1 and f2.
Spencer Graves [EMAIL PROTECTED] writes:
The Battelle Institute surely should have access to a library with
such popular and prestigious journals as Biometrika and The American
Statistians. If you don't have time for that, you surely should have
money to purchase a copy from, e.g.,
Thomas W Blackwell [EMAIL PROTECTED] writes:
help(Subscript)
Or, princomp or factanal, depending on what was meant...
I have a 12(cols) by 9(rows) matrix X. I need to reduce this matrix so that
it contains 'n' columns (eg. reduce X into a 3 by 9 matrix). What is the
best method to do
Isaac Neuhaus [EMAIL PROTECTED] writes:
We have a loan SGI Altix system runing linux with the IA-64 chip. I
was wondering if anybody has compiled R in such a system.
Seems so...:
http://buildd.debian.org/fetch.php?pkg=r-basever=1.7.1.cvs.20030814-1arch=ia64stamp=1060920704file=logas=raw
--
[EMAIL PROTECTED] writes:
I tried to use substitute in legend as follows:
pval - 0.04
plot(0)
legend(1,0.5,substitute(hat(theta) == p, list(p = pval)))
For some reason the legend is repeated 3 times.
Any suggestions or is this a bug?
It's a bug. The code is looking at
Deepayan Sarkar [EMAIL PROTECTED] writes:
legend(1,0.5,as.expression(substitute(hat(theta) == p, list(p = pval
Just out of curiosity: ?legend says
legend: a vector of text values or an 'expression' of length = 1 to
appear in the legend.
Is an object of mode call
Lun Li [EMAIL PROTECTED] writes:
Dear All,
Does anyone know how to compare two character strings in R? For eample, how
to compare A-1-B with cc-1 in logical ?
Ehhh... Is it this that you want?
A-1-B == cc-1
[1] FALSE
A-1-B = cc-1
[1] TRUE
A-1-B cc-1
[1] FALSE
Lun Li [EMAIL PROTECTED] writes:
Thanks.
But, my question is:
pro.name[1]
V1
1 B-29BT3
temp.name[1]
V1
1 A-1H
if(tem.name[1]==pro.name[1]){cat(ok)}
Error in Ops.factor(left, right) : Level sets of factors are different
How to compare them?
Then they are factors,
Ross Boylan [EMAIL PROTECTED] writes:
Thanks to Thomas Lumley, Spencer Graves, and Steve Sullivan for their
advice on and off list to modify the formula with +0 (which I did; it
worked) or -1 at the end.
...
Also, I noticed that S-Plus but not R has a glim routine that uses
maximum
Murray Jorgensen [EMAIL PROTECTED] writes:
I want to interlace two vectors. This I can do:
x - 1:4
z - x+0.5
as.vector(t(cbind(x,z)))
[1] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
but this seems rather inelegant. Any suggestions?
Well, there's as.vector(rbind(x,z)) at least... I don't
Laurent Faisnel [EMAIL PROTECTED] writes:
Hi all,
Trying to install R-1.7.1 on a RedHat 8.0 platform, I have a few problems.
R gets installed without default packages (but base and ctest) : make
script fails at the end of the procedure (configure is made
successfully). This is the
Tito de Morais Luis [EMAIL PROTECTED] writes:
Dear listers,
The following command (derived from the example in the ?stars help page)
works :
data(mtcars)
stars(mtcars[, 1:7])
But the following gives an error:
stars(mtcars[1, 1:7])
Error in s.y[i, ] : incorrect number of dimensions
Angel [EMAIL PROTECTED] writes:
Why do
x-b%*%ginv(A)
and
x-solve(A,b)
give different results?. It seems that I am missing some basic feature of
matrix indexing.
e.g.:
A-matrix(c(0,-4,4,0),nrow=2,ncol=2)
b-c(-16,0)
x-b%*%ginv(A);x
x-solve(A,b);x
[ginv() is from MASS, please remember
Uwe Ligges [EMAIL PROTECTED] writes:
Dennis Fisher wrote:
In labeling axes, I want to combine symbols and text/superscripts.
Examples include:
m2 (m, followed by a superscripted 2)
*
µg (micrograms)
How can I accomplish this in R?
See
Hedderik van Rijn [EMAIL PROTECTED] writes:
Some more information about my read.spss issues: on a cleanly
installed Debian Woody machine:
~/TMP % R
R : Copyright 2003, The R Development Core Team
Version 1.7.1 (2003-06-16)
[...]
library(foreign); x - read.spss(dataDef.sav)
[EMAIL PROTECTED] writes:
Hello everybody,
could anybody give me a hint, who I can use rbind on a list of data.frames, please?
I have a list with a large number of data.frames of the same structure, like:
LIST - list(X1=data.frame(a=1,b=2), X2=data.frame(a=3,b=4),
Simon Fear [EMAIL PROTECTED] writes:
Peter, thanks for quick reply, but I can't drop it:
system(dir /b, intern=T, show.output.on.console=T)
Error in system(dir /b, intern = T, show.output.on.console = T) :
dir not found
Oh, dir is a built-in, right?... Better wait for the real
Martin Maechler [EMAIL PROTECTED] writes:
jane == jane murray [EMAIL PROTECTED]
on Wed, 06 Aug 2003 09:16:48 +0100 writes:
jane Hi Can anyone help with the technique of obtaining
jane leverages from a conditional logistic regression
jane model? The code lm.influence does
Simon Fear [EMAIL PROTECTED] writes:
Prof Brian Ripley writes:
If you insist on using a long-obselete OS, please help us continue ot
support it by supplying patches and workarounds. I think we (like
Microsoft) are going to have to think seriously about withdrawing
support
from
Simon Fear [EMAIL PROTECTED] writes:
The following code works, to gsub single quotes to double quotes:
line - gsub(', '', line)
(that's a single quote within doubles then a double within singles if
your
viewer's font is not good).
But The R Language Manual tells me that
Quotes and
Dibakar Ray [EMAIL PROTECTED] writes:
I am very new to R. I was trying to load some publicly available Expression
data in to R.
I used the following commands
mydata-read.table(dataALLAMLtrain.txt, header=TRUE, sep
=\t,row.names=NULL)
It reads data without any error
Now if I use
Jim Lemon [EMAIL PROTECTED] writes:
Johnathan Williams wrote:
Would anyone be so kind as to write a routine to time
mouse button presses in R to the nearest millisecond?
If R had a timer of this kind and a few basic screen
handling routines (to write characters or graphics of
Simon Fear [EMAIL PROTECTED] writes:
system(COMMAND.COM /c dir /b, intern=T, show.output.on.console=T)
character(0)
print(system(COMMAND.COM /c dir /b, intern=F,
show.output.on.console=T))
[1] 0
In both cases here a DOS window opens and lists a couple of hundred
files,
before giving
Simon Fear [EMAIL PROTECTED] writes:
Peter is right, but there is a point here: it would be nice if
attach(my.dframe) would do nothing - or at least warn - if my.dframe was
already in the search list. Like library(). Attaching twice is almost
bound
to be an error?
The problem is, how can
Simon Fear [EMAIL PROTECTED] writes:
It's one of the many situations in which I would very much like to
get a warning or error message, pointing out to me that I had
absolutely no idea what I was doing.
Surely that's what warnings are for? For those of us who wonder why
our code doesn't do
Andreas Eckner [EMAIL PROTECTED] writes:
Hi,
I just encountered a problem in R that may easily be fixed: If one uses
attach for a data.frame e.g. 1 times and forgets detach, then R gets
incredibly slow (less then 10% of the original speed).
R also gets incredibly slow if you create
Thomas Lumley [EMAIL PROTECTED] writes:
On Wed, 30 Jul 2003, Damon Wischik wrote:
Why does
curve(function(y) y^2, from=0,to=1)
not work, whereas
myf - function(y) y^2
curve(myf, from=0,to=1)
work?
Because someone was trying to be too clever in writing curve(). The first
Prof Brian Ripley [EMAIL PROTECTED] writes:
On Tue, 29 Jul 2003, p.b.pynsent wrote:
What I had hoped for was what I get with srt=0, rotated as whole so it
looked the same but was just rotated to the specified angle.
But you got what is documented, and my example addessed that.
The docs
Steve Sullivan [EMAIL PROTECTED] writes:
neighbor.id - vector(length=length(D1$lat))
dist.geo - vector(length=length(D2$lat))
for(i in 1:length(neighbor.id)){
for(j in 1:length(dist.geo)){
dist.geo[j] - D1$lat[i]-D2$lat[j]}
# Yes, I know that isn't the right formula, this is just a test
Buchsbaum, Brad (NIH/NIMH) [EMAIL PROTECTED] writes:
Hi,
Suppose I have a multidimensional array:
tmp - array(1:8, c(2,2,2))
is there a function out there that, given a one-dimensional array index,
will
return the separate indices for each array dimension?
for instance, tmp[8] is
John Chambers [EMAIL PROTECTED] writes:
- not evaluating the left side of the operator. This is no problem in
principle, since generic functions can have a signature (see
?setGeneric); if the signature omits the first argument, methods can be
dispatched without evaluating that argument.
Peter Dalgaard BSA [EMAIL PROTECTED] writes:
Prof Brian Ripley [EMAIL PROTECTED] writes:
What are you trying to do with this? Assignment (-) is not a function,
and the language grammar does not convert a - b into -(a, b) (as it
would with the binary operator functions). You could call
Prof Brian Ripley [EMAIL PROTECTED] writes:
x y
0 28
0.0321
0.1 11
0.3 15
1 5
3 4
10 1
30 0
100 0
Where X is dose and Y is response.
the relation is linear for log(response) =
Hotz, T. [EMAIL PROTECTED] writes:
Dear Tobias,
The trick is Programming on the Language, see e.g. the R Language Manual.
Construct the expression you want, and have it explicitly parsed and evaluated.
toy - function(b=.95){
toyframe - eval(parse(text=paste(data.frame(lion, b, = c(1,
Spencer Graves [EMAIL PROTECTED] writes:
The problem with nls is that it is NOT maximum likelihood for the
Poisson distribution. For the Poisson, the standard deviation is the
square root of the mean, while nls assumes constant standard
deviation. That's why I stayed with glm. The answers
Bear F. Braumoeller [EMAIL PROTECTED] writes:
I'm starting R with
xterm -sb -rightbar -sl 1000 -bg black -fg blue -title R -e
/usr/local/bin/R
-- but it also happens if I just start a vanilla terminal and type R.
As to the other questions,
Sys.getenv(TERM)
TERM
xterm
Bear F. Braumoeller [EMAIL PROTECTED] writes:
I recently reformatted my TiBook and reinstalled everything from the
system up. After installing Apple's version of XWindows and the
XWindows version of R 1.7.0, I ran the two and got a strange warning:
-
R : Copyright 2003, The R
Douglas Bates [EMAIL PROTECTED] writes:
I happened to look up Frank Harrell's book Regression Modeling
Strategies on Amazon.com today. I was surprised to see that in
addition to the typical links to related books they had sponsored
links to sites about How to become a model, Try out for
Spencer Graves [EMAIL PROTECTED] writes:
The reference to Venerable and Ripley reminds me that in the late
1940s, George Box worked for Imperial Chemicals Industries, which
appeared in at least one brochure I saw as Empirical Chemicals
Industries. George invented response surface methods
Adi Humbert [EMAIL PROTECTED] writes:
Dear list,
I'm migrating a project from Matlab to R, and I'm
facing a relatively complicated problem for nls. My
objective function is below:
objFun - function(yEx,xEx,tEx,gamma,theta,kappa){
yTh - pdfDY(xEx,tEx,gamma,theta,kappa)
J. P. Morgan [EMAIL PROTECTED] writes:
Does anyone know how to execute the following sort problem in R? Matrix X
has positive integer entries, and each column has been sorted in ascending
order. The problem is now to order the columns lexicographically. For
instance if the matrix X is
[Argh. Forgot to cc: the list again...]
Ted Harding [EMAIL PROTECTED] writes:
mu[M[,1:3]]
should do what you want, I think. See chapter 1 of S Poetry.
Patrick Burns
Thanks Patrick (and honoured to hear from the author of S-Poetry!).
This doesn't do what I want (I was probably
[Oops. Accidentally sent only to D.Trainor first time around, even
though it was intended mainly for M.Kondrin and the list]
Douglas Trainor [EMAIL PROTECTED] writes:
Uwe Ligges wrote:
M.Kondrin wrote:
Hello!
Some CS-guys (the type who knows what Church formalism is) keep
asking me
Murad Nayal [EMAIL PROTECTED] writes:
[Sensible stuff, except:]
what sort of garbage collector it uses?
No garbage collector. uses reference counting to discard objects that
are no longer needed.
No, R uses garbage collection with a generational scheme which looks
more frequently at
Thomas W Blackwell [EMAIL PROTECTED] writes:
Runge -
I haven't tried it out, but half a guess says that R might not like
using underscore in a variable name. Please try exactly the same
command without quotes and without the underscore:
fraktil.df -
Ross Boylan [EMAIL PROTECTED] writes:
I want to maximize a conditional likelihood function that is basically
logistic conditional on the number of successes within strata. What
would be a good starting place for this? A complication is that the
denominator includes a term that is the sum
Shi, Tao [EMAIL PROTECTED] writes:
Hi, Ted:
I guess this problem is platform-dependent. I just tied it on a R
1.6.1 runing on Win2K, it gave me two different p values. But when I
tried it on R1.7.0 on a Linux Server, I got the similar result as
you did. I have filed a bug-report as Peter
Jesper Runge Madsen [EMAIL PROTECTED] writes:
Dear R helpers
I am trying to convert a list into a data frame but when I try, I get a
stack overflow error (Error: protect(): stack overflow). My list contains
about 17000 rows and looks like shown at the bottom. The reason that I
want to
Shi, Tao [EMAIL PROTECTED] writes:
Hi, Ted and Dennis:
Thanks for your speedy replies! I don't think this happens just randomly, rather,
I'm thinking it may be due to the way chisq.test function handles simulation. Here
shows why: (Ted, I think there is an error in your code, tx should
1 - 100 of 203 matches
Mail list logo
