On Thu, 31 May 2007, Saptarshi Guha wrote:
Hi,
I am writing a package which requires the ANN(Approx N.Neighbors)
library (found here http://www.cs.umd.edu/~mount/ANN/).
My package directory looks something like this
atry/
atry/src/ann/ --ann source is here
Dear R users
I have 3 classification variables (string type) and one string variable in a
data.frame (x).
the first factor is a region (north, south, west, and east) (4 levels)
the second factor is a market type (A, B, C, and D) (4 levels)
the third factor is a item type (w, x, y, z) (4 levels)
On 2007-June-01 , at 01:03 , Tom La Bone wrote:
The function wtd.var(x,w) in Hmisc calculates the weighted variance
of x
where w are the weights. It appears to me that wtd.var(x,w) = var
(x) if all
of the weights are equal, but this does not appear to be the case. Can
someone point out
Charles C. Berry wrote:
On Thu, 31 May 2007, Jacques Wagnor wrote:
Dear List,
I have a simple model as follows:
x - rnorm(500)
library(logspline)
fit - logspline(x)
n - 100
y - replicate(n, sum(rlogspline(rpois(1,10), fit))) # last line
What you do not see in RGui is this:
Dears, I have the below code for metropolis of the GLM logit (logistic
regression) using a flat prior. Can someone help me modify the prior so that
the model becomes hierarchical by using a flat prior for mu and sigma, the
derived density for beta ~ N(mu, sigma^2)? Actually I took my code from a
look at the `arr.ind' argument of ?which(), e.g.,
x - matrix(rnorm(9), 3, 3)
x
which(x == max(x), arr.ind = TRUE)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35,
zhijie zhang wrote:
Dear friends,
I need the get the 2.5% and 97.5% quantile from each row of a matrix, how
should i get it?
BTW, i can get the min/max value from each row of a matrix, using the
following programs, is there an easy function to do it?
simmin-matrix(NA,nrow=47,ncol=1)
The most 'relevant document' is ?download.file, and there is AFAIK
nothing specific to MacOS X about the code (and I did write it).
You can use
Sys.setenv(http_proxy=http://un:[EMAIL PROTECTED]:port)
at the start of your R session, or use
http_proxy=http://un:[EMAIL PROTECTED]:port
in your
Hi,
I have this dataset where columns z1.3 and z1.4 are full of NA's. I want
to perform some calculations in the remaining columns, but after doing
this, I want to recontruct the original matrix. I can with:
out - which( colMeans( is.na( z ) ) == 1 )
gd-z[, - out]
select the columns full of
Dear R Users,
I would like reduce font size of the plots of my fanny result, but I
count not. Is there a way to make changes on font sizes?
Thanks.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Dear friends,
I need the get the 2.5% and 97.5% quantile from each row of a matrix, how
should i get it?
BTW, i can get the min/max value from each row of a matrix, using the
following programs, is there an easy function to do it?
simmin-matrix(NA,nrow=47,ncol=1)
for (i in 1:47) {
Hi
I have a problem with computing the logrank test using R, can someone give
me the relavent code or help me otherwise please!
Thanks
--
View this message in context:
http://www.nabble.com/Logrank-test-tf3851514.html#a10910383
Sent from the R help mailing list archive at Nabble.com.
I think nearly this same question was just asked by someone within
the last few days:
Anyways, this multiplies gd by 10 and then places it back into the
original columns:
z[, -out] - 10*gd
On 6/1/07, antonio rodriguez [EMAIL PROTECTED] wrote:
Hi,
I have this dataset where columns z1.3 and
Fränzi Korner fraenzi.korner at oikostat.ch writes:
how do I obtain standard errors of the predicted values from a lme (or
lmer)-object?
Not totally clear what you mean. intervals(lmeresult) gives the confidence
intervals for the coefficients. Otherwise, you can do some calculations with
I use:
foo - function(...) {
args - list(...);
names(args);
}
/Henrik
On 6/1/07, Mike Meredith [EMAIL PROTECTED] wrote:
Is there a tidy way to get the names of objects passed to a function via the
... argument?
rbind/cbind does what I want:
test.func1 - function(...) {
nms -
Hello, I want to conduct normality test to a series of data and get the
p-value for each subset. I am using the following codes, but it does not
work.
tapply(re, list(reg, ast), pvalue(shapiro.test))
Could anyone give me some advice? Many thanks.
--
View this message in context:
Thanks, Henrik, but 'foo' doesn't do what I want:
x - some stuff
second - more stuff
foo(first=x, second)
[1] first
Brian's right:
...he wants the argument name if there is one otherwise the
deparsed argument value, but clarification would be helpful.
The function using this compares
Dear all,
I would like to generate a regular time serie, i.e. a list of dates and
time for each our of the period 2002-2004.
the time format should be
2002-01-01 12:00:00 (year-month-day hour:min:sec)
so the list should contain all hours of the period 2002-2004
2002-01-01 00:00:00
2002-01-01
On Fri, 1 Jun 2007, Henrik Bengtsson wrote:
I use:
foo - function(...) {
args - list(...);
names(args);
}
But that does not do what was asked: it gets the argument names, not the
object names. (Did you actually try it?) It looks from the example that
he wants the argument name if
Hi,
Does R runs under citrix version 4 environnement or terminal services?
Regards
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On 6/1/07, Alfonso Sammassimo [EMAIL PROTECTED] wrote:
Hi R-experts,
Thanks very much to Jim Holtman and Gabor on my previous question.
I am having another problem with data manipulation in zoo. The following is
data (Z) for first business day of every month in zoo format. I am trying to
I'm not sure what is the 'pvalue' function (it's not found in base nor
stats packages) but
this should give you what you want:
# some example
re - rnorm(100)
reg - rep(1:3, length=100)
ast - rep(1:2, length=100)
tapply( re, list(reg, ast), function(v) shapiro.test(v)$p.value )
# or neater by
try this:
tapply(re, list(reg, ast), function(x) shapiro.test(x)$p.value)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax:
Hi
If you have tried to go through help pages of max you could find out which
function, which can tell you position of your maximum.
which(x==max(x), arr.ind=T)
Regards
Petr Pikal
[EMAIL PROTECTED]
[EMAIL PROTECTED] napsal dne 01.06.2007 11:09:09:
Dear UseRs,
I have a very simple
Dear UseRs,
I have a very simple question. I have a big matrix (say x) including
probabilities (values in (0,1)).
I have to store in a list the names of the row and the column where max(x)
is located. How can I proceed?
Thanks in advance for your assistance!
mirko
Hello
i have a model with two random effects. i have used lmer (package lme4)
on the model. To find the predictions of the random effects used
ranef(). How can i get p-values on theese predictions??
Thanks
Rina
[[alternative HTML version deleted]]
On 2007-June-01 , at 13:00 , Tom La Bone wrote:
The equation for weighted variance given in the NIST DataPlot
documentation
is the usual variance equation with the weights inserted. The
weighted
variance of the weighted mean is this weighted variance divided by N.
There is another
jiho wrote:
On 2007-June-01 , at 01:03 , Tom La Bone wrote:
The function wtd.var(x,w) in Hmisc calculates the weighted variance
of x
where w are the weights. It appears to me that wtd.var(x,w) = var
(x) if all
of the weights are equal, but this does not appear to be the case. Can
On 31/05/07, Anup Nandialath [EMAIL PROTECTED] wrote:
oops forgot the example example
try this line
sqrt(mahalanobis(all, colMeans(predictors), cov(all), FALSE)
Hi and thanks for the reply Anup. Unfortunately, I had a look on the
example before posting but not much of a help... I did some
Hello
i have a model with two random effects. i have used lmer (package lme4)
on the model. To find the predictions of the random effects used
ranef(). How can i get p-values on theese predictions??
Thanks
Rina
[[alternative HTML version deleted]]
On Fri, 1 Jun 2007 12:51:35 +0200 [EMAIL PROTECTED] wrote:
Dear all,
I would like to generate a regular time serie, i.e. a list of dates
and time for each our of the period 2002-2004.
use the seq() method for POSIXct objects:
seq(from = as.POSIXct(2002-01-01 00:00:00),
to =
Look here for the answer:
https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html
:
::
::: Note that my e-mail address has changed to m.j.bojanowski at uu dot
nl
::: Please update you address books accordingly. Thank you!
::
:
Michal Bojanowski
ICS /
library(chron)
dt1 - chron(1/1/2002, 00:00:00)
dt2 - chron(1/2/2002, 00:00:00)
chron(seq(as.numeric(dt1), as.numeric(dt2), by = 1/24))
Also (but see caveat in R News 4/1 help desk article):
dt1 - as.POSIXct(2002-01-01 00:00:00)
dt2 - as.POSIXct(2002-01-02 00:00:00)
seq(dt1, dt2, by = hour)
On
Thanks. I have another related question:
The equation for weighted variance given in the NIST DataPlot documentation
is the usual variance equation with the weights inserted. The weighted
variance of the weighted mean is this weighted variance divided by N.
There is another approach to
It certainly runs under Windows terminal services. Since it is just a
standard Windows application I see not reason why it should not run under
Citrix, but have no first-hand experience.
On Fri, 1 Jun 2007, Patrick Haurie wrote:
Hi,
Does R runs under citrix version 4 environnement or
Wonderful! Thanks for the assistance.
Tom La Bone
-Original Message-
From: jiho [mailto:[EMAIL PROTECTED]
Sent: Friday, June 01, 2007 8:17 AM
To: [EMAIL PROTECTED]
Cc: 'R-help'
Subject: Re: [R] Problem with Weighted Variance in Hmisc
On 2007-June-01 , at 13:00 , Tom La Bone wrote:
See:
http://tolstoy.newcastle.edu.au/R/e2/help/06/10/2242.html
which we can modify slightly for the case in question like this:
f - function(...) {
x - list(...)
if (is.null(names(x))) names(x) -
names(x)[names(x) == ] - NA
mc - match.call()[-1]
See ?survdiff in the survival package.
Mbini wrote:
Hi
I have a problem with computing the logrank test using R, can someone give
me the relavent code or help me otherwise please!
Thanks
--
Kevin E. Thorpe
Biostatistician/Trialist, Knowledge Translation Program
Assistant Professor,
Rina
p-values are not generated as a part of the lmer object returned. If you read
the ranef help it will show you two things. First, how you can get a
caterpillar plot and second how you can get the posterior variances of the
random effects. Then, you can do your test to see if they are
Dear members,
I would like to pass the histogram settings to each subset of the dataframe,
and generate a multiple figures graph.
First, can anyone tell me how to generate a multiple figures environment? I
am trying
mfrow=c(2,4) and nothing appears.
Secondly, I want to pass the following
use lattice graph
At 08:00 AM 6/1/2007, livia wrote:
Dear members,
I would like to pass the histogram settings to each subset of the dataframe,
and generate a multiple figures graph.
First, can anyone tell me how to generate a multiple figures environment? I
am trying
mfrow=c(2,4) and
How about
foo - function(...)
{
m - as.list(match.call(expand.dots=TRUE))[-1]
nm - names(m)
for(i in seq_along(m)) if(!nchar(nm[i])) nm[i] - deparse(m[[i]])
nm
}
Such things are hard to do from R level, hence the use of match.call to do
it at C level.
On Fri, 1 Jun 2007, Mike
Dear Users,
I would also like to learn how to change font size of the these two
components... text in the fanny plot, or to remove it? (It maybe so
easy, but you know sometimes can things be hard to find.)
Thanks.
__
R-help@stat.math.ethz.ch mailing
I wrote some rough functions for principal factor,
principal-components factor, and iterated principal factor analysis.
I think they are workable, the same results as stata can be retained.
In addition, functions for gls and uls factor analysis is in progress,
which is based on the algorithms of
http://yenibiris.sendeyolla.com/medyadetay.aspx?tid=3cid=57id=61365
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented,
http://yenibiris.sendeyolla.com/medyadetay.aspx?tid=3cid=57id=61365
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented,
Thanks very much to all of you. It looks like 'match.call' is the key, and
both Brian's and Gabor's solutions work fine. --- Mike.
--
View this message in context:
http://www.nabble.com/Getting-names-of-objects-passed-with-%22...%22-tf3850318.html#a10913978
Sent from the R help mailing list
Dear all,
Concerning the function density from the stats package, I don´t know what is
the dependency between the number of equally spaced points at which the density
is to be estimated (the ´n´ argument) and the data from which the kernel
density estimate is to be computed (the ´x´
On 6/1/07, Mike Meredith [EMAIL PROTECTED] wrote:
Thanks very much to all of you. It looks like 'match.call' is the key, and
both Brian's and Gabor's solutions work fine. --- Mike.
Brian's is shorter but I think the one in my post is a bit more robust:
f1 - function(...) {
+m -
Hi folks,
Is it possible to have Excel call a R function. If not, how about making Excel
send off a command to call a R script and then read the result back into Excel.
I know, I know, this should belong to some Excel forum, but i just try my luck
here.
Thanks in advance.
Horace W. Tso
On Fri, 2007-06-01 at 06:00 -0700, livia wrote:
Dear members,
I would like to pass the histogram settings to each subset of the dataframe,
and generate a multiple figures graph.
First, can anyone tell me how to generate a multiple figures environment? I
am trying
mfrow=c(2,4) and
You might consider having R do everything: R can read the Excel sheet, do
what needs to be done, and write the results to an Excel sheet.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From:
Charles, thanks for the thought. In my case that's quite impractical/impossible
because I'm working from a monstrous spreadsheet having dozens of sheets and
numerous links. In some way that's the irreplaceable aspect of a spreadsheet
program like Excel.
H.
Charles Annis, P.E. [EMAIL
On Fri, 2007-06-01 at 08:50 -0700, Horace Tso wrote:
Hi folks,
Is it possible to have Excel call a R function. If not, how about
making Excel send off a command to call a R script and then read the
result back into Excel.
I know, I know, this should belong to some Excel forum, but i just
Dear R-help list members,
I am new to R, and having problems with plotting part of a time series. I
have read in my data using read.table, and my 'time' column is recognised
as a numeric variable. When I convert this to a time format, I am no longer
able to use its values on my x-axis. The
Horace Tso wrote:
Is it possible to have Excel call a R function. If not, how about
making Excel send off a command to call a R script and then read the
result back into Excel.
I know, I know, this should belong to some Excel forum, but i just
try my luck here.
You can always put R
Marc Schwartz wrote:
On Fri, 2007-06-01 at 08:50 -0700, Horace Tso wrote:
Hi folks,
Is it possible to have Excel call a R function. If not, how about
making Excel send off a command to call a R script and then read the
result back into Excel.
I know, I know, this should belong to some
Dear R users,
I'm pretty new on using lmer package. My response is binary and I have fixed
treatment effect (2 treatments) and random center effect (7 centers). I want
to test the effect of treatment by fitting 2 models:
Model 1: center effect (random) only
Model 2: trt (fixed) + center
HI Emine,
The parentheses specify a variable as random or fixed. You are using
center for random intercepts, but fixed for interactions.
You may want
mer(cbind( yvect, nvect-yvect) ~ trt + ( 1 | center) + (1 |
trt:center),
family = binomial, niter = 25, method = Laplace, control = list
And just for completeness for the thread,
there is the RDCOMServer package which has
a more general mechanism for making R
functionality available to other applications
that does not require writing R commands
but uses the language independent spirit of
DCOM. And it also allows event handlers
for
Hi Patrick,
[EMAIL PROTECTED] wrote:
Dear R-help list members,
I am new to R, and having problems with plotting part of a time series. I
have read in my data using read.table, and my 'time' column is recognised
as a numeric variable. When I convert this to a time format, I am no longer
On Thu, 31 May 2007, Strickland, Matthew (CDC/CCHP/NCBDDD) (CTR) wrote:
Thanks for your reply Charles. I do indeed have other variables. I
apologize for being vague, here is my study in more detail:
I have a cohort of births. My outcome is a dichotomous variable for
presence/absence of a
On 6/1/07, emine özgür Bayman [EMAIL PROTECTED] wrote:
Dear R users,
I'm pretty new on using lmer package. My response is binary and I have fixed
treatment effect (2 treatments) and random center effect (7 centers). I want
to test the effect of treatment by fitting 2 models:
Model 1: center
Hi,
I want to make use of some C routines from other packages to write extensions
in C.
In Writing R Extensions, it says there is an experimental interface to
support this in (or from ?) R 2.4.0.
I had a dummy library containing src/dummy.cpp and R/zzz.R:
== src/dummy.cpp
Dear all,
I'm completely new to R and at first I must say that it is a great program!
But I have a problem with the function isoMDS from the MASS package. I
have this code which I load with source() from a file:
x - c(163.59514923926784, 150.01448475257115, .. {here are some
more values})
I was reading the help but just did not get how to specify starting values for
varIdent() of the lme() function, although I managed to do it for corSymm().
Do I specify the values just as they are printed out in an output, like c(1,
1.3473, 1.0195). Or do I need to take the residual and
It looks like your C++ compiler is mangling the name
R_RegisterCCallable. Move the
extern C
before the
#include R_ext/Rdynload.h
On 6/1/07, Alex Chen [EMAIL PROTECTED] wrote:
Hi,
I want to make use of some C routines from other packages to write extensions
in C.
In Writing R
On 6/1/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
I was reading the help but just did not get how to specify starting values for
varIdent() of the lme() function, although I managed to do it for corSymm().
Do I specify the values just as they are printed out in an output, like c(1,
Hi All,
I would like to extend the standard biplot of Grabriel (1971), which is done
by the R function biplot(). Rather than just showing the variables as
vectors (with arrows at end point), I would like to show them as line
segments, where the segment runs from minimum value of the variable
Thanks to Marc, Duncan, Tobias, Alberto. I'm aware of rcom and have in fact
downloaded it somewhere. But unfortunately it won't work for me as I'm in a
corporate environment and don't have admin rights to set up the R(D)COM
(background) server. If there is a way to work around it, pls let me
Hi, all,
This is for R-2.5.0 on WinXP and in particular RGui.
I'm trying to teach some colleagues of mine R and rather than impose
Xemacs/ESS upon them I decided to simply start by showing them RGui.
When R is installed, R workspaces (.RData) are automatically registered
so that I can
hi, Rers
I tried to find out the max size (# of rows, # of
columns) of a matrix that is allowed by R but failed.
Can anybody let me know?
Thanks!
Guanrao
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
AFAIK, it only depends on how much free memory you have.
b
On Jun 1, 2007, at 5:05 PM, Guanrao Chen wrote:
hi, Rers
I tried to find out the max size (# of rows, # of
columns) of a matrix that is allowed by R but failed.
Can anybody let me know?
Thanks!
Guanrao
There is no maximum size. This will be driven by (at least) two issues.
First, how much memory you have on your own computer and second what
data you have in each cell. For instance, an integer takes less memory
than a floating point.
Other spreadsheet programs like excel limit the number of rows
Harold,
Actually there is a maximum size, even if one had sufficient memory.
Since a matrix is a vector with a dim attribute, and these objects are
indexed using integers, the maximum sized vector one 'could' create is:
.Machine$integer.max
[1] 2147483647
which is:
2^31 - 1
[1] 2147483647
Hi!
When running the following code I get the message: Error in eval(expr,
envir, enclos) : object A not found
fm - function(p,ydata, env = parent.frame()) {
#fit model y ~ (b0+b1x1+...+bpxp)*exp(g2plus*z2plus) where bi and g2plus
are parameters
#create design matrix for linear part
fo00
I have a windows source file.r, with the default charset of windows.
I can't use it in Linux as source(file.r), because Linux's default is
Unicode. How can I read it?
Alberto Monteiro
__
R-help@stat.math.ethz.ch mailing list
Sorry, I responded a bit too hastily last night, without testing the two
functions properly.
Brian's is shorter but I think the one in my post is a bit more robust:
Indeed! 'f1' only works if at least one of the arguments is named. Otherwise
'nm' is NULL and 'nchar(nm[i])' fails. 'f2' seems to
iconv on your linux box should do the work.
b
On Jun 1, 2007, at 7:19 PM, Alberto Vieira Ferreira Monteiro wrote:
I have a windows source file.r, with the default charset of windows.
I can't use it in Linux as source(file.r), because Linux's
default is
Unicode. How can I read it?
Alberto
Dear R-users,
I'm trying to do some spatial simulation. I have two covariates, Z and C. I
want to examine their relationship under different spatial distribution.
I have no problem simulating completely spatial random process but I'm totally
stuck on poisson (cluster) pattern. I already have
Benilton Carvalho wrote:
iconv on your linux box should do the work.
I have a windows source file.r, with the default charset of windows.
I can't use it in Linux as source(file.r), because Linux's
default is
Unicode. How can I read it?
Soon after I posted, I found out how to do it.
Hi - I tried to get some answer before but there seems to have no one response.
My question is that I have a class like below definition:
setClass(DBManager,
representation(channel=RODBC))
My purpose of the conn slot is to hold a channel to database connection which
is returned by a
Hi,
I have a couple of questions about johansen's test, in general:
1. I was able to obtain error correction term (ect) from cajorls$rlm$model
properly. According the my ca.jo object on 2-variate series, the test suggests
that the integration rank is 1. Which means that my ect should be
83 matches
Mail list logo