On Tue, 9 Sep 2008, [EMAIL PROTECTED] wrote:
As a sys admin, how do I install packages so that there is one common
library for all users of the MS Windows computer, instead of the
default individual location for each user?
I've done this for Linux.
The same way. You set R_LIBS_SITE in
This isn't accurate. You are talking about link functions *known by name*.
link: a specification for the model link function. This can be a
name/expression, a literal character string, a length-one
character vector or an object of class 'link-glm' (provided
it
Rajdeep Das wrote:
Hi,
I am trying to calculate probability density of normal inverse gaussian
distribution. I am using dnig function of fBasics package. However, I am
getting following result. The density at x = 0.003042866 is:
dnig(x= 0.003042866, alpha=5.184868, beta= 0.11841, delta=
[EMAIL PROTECTED] wrote:
Does anyone know why I get the following error when trying tsdiag?
Error in UseMethod(tsdiag) : no applicable method for tsdiag
I am invoking it as: tsdiag(mar).
tsdiag is a generic function that tries to dispatch for mar, but there
is no appropriate method
pgseye wrote:
After doing a PCA using princomp, how do you view how much each component
contributes to variance in the dataset. I'm still quite new to the theory of
PCA - I have a little idea about eigenvectors and eigenvalues (these
determine the variance explained?). Are the eigenvalues
Andreas Wittmann wrote:
Dear R useRs,
i have the following code to compute values needed for a contour plot
myContour - function(a, b, plist, veca, vecb, dim)
{
tmpb - seq(0.5 * b, 1.5 * b, length=dim)
tmpa - seq(0.5 * a,
Dear R-users,
I have a couple of questions about the relsurv package:
1) when I try to run the example:
fit -
rsmul(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365.24,sex=sex,year=year),ratetable=slopop,data=rdata)
with the datasets in the package (rdata and slopop) it gives me
Many thanks Luke, much appreciated,
Tolga
Luke Tierney [EMAIL PROTECTED]
09/09/2008 21:21
To
[EMAIL PROTECTED]
cc
Prof Brian Ripley [EMAIL PROTECTED], r-help r-help@r-project.org
Subject
Re: [R] Information on the number of CPU's
The wmic command line utility can also be used to query
Raj has mixed up the concept of density function with CDF
(Cumulative Distribution Function) - the former one doesn't need to be
smaller than 1, while the latter one is certainly no greater than 1.
Yihui
On Wed, Sep 10, 2008 at 3:55 PM, Uwe Ligges
[EMAIL PROTECTED] wrote:
Rajdeep Das wrote:
Dear R community
The following code gives me the most repeated sequence of values. i.e.
x=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,2,2,2,3,3,0,0,0,0,0,0,0,0,0,0,1,1,1,2,2,3,3,3,4,4,4,0,0,0,0,0,0,0,1,1,1,2,2,2,3,3,3,4,4,4,
For the second question, try something like:
strptime('12/4/1996','%M/%d/%Y')
[1] 1996-09-04 00:12:00
# only date?
as.Date(strptime('12/4/1996','%M/%d/%Y'))
[1] 1996-09-04
Yihui
On Wed, Sep 10, 2008 at 4:22 PM, Giulia Barbati [EMAIL PROTECTED] wrote:
Dear R-users,
I have a couple of
DScottNZ == David Scott [EMAIL PROTECTED]
on Wed, 10 Sep 2008 11:29:31 +1200 (NZST) writes:
DScottNZ On Tue, 9 Sep 2008, Edna Bell wrote:
Dear R Gurus:
I want to look at the code for the t.test function. I did the following:
t.test
function (x, ...)
Nic,
I'd buy a Mac Power server.
Not that's much faster, but while one of the cores toils away
at R you can play games on the others :-)-O
el
on 9/9/08 8:54 PM Henrik Bengtsson said the following:
On Tue, Sep 9, 2008 at 6:31 AM, Nic Larson [EMAIL PROTECTED] wrote:
Need to buy fast computer
Dear Mr./Mrs.,
I would like to inform you about the new course Pharmacokinetic and pharmacodynamic
modeling and simulation by Dr. Jan Freijer in Amsterdam, November 28th.
The course is meant for users of R or S-PLUS in the bio-pharmaceutical sciences
who would like to use R for clinical trial
Dear R Users,
Apologies for a somewhat basic cluster question. I am trying to come to
grips with how variables are instantiated within the environment of each
node in a cluster. When I run the following code, I expect to see the
variable x in the environment of each node, but fail to do so.
On Wed, 10 Sep 2008, Yihui Xie wrote:
For the second question, try something like:
strptime('12/4/1996','%M/%d/%Y')
[1] 1996-09-04 00:12:00
# only date?
as.Date(strptime('12/4/1996','%M/%d/%Y'))
[1] 1996-09-04
Why not
as.Date('12/4/1996', format = '%m/%d/%Y')
[1] 1996-12-04
? (And
Hi Tolga,
I am not sure why ls() is showing only args and fun.
But x is at the nodes. Try this:
clusterEvalQ(cl,x-x+3)
[[1]]
[1] 4
[[2]]
[1] 4
clusterCall(cl,function()x+5)
[[1]]
[1] 9
[[2]]
[1] 9
Best
Markus
[EMAIL PROTECTED] schrieb:
Dear R Users,
Apologies for a somewhat basic
Whats the R equivalent for Proc logistic in SAS ? Is there a stepwise
method there ?
How to create scoring models in R , for larger datasets (200 mb), Is
there a way to compress and use datasets (like options compress=yes;)
Ajay
On Wed, Sep 10, 2008 at 11:12 AM, Peter Dalgaard
[EMAIL PROTECTED]
Dear R community
Thanks. The problem has been solved by just using return(rows)
Muhammad Azam
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hello,
[...]
--
Message: 41
Date: Tue, 9 Sep 2008 9:44:34 -0700
From: [EMAIL PROTECTED]
Subject: [R] Help with 'spectrum'
To: r-help@r-project.org
Message-ID: [EMAIL PROTECTED]
Content-Type: text/plain; charset=utf-8
For the command 'spectrum' I read:
The
Ajay ohri wrote:
Whats the R equivalent for Proc logistic in SAS ?
glm with the appropriate family (binomial) and link, I guess.
There is a book 'R for SAS and SPSS users' forthcoming
http://www.springer.com/statistics/computational/book/978-0-387-09417-5
Is there a stepwise
method there
It is a mixture of both. The data is so notorious excel cant format properly.
Therefore I thought whether R can do something otherwise I have to do manually.
--- On Tue, 9/9/08, Dr Eberhard W Lisse [EMAIL PROTECTED] wrote:
From: Dr Eberhard W Lisse [EMAIL PROTECTED]
Subject: Re: [R]
Hi Bert,
## returns TRUE only if the distributions are the same ...
Thanks for the elegant code. Problem is, the result is elusive, constantly
slipping from view and then rolling into the dip. Perhaps it's broken in
this version of R (2.7.2.Patched). A fix would therefore be much
appreciated,
I have following
1975 01 7711.16
Here I need to identify where the space is there and then concatenate rest of
the digits without space, i.e. I want to have 1975017711.16. Is there any R
function?
Regards,
__
R-help@r-project.org mailing list
Why don't you use 'try' in this fashion:
f.error - function(x) if (x == 1) stop('error')
value - try(f.error(1))
Error in f.error(1) : error
if (inherits(value, 'try-error')) cat(Got this error:, value) else print
(no error)
Got this error: Error in f.error(1) : error
value -
Is this what you want:
gsub( , , 1975 01 7711.16)
[1] 1975017711.16
On Wed, Sep 10, 2008 at 8:07 AM, Megh Dal [EMAIL PROTECTED] wrote:
I have following
1975 01 7711.16
Here I need to identify where the space is there and then concatenate rest
of the digits without space, i.e. I want
Many thanks. Yes, weird, somehow ls does not show the variables in the
environment, even though they are clearly there. It would be nice to
understand why that is the case, just to be clear about what is going on.
Thanks Markus,
Tolga
Markus Schmidberger [EMAIL PROTECTED]
10/09/2008 10:45
on 09/10/2008 07:07 AM Megh Dal wrote:
I have following
1975 01 7711.16
Here I need to identify where the space is there and then
concatenate rest of the digits without space, i.e. I want to have
1975017711.16. Is there any R function?
See ?gsub:
# Strip out blanks
gsub( , , 1975 01
Hello. I've being having problems to create directories in a windows
server environment . It seems that the recursive argument is not working
properly on the intranet, as it does in a local path. For example:
dir.create(server/directory1/directory2) ,
this works fine, and creates the
hello,
thanks a lot for your help.
@ Stephen: In my opinion the proscrutes rotation is more used to
compare the
alignment of different ordination methods. But maybe I'm wrong.
@ Bill: Thanks for the comprehensive description. That was the
information I was
looking for.
Greetings,
Bernd
Thanks everyone,
Paul
--
View this message in context:
http://www.nabble.com/PCA-and---variance-explained-tp19388970p19410675.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Try this:
str - 1975 01 7711.16
# Find the position of the space in string
gregexpr([[:space:]], str)
# Remove the space
gsub([[:space:]], , str)
On Wed, Sep 10, 2008 at 9:07 AM, Megh Dal [EMAIL PROTECTED] wrote:
I have following
1975 01 7711.16
Here I need to identify where the space
Hi Kate,
Things that can reasonably be combined will appear in the combined
randomForest object. Those that cannot are set to NULL. I believe
variable importance is kept, and possibly (pseudo) r-squared as well.
Best,
Andy
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL
Use ksmooth(x, y, kernel=normal)
Andy
From: dxc13
useR's,
Does anyone know what function or package allows one to
specify use of the
Gaussian kernel to smooth a data set, or is this just the
default for most
functions?
Thanks,
dxc13
--
View this message in context:
ls() looks in its calling environment, which in this case would be the
internals of the snow mechnism. If you want ls of the global
environment use one of
clusterEvalQ(cl, ls())
clusterCall(cl, ls, .GlobalEnv)
luke
On Wed, 10 Sep 2008, [EMAIL PROTECTED] wrote:
Many thanks. Yes,
Dear R-gurus,
I know that it is not a R problem, but as need
to implement some solution for this on R, I would
like to know if someone help me on how to isolate
X on the expression below.
K=(sen A * sen B ) + ( cos A * cos B * cos ( X - Y ) )
Any help are welcome
mitinho astronauta
brazil
Hi,
We released individual birds into a room with 2 trees. We counted the number
of visits to each of the 2 tree. One of the trees is always a control tree
and the other tree is either treatment 1, treatment 2 or treatment3 or
treatment 4.
Ind Treat ContrTree ExpTree Total visits
1
So, am I correct that each datum is either of the form mm/dd/yy or
of the form dd/mm/? If that is correct, then the following
should work, and takes care of converting 99 to 1999 instead of 2099:
dates - c(06/15/07,04/09/99,20/03/2008)
short - grep(\\d\\d/\\d\\d/\\d\\d$, dates,
Hi all,
I try to write a simple function in a script. The script is as follows
yo-function(Xdata)
{
n-length(Xdata[,1])
Lgm-nls(formula=LgmFormula,
data=Xdata,
start=list(a=1500,b=0.1),weights=Xdata$Qe)
return(Lgm)
}
After the execution of the script, when I call the function yo on
Dear All
I have a question about the hypergeomteric distribution.
Example 1: I have a universe of 6187 objects, and 164 have a particular
attribute, therefore 6187-164 do not have that attribute. I sample 249
of those objects, and find that 19 have that attribute. I get a p-value
here (looking
The partial dependence at X1=x is simply the average of predicted
response at all data points, holding X1=x. No additional
`normalization' is done.
That web link says the image does not exist.
Andy
From: Christopher Schwalm
Hello ALL,
I'm unsure how to interpret partial dependence plots
Yahoo and OandA provide only daily data, but the history goes back quite a ways:
msft - get.hist.quote(instrument=MSFT, start=1986-03-31,
end=2008-09-10, quote=c(O,H,L,C,A,V), provider=yahoo,
retclass=zoo)
NROW(msft)
[1] 5663
That's 22 years of data.
Not sure what you were expecting.
David
Its looking in environment(LgmFormula) for Xdata in order to evaluate
Xdata$Qe in the weights argument so just specify weights = Qe.
Alternately insert the following as the beginning of yo.
It creates a new local LgmFormula but with its environment set
to the current environment:
Take a look again at help(getVarCov).
Mark Lyman, Statistician
ATK Launch Systems
[EMAIL PROTECTED]
From: huang min [mailto:[EMAIL PROTECTED]
Sent: Tuesday, September 09, 2008 10:02 PM
To: r-help@r-project.org
Cc: [EMAIL PROTECTED]
Subject: Re: extract variance components
Hi,
Try Ryacas:
library(Ryacas)
A - Sym(A); B - Sym(B); X - Sym(X); Y - Sym(Y); K - Sym(K)
Solve(K==(sin(A) * sin(B) ) + ( cos(A) * cos(B) * cos ( X - Y ) ), X)
[1] Starting Yacas!
expression(list(X == Y + acos((K - sin(A) * sin(B))/(cos(A) *
cos(B))), X == Y - acos((K - sin(A) *
Benoit Boulinguiez benoit.boulinguiez at ensc-rennes.fr writes:
yo-function(Xdata)
{
n-length(Xdata[,1])
Lgm-nls(formula=LgmFormula,
data=Xdata,
start=list(a=1500,b=0.1),weights=Xdata$Qe)
return(Lgm)
}
After the execution of the script, when I call the function yo on data
called
On 10 Sep 2008, at 15:19, michael watson (IAH-C) wrote:
Example 1: I have a universe of 6187 objects, and 164 have a
particular
attribute, therefore 6187-164 do not have that attribute. I sample
249
of those objects, and find that 19 have that attribute. I get a p-
value
here (looking
ok, the model does run now!
but, dont i need the method=ML when i want to compare this model with a
reduced model using anova(model1, model2)? The R-Book tells me that REML is
not good for that (p.635)
So, besides that... how do i now do a sort of posthoc test to see 1)
estimates of all
Hi,
I have been dealing with some problems finding a fast way of getting to
know in what dimension a specific value is located out of an array (like
the 'which' function for a vector returns its position).
Unable to find anything about this on the internet I wrote a function
myself.
Could you
Is this what you were looking for:
a - array(100:1,dim=c(5,4,5))
which(a == 7)
[1] 94
which(a == 7, arr.ind=TRUE)
dim1 dim2 dim3
[1,]435
On Wed, Sep 10, 2008 at 10:22 AM, Hintzen, Niels [EMAIL PROTECTED] wrote:
Hi,
I have been dealing with some problems finding a fast
Thanks for your response. However, I need more guidance.
1) I searched the archives of this list back to June 1. I found no
discussion on this subject as I examined all the subject titles for the
string [Ww]in.
2) I tried (what I understand) to be you suggestion.
I created etc/Renviron.site:
The graphic design elements in ggplot2 are so well done for color
graphics. Thank you! For BW publication graphs, there are a few
options I would like to be able to tweak, but can't seem to find the
grob handles necessary OR the right plotting strategy. Here's one
problem I've been
Hi all,
I am trying to use plm to estimate coefficients in a model consisting of a
system of equations. So far I used mle2 from the package bbmle, but now I
need to test for autocorrelation and mle2 does not provide for the necessary
tests. mle2 needs a function as input that might as well
Prof Brian Ripley wrote:
This isn't accurate. You are talking about link functions *known by name*.
link: a specification for the model link function. This can be a
name/expression, a literal character string, a length-one
character vector or an object of class
Hi Adam,
The next version of ggplot2 supports the build-in R plotting symbols
that have different fills and borders, so you can do something like:
geom_point(aes(colour=TRT), fill=red, colour=black, shape=21)
Otherwise in the current version you can do:
ggplot(data.frame(x = runif(20), y =
Many thanks, that is very helpful.
Tolga
Luke Tierney [EMAIL PROTECTED]
10/09/2008 14:05
To
[EMAIL PROTECTED]
cc
[EMAIL PROTECTED], r-help@r-project.org
Subject
Re: [R] cluster/snow question
ls() looks in its calling environment, which in this case would be the
internals of the snow
This sounds vaguely familiar to me. This could be reason for why I
wrote mkdirs() of the R.utils package. mkdirs() does it's own
(manual) recursive creation of directories. See if that does it for
you.
/Henrik
On Wed, Sep 3, 2008 at 1:14 AM, Francisco Jose Sastre [EMAIL PROTECTED] wrote:
Hi R users, I have a question which is some what related to this thread.
Here is the setup:
z is a list of lists, with no names.
zz is same is z, but with names.
with names on zz, I can crawl on the data with lapply, but I don't know how do
the same on z (without names). Can someone help?
###
Thanks Jim and S Ellison for your help
This should do what you want.
#--x - read.table('clipboard', header=TRUE,
as.is=TRUE)
# convert dates
x$date - as.POSIXct(strptime(x$SampleDate,
%m/%d/%Y))
# put ForkLength into bins
x$bins - cut(x$ForkLength, breaks=c(32, 34, 37, 40),
Dear R community
I have stored the results of arrays in a list consist of J-components (say 200
components). Each component containing same no of columns but may be different
no of rows. e.g
[[1]]
[,1] [,2] [,3] [,4] [,5]
[1,]40000
[2,]43400
Is this what you want:
lapply(z, '[[', 1)
[[1]]
[1] 1 2 3
[[2]]
[1] 11 12 13
lapply(z, '[[', 2)
[[1]]
[1] 4 5 6
[[2]]
[1] 14 15 16
On Wed, Sep 10, 2008 at 11:43 AM, Ling, Gary (Electronic Trading)
[EMAIL PROTECTED] wrote:
Hi R users, I have a question which is some what related to this
If want you want is the summary from all of them, then 'rbind' the
data together into one matrix and analyze it:
totalMat - do.call(rbind, listOfMatrices)
On Wed, Sep 10, 2008 at 11:49 AM, Muhammad Azam [EMAIL PROTECTED] wrote:
Dear R community
I have stored the results of arrays in a list
YES! many many thanks, Jim.
I know it'd be simple, yet I can't believe it's that simple. =)
-gary
-Original Message-
From: jim holtman [mailto:[EMAIL PROTECTED]
Sent: Wednesday, September 10, 2008 11:55 AM
To: Ling, Gary (Electronic Trading)
Cc: Peter Dalgaard; r-help@r-project.org
If you are trying to test the hypotheses:
H0: sigma^2 == sigma[0]^2
Ha: sigma^2 != sigma[0]^2
i.e. is the variance of the distribution from which the data came different
from the hypothesized value?
Then the normal theory test can be done fairly simply (I don't know of any
prepackaged
hi: you should probably send below to R-Sig-Finance because there are
some econometrics people over there who could also possibly give you
a good answer and may not see this email ? Also, there's package called
mar ( I think that's the name ) that may do what you want ?
Finally, I don't
Are you really going to look at all the rows?
Probably better options for what you want are to look at functions like 'head',
'tail', 'subset', and 'View' that will let you look at specific parts of the
data rather than printing out much more information than you are likely to care
about, then
Hi everyone,
I am performing the time series regression analysis on a series of data sets. A
few data sets followed an ARMA(1,1) process. However, they all had a same value
of moving average MA coefficients = -1, constantly, from output of function
arima .
Example:
arima(residuals,
hi,
I have list of length 5453 and vector of length 14318.I need to compare the
vector with the list and return the list name if matched.I am thinking of
using an lapply but how to retrive the listname is wat i am puzzled abt.
kindly let me know how to go abt it.
Ramya
--
View this message in
Dear R Gurus:
I would like to set the number of commands in the history file as 50
rather than 25.
I would do this in the .First function.
How do I set this, please? I tried
history(max=50)
and it prints what it has.
I just want to set the size in .First.
Thanks in advance,
Sincerely,
Edna
If I understand, try this:
names(l[l %in% v])
where l is your list and v is your vector.
On Wed, Sep 10, 2008 at 1:38 PM, Rajasekaramya [EMAIL PROTECTED]wrote:
hi,
I have list of length 5453 and vector of length 14318.I need to compare the
vector with the list and return the list name if
Others showed you how to find the mode in a dataset, I just want to point out
that if your data is from a continuous distribution (or near continuous), then
the mode of the data is more likely to be the result of a quirk of rounding
than representative of anything useful. If the data is
Hi,
that's the information about R and R.matlab:
*R version 2.6.2 (2008-02-08)
x86_64-pc-linux-gnu
locale:
C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] R.matlab_1.2.3R.oo_1.4.5R.methodsS3_1.0.3
On Wed, 10 Sep 2008, Edna Bell wrote:
Dear R Gurus:
I would like to set the number of commands in the history file as 50
rather than 25.
It is neither of those. That is the default value of max.show, not of the
size of the file.
I would do this in the .First function.
How do I set
Hi all,
I have run a pca and want to plot the samples multivariate space using the
sample numbers. In addition, I would like to color code the sample numbers
by group, but can't find a way to color-code the sample numbers by group.
So, I have to do it in two separate steps:
plot(bulbil.pca,
Stephen and Jorge,
Perhaps a simpler solution is to use the which function
test.data - Harman74.cor$cov#a test data set
td - test.data * lower.tri(test.data) #this will examine only the
lower off diagonal elements
td.1 - which(abs(td).6,arr.ind=TRUE) # the critical pairs
td.2 -
There are at least 2 more packages which can handle the conversion.
xlsReadWrite (only available on Windows) allows you to read Excel
files directly in R. At the moment xls, but not xlsx and xlsm, AFAIK.
RExcel (available via RExcelInstaller on CRAN)
is an Excel addin which allows you to
I am using R-2.7.2 on Windows XP.
Sorry for the problems.
Sincerely,
Edna
On Wed, Sep 10, 2008 at 11:54 AM, Prof Brian Ripley
[EMAIL PROTECTED] wrote:
On Wed, 10 Sep 2008, Edna Bell wrote:
Dear R Gurus:
I would like to set the number of commands in the history file as 50
rather than 25.
Hi,
I've managed to reproduce the error with R.matlab v1.2.2 and Matlab
v7.4.0.287 (R2007a).
There seems to be a bug introduced in R.matlab v1.2.2 (2008-07-12). In
this version we tried to optimize the speed of writeMat(), and
something failed (although all tests show the generated files loads
I have following dataset:
res
[,1] [,2] [,3]
[1,] 19464 1.27
[2,] 19465 1.27
[3,] 19466 1.27
[4,] 19467 1.27
[5,] 19468 1.52
[6,] 19469 1.52
[7,] 1946 10 1.52
[8,] 1946 11 1.52
[9,] 1946 12 1.62
[10,] 19471 1.62
[11,] 19472 1.62
[12,] 1947
The main difference is that read.xls and xls2csv in the
gdata package use a perl program which works on all
platforms whereas the others require Windows but may
have advantages there.
On Wed, Sep 10, 2008 at 1:40 PM, Erich Neuwirth
[EMAIL PROTECTED] wrote:
There are at least 2 more packages
I can't reproduce that. When I run it I get no warning. Check
your versions:
Lines - 19464 1.27
+ 19465 1.27
+ 19466 1.27
+ 19467 1.27
+ 19468 1.52
+ 19469 1.52
+ 1946 10 1.52
+ 1946 11 1.52
+ 1946 12 1.62
+ 19471 1.62
+ 19472 1.62
+ 19473
Hi all,
I realized what I was doing wrong and aplogize about posting before thinkg
about it for long enough. If I try to plot the pca results directly
(bulbil.pca), I can't do it. But if I extract the axes 1 and 2 and plot
directly, I can add col=as.numeric(group.bulbil[,3]) to the text()
On Wed, 2008-09-10 at 14:24 +0200, Bernd Panassiti wrote:
hello,
thanks a lot for your help.
@ Stephen: In my opinion the proscrutes rotation is more used to
compare the
alignment of different ordination methods. But maybe I'm wrong.
@ Bill: Thanks for the comprehensive description.
RFTW l.temarvelde at nioo.knaw.nl writes:
ok, the model does run now!
but, dont i need the method=ML when i want to compare this model with a
reduced model using anova(model1, model2)? The R-Book tells me that REML is
not good for that (p.635)
maximum likelihood (ML) is the default
Hi,
I've been using R for many years and have always tried to keep my R
version up to date, and when I switched from 2.7.1 to 2.7.2 I'm
suddenly having trouble printing from the graphics device. I have
Windows XP, and installed both versions of R from the binaries. If I
start 2.7.2 and
Hi to all,
there are some binary .dmg packages to install easly rpy under mac os x?
I'm tryin' to install Rpy under Mac Os X but I receive the following
message:
host46-209-dynamic:rpy usernone$ python setup.py install
running install
running build
running build_py
running build_ext
Suppose x and y are numeric vectors of the same length.
plot(x,y) #scatterplot
lmObj1 - lm(y~x) # best fit line
abline(lmObj1) # good
lmObj2 - lm(x~y) #get best fit but with axes interchanged
abline(lmObj2) # not what I want. I want the correct line, drawn on the same
graph, but with
Hi David,
I don´t know if I got what you are looking for.
But see the code below.
x-1:100
y-x+(runif(100)*x)
plot(y~x)
mymod-glm(y~x)
my.coefs-coef(mymod)
my.coefs
curve(my.coefs[1]+my.coefs[2]*x, lwd=2, col=red, add=T)
Cheers,
mitinho astronauta
brazil
On Wed, Sep 10, 2008 at 6:06 PM,
I tried to search for a string of words ``as one entity'' following the
example in the help file:
RSiteSearch({logistic regression})
and got the error message:
2008-09-11 08:55:41.356 open[823] No such file: /Users/rturner/http:/
Hi R users,
I am trying to do two-way clustering (using information of both observation and
variables). Is there any package available in R.
Another querry, if somebody can provide related information (website) regarding
this statistics, it will be great.
Thanks in advance.
Regards,
Ram
Seems that one of the slashed is dropped after the protocol prefix,
i.e. http:/search.r-project.org/... should be
http://search.r-project.org/... Someone else has to take it from
here.
/Henrik
On Wed, Sep 10, 2008 at 2:03 PM, Rolf Turner [EMAIL PROTECTED] wrote:
I tried to search for a
on 09/10/2008 04:03 PM Rolf Turner wrote:
I tried to search for a string of words ``as one entity'' following the
example in the help file:
RSiteSearch({logistic regression})
and got the error message:
2008-09-11 08:55:41.356 open[823] No such file:
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of David Epstein
Sent: Wednesday, September 10, 2008 1:06 PM
To: r-help@r-project.org
Subject: [R] re flecting a line
Suppose x and y are numeric vectors of the same length.
plot(x,y) #scatterplot
On Wed, 10 Sep 2008, Henrik Bengtsson wrote:
Seems that one of the slashed is dropped after the protocol prefix,
i.e. http:/search.r-project.org/... should be
http://search.r-project.org/... Someone else has to take it from
here.
/Henrik
On Wed, Sep 10, 2008 at 2:03 PM, Rolf Turner [EMAIL
On Wed, 10 Sep 2008, Prof Brian Ripley wrote:
As to why the list of links known by name is as it is, that seems history. in
part the White Book history of S. I've always thought it an error that 'log'
was a standard link for binomial, as the range does not match the
specification of
Hi Erika,
I have not tried this before, and I hope that somebody will correct
me if I'm wrong, but the glmer function in the lme4 library appears to do
what you want. From examples(lmer):
lmer (gm1 - glmer(cbind(incidence, size - incidence) ~ period + (1 |
herd), family = binomial, data
On 11/09/2008, at 9:29 AM, Marc Schwartz wrote:
on 09/10/2008 04:03 PM Rolf Turner wrote:
I tried to search for a string of words ``as one entity''
following the
example in the help file:
RSiteSearch({logistic regression})
and got the error message:
2008-09-11 08:55:41.356 open[823]
Hi,
How can I make the cor(x, y, method=spearman) call to produce an
error when the input to it (x, y) produces an error? Here is a simple
example:
a - c(0, 1, 2)
b - c(100, 2, 4)
## error:
log(a)
[1] -Inf 0.000 0.6931472
## error, as expected:
cor(log(a), log(b), method=pearson)
Thomas Lumley wrote:
As a side note, in (some versions of) S, due to partial matching,
binomial(log) is valid -- it just does logistic regression.
ouch!
signature.asc
Description: OpenPGP digital signature
__
R-help@r-project.org mailing list
Hi all,
I've uploaded STAR (Spike Train Analysis with R) on CRAN two days ago.
The package is designed to analyze neuronal spike (action potential)
trains. It uses S3 classes and methods and makes heavy use of other
CRAN packages like gss, R2HTML, mgcv, survival.
* Analysis of both spontaneous
1 - 100 of 129 matches
Mail list logo