Hi all!
I'm doing a coxph analyses on some 50.000 subjects. Is there a simple
function in R that provide some general model information like:
Summary of the number of event and Cencored values like in SAS?
I'm working together with someone using SAS and the summary function in
R does not
This is an interesting question.
What is the problem you are trying to solve and how do the
boundary conditions function as part of this system?
I ask because the asymptotic theory behind your formula for 's.e.'
breaks down with parameters at boundaries. It assumes that you
Hello,
Today I run a R script that I didn't use for a while (one or two months)
and I got an error which I don't understand.
source(run.R)
*** glibc detected *** double free or corruption (!prev): 0x0809ea98 ***
Abandon (core dumped)
Does anybody have a clue about the possible
I'm trying to make a loop with many files...
library(dplR)
files - system(ls *.rwl, intern=TRUE)
files
[1] cimfasy.rwl rocquce.rwl
for (i in files) {a - read.rwl(i,header=0)}
There are 70 series
There are 21 series
class(a)
[1] data.frame
This loop import all the files rwl in a single
Hi,
When I run the following code,
r - c(3,4,4,3,5,4,5,9,8,11,12,13)
n - rep(15,12)
x - c(0, 1.1, 1.3, 2.0, 2.2, 2.8, 3.7, 3.9, 4.4, 4.8, 5.9, 6.8)
x - log10(x)
fr - function(c, alpha, beta) {
P - c + (1-c) * pnorm(alpha + beta * x)
P - pmax(pmin(P,1),0)
-(sum(log(choose(n,r))) + sum(r *
Alfredo Alessandrini wrote:
I'm trying to make a loop with many files...
library(dplR)
files - system(ls *.rwl, intern=TRUE)
files
[1] cimfasy.rwl rocquce.rwl
for (i in files) {a - read.rwl(i,header=0)}
There are 70 series
There are 21 series
class(a)
Is it possible to temporarily suspend output? I am using a function that gives
lots of output using cat, that is not very informative. Is it possible to turn
output off, run the function, then turn output back on again? I've looked at
?options and done various RSiteSearch()s but couldn't find
Alfredo Alessandrini wrote:
can you *explain* what the goal is?
I want import any rwl files (cimfasy.rwl, rocquce.rwl, ...), in a
data.frame with the name like to name of file rwl:
cimfasy.rwl - cimfasy
rocquce.rwl - rocquce
with this loop:
library(dplR)
files -
Richard Pearson wrote:
Is it possible to temporarily suspend output? I am using a function
that gives lots of output using cat, that is not very informative. Is
it possible to turn output off, run the function, then turn output
back on again? I've looked at ?options and done various
On Tue, 24 Jun 2008, Jinsong Zhao wrote:
Hi,
When I run the following code,
r - c(3,4,4,3,5,4,5,9,8,11,12,13)
n - rep(15,12)
x - c(0, 1.1, 1.3, 2.0, 2.2, 2.8, 3.7, 3.9, 4.4, 4.8, 5.9, 6.8)
x - log10(x)
So x[1] = -Inf
fr - function(c, alpha, beta) {
P - c + (1-c) * pnorm(alpha + beta * x)
dear Xu,
does:
library(urca)
example(ur.ers)
ers.gnp
str(ers.gnp)
[EMAIL PROTECTED]
,do what you want?
(this reminds me that I have to learn S4 sometime)
best,
Gustaf Rydevik
On Tue, Jun 24, 2008 at 3:52 AM, Xu, Ke-Li [EMAIL PROTECTED] wrote:
Dear Sir/Madam,
I found your email address and
On Tue, 24 Jun 2008, pqs wrote:
Hello,
Today I run a R script that I didn't use for a while (one or two months)
and I got an error which I don't understand.
source(run.R)
*** glibc detected *** double free or corruption (!prev): 0x0809ea98 ***
Abandon (core dumped)
Does anybody
Dear list:
Given
str(b3)
'data.frame': 159 obs. of 6 variables:
$ index_pollution : num 8.228 10.513 0.549 0.915 10.416 ...
$ position_descrip: chr 2 2 2 NA ...
$ position_geo: chr 3 0 3 3 ...
$ institution : Factor w/ 3 levels digesa,mem,..: 3 3 3 3 3 3
3 3 3 3 ...
$
Dear List Members
As a relatively new R user I am extremely glad for there being a nice
tool such as ggplot for producing easily nice graphs in R. Now I want to
produce graphs with two y-axis. I know that this has been asked a
hundred times with respect to standard plots in R. (The answer, to
A comment from the peanut gallery- It is always good to write an explicative
title because it facilitates searching in the forums.
thanks
Stephen
On Tue, Jun 24, 2008 at 6:49 AM, Gustaf Rydevik [EMAIL PROTECTED]
wrote:
dear Xu,
does:
library(urca)
example(ur.ers)
ers.gnp
str(ers.gnp)
Another question related with the same code:
r - c(3,4,4,3,5,4,5,9,8,11,12,13)
n - rep(15,12)
x - c(0, 1.1, 1.3, 2.0, 2.2, 2.8, 3.7, 3.9, 4.4, 4.8, 5.9, 6.8)
x - log10(x)
fr - function(c, alpha, beta) {
P - c + (1-c) * pnorm(alpha + beta * x)
P - pmax(pmin(P,1),0)
-(sum(log(choose(n,r)))
Wacek, many thanks! I'm wondering however whether this will be problematic on
Windows (I have no windows box to hand to check this, but am creating a package
that I would like to be cross-platform)?
Richard.
Wacek Kusnierczyk wrote:
Richard Pearson wrote:
Is it possible to temporarily
Agustin Lobo wrote:
Dear list:
Given
str(b3)
'data.frame':159 obs. of 6 variables:
$ index_pollution : num 8.228 10.513 0.549 0.915 10.416 ...
$ position_descrip: chr 2 2 2 NA ...
$ position_geo: chr 3 0 3 3 ...
$ institution : Factor w/ 3 levels digesa,mem,..: 3 3
Instead of prepend or append new columns to a matrix, how to insert them to a
matrix ? For example, I would like to insert 3 new columns after the 5th column
of matrix m.
_
[[elided Hotmail spam]]
[[alternative HTML
On Tue, 24 Jun 2008, Jinsong Zhao wrote:
Another question related with the same code:
r - c(3,4,4,3,5,4,5,9,8,11,12,13)
n - rep(15,12)
x - c(0, 1.1, 1.3, 2.0, 2.2, 2.8, 3.7, 3.9, 4.4, 4.8, 5.9, 6.8)
x - log10(x)
fr - function(c, alpha, beta) {
P - c + (1-c) * pnorm(alpha + beta * x)
P -
Hi all,
My application usually works on huge amounts of data in the size of MB's .So
the performance of R in this case will increase if Huge pages are used
instead of malloc'ing (as in memory.c) .
So Are there any methods to support Huge Pages in R ?
Can any one give some inputs to help me
On Tue, 24 Jun 2008, Agustin Lobo wrote:
Dear list:
Given
str(b3)
'data.frame': 159 obs. of 6 variables:
$ index_pollution : num 8.228 10.513 0.549 0.915 10.416 ...
$ position_descrip: chr 2 2 2 NA ...
$ position_geo: chr 3 0 3 3 ...
$ institution : Factor w/ 3 levels
On Tue, 24 Jun 2008, Richard Pearson wrote:
Wacek, many thanks! I'm wondering however whether this will be problematic on
Windows (I have no windows box to hand to check this, but am creating a
package that I would like to be cross-platform)?
It will fail. But
sink(tempfile())
...
sink()
Try something like this:
m - matrix(1:20, ncol = 5)
cbind(m[, 1:3], matrix(1:16, ncol = 4), m[,4:ncol(m)])
On Tue, Jun 24, 2008 at 8:02 AM, Daren Tan [EMAIL PROTECTED] wrote:
Instead of prepend or append new columns to a matrix, how to insert them to
a matrix ? For example, I would like to
Dear all,
I am writing you for two problems I can not solve with R. I used both
the 2.6.1 version and the 2.7.0 version.
I run a sintax written by a collegue of mine, where there are involved
both sparse matrix and quantile estimations. My data begins with
10008*14 observations and at maximum
On Tue, Jun 24, 2008 at 1:15 PM, Prof Brian Ripley
[EMAIL PROTECTED] wrote:
On Tue, 24 Jun 2008, Richard Pearson wrote:
Wacek, many thanks! I'm wondering however whether this will be problematic
on Windows (I have no windows box to hand to check this, but am creating a
package that I would
Dear Maria,
On 23-Jun-08, at 11:40 PM, Maria Gavilanez wrote:
Hi.
is there a way to transform a table that has numeric and categorical
variables into a matrix?
No, because all of the elements of a matrix must be of the same mode,
but a dataframe is a data structure that can have
On Tue, 24 Jun 2008, hadley wickham wrote:
On Tue, Jun 24, 2008 at 1:15 PM, Prof Brian Ripley
[EMAIL PROTECTED] wrote:
On Tue, 24 Jun 2008, Richard Pearson wrote:
Wacek, many thanks! I'm wondering however whether this will be problematic
on Windows (I have no windows box to hand to check
On Tue, 24 Jun 2008 11:50:44 +0100 (BST), Prof Brian Ripley
[EMAIL PROTECTED] wrote:
Some compiled code used by your script screwed up the memory allocator,
directly or indirectly. But we have no idea what is in 'run.R', and little
about your OS, version of R ... see the posting guide.
Or:
tc - textConnection(NULL, w)
sink(tc)
...
sink()
close(tc)
That is a lot more expensive and subject to
o Use of textConnection(NULL, mode=w) could segfault.
(NEWS for 2.7.1). Output textConnections are convenient, but have quite a
lot of overhead since they need to keep the
Thorsten Vogel thorsten.vogel at staff.hu-berlin.de writes:
Dear List Members
As a relatively new R user I am extremely glad for there being a nice
tool such as ggplot for producing easily nice graphs in R. Now I want to
produce graphs with two y-axis. I know that this has been asked a
If anyone remains interested, the solution in base graphics is to modify
stripchart.default, the last couple of lines where the coloring of points
defaults in a way that depends on groups. In my example, the groups are
being handled collectively with the coloring. Code is below.
Deepayan has
Dear all
I just encountered a strange problem with nls formula. I tried to use nls
in cycle but I was not successful. I traced the problem to some parse
command.
Here is an example
DF-data.frame(x=1:10, y=3*x^.5+rnorm(10))
coef(lm(log(DF[,2])~log(DF[,1])))
(Intercept) log(DF[, 1])
Hello,
I am constructing a boxplot but have a very wide range of values (zero -
~28000). I have placed an axis break at the 8000 mark but would like to have a
different scale above the break that ranges from 8000-28000. Right now my axis
is so large that the boxplots are only represented as
just cbind the cols in the appropriate order:
m.2 = cbind( m.1[,1:5], yourthreecolumns, m.1[,6:ncol(m.1)] )
on 06/24/2008 07:02 AM Daren Tan said the following:
Instead of prepend or append new columns to a matrix, how to insert them to a matrix ? For example, I would like to insert 3 new
Hi all,
Suppose I have 2 matrices A and B.
And I want to measure how good each of this matrix is.
So I intend to compare A and B with another gold standard
matrix X. Meaning the more similar a matrix to X the better it is.
What is the common way in R to
measure matrix similarity (ie. A vs X,
Gundala Viswanath wrote:
Hi all,
Suppose I have 2 matrices A and B.
And I want to measure how good each of this matrix is.
So I intend to compare A and B with another gold standard
matrix X. Meaning the more similar a matrix to X the better it is.
What is the common way in R to
measure
on 06/24/2008 09:41 AM Gundala Viswanath wrote:
Hi all,
Suppose I have 2 matrices A and B.
And I want to measure how good each of this matrix is.
So I intend to compare A and B with another gold standard
matrix X. Meaning the more similar a matrix to X the better it is.
What is the common way
I would like to add error bars to a bar chart, I have created in R. I am able
to add error bars to a bar plot, but the same method does not seem to work
for my bar chart version.
Â
Is there a way to add error bars to bar charts?
I followed the instruction of building R package on PC by Peter Rossi. I
got
problem when I tried to build a R package from his example, test.zip. The
problem occurs when I tried to
R CMD check test
the screen stops at
D:\Rpackagebuild\Rossi_testR CMD check test
* checking for working
Hi
I created a new table by extracting only certain rows of table, but the row
numbers in my new table correspond to those from the original table. Is there a
way to reset the row numbers in my new table so that they start from one? like
below:
my table:
COL1 COL2
17 v 45
18 b
Many thanks Brian and Hadley, the quiet function works well for me. Maybe a
candidate for utils?
Best wishes
Richard.
hadley wickham wrote:
Or:
tc - textConnection(NULL, w)
sink(tc)
...
sink()
close(tc)
That is a lot more expensive and subject to
o Use of textConnection(NULL, mode=w)
I don't know the problem per se, but I can say this. In the past, I
used the Rossi example and wasn't able to build an R package for
Windows. I also explored other resources on the web and came to a
similar halt.
The only solution for me was to read the R manual and follow the Murdoch
is it possible to change the position of the ylab in a lavelplot from
the default (left) to the right?
thanks
sw
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Nina, these are not row NUMBERS, but row NAMES. Numbers are actually
reset, they always start with 1 and they are continuous. Just try
doing
T[1,]
on your table. If you want to reset row names, you can do this:
rownames(T) - seq(length=nrow(T))
or you can even remove them:
rownames(T) -
I have a set of data in the form
x1, y1, z1
x1, y2, z2
...
x1, yN, zN
x2, y1, z(N+1)
x2, y2, z(N+2)
...
x2, yN, z(2N)
...and so on...
xM, yN, val(M*N)
I have been trying to figure out how to get R to use this data in a
persp plot. So far the only thing that I can figure out to do is to
Megan J Bellamy bellammj at gov.ns.ca writes:
Hello,
I am constructing a boxplot but have a very wide range of values (zero -
~28000). I have placed an axis break at
the 8000 mark but would like to have a different scale above the break that
ranges from 8000-28000. ...
See gap.plot
I assume that by table you mean data frame. Using the built in
data frame, BOD, try this:
ix - 2:4
BOD[ix, ]
Time demand
22 10.3
33 19.0
44 16.0
data.frame(BOD[ix, ], row.names = seq_along(ix))
Time demand
12 10.3
23 19.0
34 16.0
On Tue, Jun 24, 2008
JohnF == John Fox [EMAIL PROTECTED]
on Tue, 24 Jun 2008 08:10:43 -0400 writes:
JohnF Dear Maria, On 23-Jun-08, at 11:40 PM, Maria
JohnF Gavilanez wrote:
Hi.
is there a way to transform a table that has numeric and
categorical variables into a matrix?
JohnF
Hello,
how can I fit a three, as opposed to a two parameter Gamma distribution,
i.e. one in which shape and rate are not tied.
Thanks,
Serguei
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Dear Gustaf and Stephen,
The lines work - Thank you very much. $B!H([EMAIL PROTECTED](J is exactly
what I want.
I think the forums are very helpful. Could you kindly let me know how to be a
member since my message was held because of $B!H(JPost by non-member to a
members-only
Dear Keli,
To be a member of the R mailing list, just go to
https://stat.ethz.ch/mailman/listinfo/r-help
HTH,
Jorge
2008/6/24 Xu, Ke-Li [EMAIL PROTECTED]:
Dear Gustaf and Stephen,
The lines work - Thank you very much. @teststat is exactly what I want.
I think the forums are very
Serguei,
You might want to look at the TransformedGamma or InverseTransformedGamma
distributions, which are part of the actuar package, see page 41, 25 here:
http://cran.r-project.org/web/packages/actuar/actuar.pdf
You can then use fitdistr of the MASS package to fit the parameters.
Regards,
I would like to filter using a periodogram produced by spec.pgram() then
pick the frequencies that I would like to filter out and then get the
inverse back into the time domain. Any thoughts.
Stephen
--
Let's not spend our time and resources thinking about things that are so
little or so
Does this example help?
my.values - c(10, 5)
err1 - c(.5,.3)
x - barplot(my.values, ylim=c(0,12))
arrows(x,my.values-.5 ,x,my.values+.5, code=3, angle=90, length=.1)
box()
--- On Tue, 6/24/08, Liat RT [EMAIL PROTECTED] wrote:
From: Liat RT [EMAIL PROTECTED]
Subject: [R] Bar charts with
Hi All,
The for-loop below stopped when error(Cannot get confidence intervals on
var-cov components: Non-positive definite approximate variance-covariance)
occurred.
I assigned a row of NA values
to the data frame m1 manually and reset j in the for-loop every time error
returned. Im
R 2.7.0, WinXP
Hi list,
I have written some code in R that I would like to share with others who
are not R-users. Is there a way to compile the code so it will run
outside of R?
The best thing would be a functionality to compile the code into a .exe
file that does not require the user to have
Hi,Dear all R experts,
I am trying to do the 2-way contingency table analysis by fitting the
loglinear models. However, I found my table has several empty cells which are
theoretically missing values.I have no idea of how to solve them coz we cannot
compute the simulated p-value with
?try
On Tue, Jun 24, 2008 at 1:59 PM, Sit Yee Kong [EMAIL PROTECTED] wrote:
Hi All,
The for-loop below stopped when error(Cannot get confidence intervals on
var-cov components: Non-positive definite approximate variance-covariance)
occurred.
I assigned a row of NA values
to the data
Hello,
Is there some equivalent function in R of Mathematica's Apply
i.e.
Apply[ f, g[a,b]] == f[a,b]
In concret, I have a matrix of indexes:
ind = cbind(1:5, 1:5+5)
ind
[,1] [,2]
[1,]16
[2,]27
[3,]38
[4,]49
[5,]5 10
and I'd like to extract
On 24/06/2008, at 8:10 AM, Daniel Folkinshteyn wrote:
on 06/23/2008 03:40 PM Thomas Frööjd said the following:
SNIP
2. Scale the data so they can be plotted on the same axis. The
reference dataset has around 20 000 observations and my data from the
clinic only around 3000
Thanks for your reply.
When I did what you said
z - t(matrix(zdata, nrow = length(ydata)))
persp(xdata, ydata, z)
I got the following errors:
Error in persp.default(xdata, ydata, z) : invalid 'x' argument
Calls: persp - persp.default
In addition: Warning messages:
1: In min(x, na.rm =
hi,
I have a data stored in table format. The data is stored in the name map.
500545 rows and 3 columns.
I want to use a for loop to scan through till the end of the row.
can i use like this.
for ( i in 1:row(map)
{
}
Kindly clarify me
Cheers
Ramya
--
View this message in context:
Hi Ramya.
Like this?
for ( i in 1:nrow(map) )
{
cat (this is the line, i, \n)
}
Miltinho.
Brazil
On 6/24/08, Rajasekaramya [EMAIL PROTECTED] wrote:
hi,
I have a data stored in table format. The data is stored in the name map.
500545 rows and 3 columns.
I want to use a for loop to
Try
apply(ind, 1, function(i) x[i[1]:i[2]])
which returns a matrix if all sequences are the same length, a list
otherwise
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Juan Pablo Romero Méndez
Sent: Tuesday, June 24, 2008 11:51 AM
To: [EMAIL
On 25/06/2008, at 6:01 AM, Roger Leenders wrote:
R 2.7.0, WinXP
Hi list,
I have written some code in R that I would like to share with
others who
are not R-users. Is there a way to compile the code so it will run
outside of R?
The best thing would be a functionality to compile the code
Dear All,
I want to plot a heat map with annotated columns. Both functions heatmap
(stats) and heatmap.2 (gplots) can plot a horizontal side bar that can
be used to visualize a categorical variable. In addition to a
categorical variable, I would like to visualize a continuous variable.
This
On 6/24/08, Bryan Hanson [EMAIL PROTECTED] wrote:
If anyone remains interested, the solution in base graphics is to modify
stripchart.default, the last couple of lines where the coloring of points
defaults in a way that depends on groups. In my example, the groups are
being handled
Hi Gundala,
Suppose I have 2 matrices A and B.
And I want to measure how good each of this matrix is.
You really want to be using Robert Escoufier's RV-coefficient (A unifying
tool for linear multivariate statistical methods: The $RV$-coefficient Appl.
Statist., 1976, 25, 257-265).
Several
I cannot modify the font size of the standard text items for the
splineplot function
found in the standard graphics package. Here is a tiny example, in
which
I try (but fail) to make the main title large.
tbl = read.table ('http://www.ggobi.org/book/data/tips.csv',
sep=',', header=T)
I need some help with a genome-wide association scatterplot. A
recently published example of such a plot is attached to this message.
Typical variables include:
Y: -logP-value (range from 0 to 8)
X: chromosomal distance
Z: chromosomes 1 through 22 and X (a factor with 23 levels)
I used
I am indebted to Deepayan Sarkar for lattice graphics, Fritz Leisch
for Sweave and the authors of the beamer package for LaTeX. For me
these are a killer app combination. I wouldn't want to go back to
creating presentations in any other way.
When I create a presentation I use a Makefile that
On 24/06/2008 5:21 PM, Paul Shannon wrote:
I cannot modify the font size of the standard text items for the
splineplot function
found in the standard graphics package. Here is a tiny example, in
which
I try (but fail) to make the main title large.
tbl = read.table
Hi,
Here's the data we have:
rs[1:5,]
probe_id f1 f2 f3 f4 MA f
1 A_68_P20002076 2 58 0 0 1.51778114 6.344453 59
2 A_68_P20002775 22 8 15 0 0.43419304 5.488819 59
3 A_68_P20005791 43 3 0 0 0.05698666 9.830594 59
4 A_68_P20005805 11 34 0 0 1.71076835 6.624038
Hi!
I am working on a line graph under xyplot. what command should I feed
to auto.key so that it uses the line pattern of series as labels?
Also where can I find doc about all options that I can feed to auto.key?
Thanks!
Mark
__
Is this what you want:
x - read.table(textConnection( probe_id f1 f2 f3 f4 M
A f
+ 1 A_68_P20002076 2 58 0 0 1.51778114 6.344453 59
+ 2 A_68_P20002775 22 8 15 0 0.43419304 5.488819 59
+ 3 A_68_P20005791 43 3 0 0 0.05698666 9.830594 59
+ 4 A_68_P20005805 11 34 0
What do you mean by A similar to X?
Do you mean norm of the difference, similar eigenvalues/vectors, anything else?
--- On Wed, 25/6/08, Gundala Viswanath [EMAIL PROTECTED] wrote:
From: Gundala Viswanath [EMAIL PROTECTED]
Subject: [R] Measuring Goodness of a Matrix
To: [EMAIL PROTECTED]
Hi Anh,
I know that it is not so elegant, but may work.
x-runif(10)
y-runif(10)*2
z-runif(10)*3
df-data.frame(cbind(x,y,z))
df
df$f-NULL
for (i in 1:nrow(df))
{
df[i,f]-max(c(df[i,x],df[i,y],df[i,z]))
}
Of course, others that know better apply family will solve if with a line
:-)
There
Dear R help,
In the middle of my session, I started receiving the error
message bad value
regardless of what I entered. I had to close R and restart
it. When the error
occurred, I was creating basic lattice plots, with no other
add-on packages running.
I am running R through John Fox's XEmacs
Thanks all!
JP
On Tue, Jun 24, 2008 at 3:30 PM, Franz Mueter [EMAIL PROTECTED] wrote:
Try
apply(ind, 1, function(i) x[i[1]:i[2]])
which returns a matrix if all sequences are the same length, a list
otherwise
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
?pmax
Anh Tran wrote:
Hi,
Here's the data we have:
rs[1:5,]
probe_id f1 f2 f3 f4 MA f
1 A_68_P20002076 2 58 0 0 1.51778114 6.344453 59
2 A_68_P20002775 22 8 15 0 0.43419304 5.488819 59
3 A_68_P20005791 43 3 0 0 0.05698666 9.830594 59
4 A_68_P20005805 11 34
Hi everyone,
I'm sorry if this turns out to be more a statistical question than one
specifically about R - but would greatly appreciate your advice anyway.
I've been using a logistic regression model to look at the relationship
between a binary outcome (say, the odds of picking n white balls
Thanks so much all.
I will try out each of your suggestions.
- Gundala Viswanath
Jakarta - Indonesia
On Wed, Jun 25, 2008 at 8:01 AM, Moshe Olshansky [EMAIL PROTECTED] wrote:
What do you mean by A similar to X?
Do you mean norm of the difference, similar eigenvalues/vectors, anything
else?
A short answer to your question is 'yes and no': R is an
interpreted language. There may have been attempts to create a
compiler for R, but I don't know if any have actually succeeded very
well.
However, there are ways to (a) call R from other languages and (b)
run R in batch
It looks like A*B*C*D is a complete, totally saturated model, (the
residual deviance is effectively zero, and the residual degrees of
freedom is exactly zero - this is a clue). So when you try to put even
more parameters into the model and even higher way interactions,
something has to give.
Roger,
Following on to Spencer's comments, we have had some success using RWui
http://rwui.cryst.bbk.ac.uk/
to build web enabled versions of (non-interactive) R programs.
Advantage is your users only need a browser. No exe to install. Downside is
you need to have a web server accessible to your
You can set up a web site and have them use that. Then they
just need a browser.
See www.rpad.org for one way or develop an internal or external site
using any technology you know and call R on the back end.
On Tue, Jun 24, 2008 at 2:01 PM, Roger Leenders [EMAIL PROTECTED] wrote:
R 2.7.0,
On 24/06/2008 3:12 PM, Toby Gass wrote:
Dear R help,
In the middle of my session, I started receiving the error
message bad value
regardless of what I entered. I had to close R and restart
it. When the error
occurred, I was creating basic lattice plots, with no other
add-on packages running.
Hello Mr. Graves,
Hello all useRs,
Many thanks for your attention.
This is an interesting question.
What is the problem you are trying to solve and how do the
boundary conditions function as part of this system?
One of the functions that I need to minimized is:
sum( ( ( hat_xi - xi
Hello Hans,
Do you need *both* a categorical *and* a continuous variable, or just a
continuous variable?
-Greg
On 6/24/08 4:51PM , Hans-Ulrich Klein [EMAIL PROTECTED] wrote:
Dear All,
I want to plot a heat map with annotated columns. Both functions heatmap
(stats) and heatmap.2 (gplots)
Hi,
I have two matrices x and y (same dimensions).
How can I compute the average of x and y for each coordinate
to form a new average matrix z ?
- Gundala Viswanath
Jakarta - Indonesia
__
R-help@r-project.org mailing list
What about z - (x + y) / 2
Gundala Viswanath wrote:
Hi,
I have two matrices x and y (same dimensions).
How can I compute the average of x and y for each coordinate
to form a new average matrix z ?
- Gundala Viswanath
Jakarta - Indonesia
__
Thank you for the response. I was doing a basic operation
with a small dataset and am (fortunately for me) unable to
reproduce the error. After restarting R, without shutting
down anything else, there was no error. I thought I would
post, however, since 2 others had reported the same
Mark Kimpel wrote:
I'm using the command below on an open3d() object to create a shaded
cube. Changes to myScalingFactor do not effect changes in the size of
the cube. What is the correct approach? Mark
shade3d(translate3d(cube3d(identityMatrix() *
myScalingFactor),-6,1,-1),col=green, alpha =
Mark Kimpel mwkimpel at gmail.com writes:
I'm using the command below on an open3d() object to create a shaded
cube. Changes to myScalingFactor do not effect changes in the size of
the cube. What is the correct approach? Mark
how about scale3d() ?
You can specify units= to be extra sure:
d - difftime(Sys.Date() + 1, Sys.Date(), units = mins)
as.numeric(d, units = days)
On Wed, Jun 25, 2008 at 12:23 AM, Richard M. Heiberger [EMAIL PROTECTED]
wrote:
I received a warning about incompatible methods that I don't understand.
The workaround
The user of your R script sees only the outputs you create. The R source is
hidden.
HTH,
Jim Porzak
Responsys, Inc.
San Francisco, CA
http://www.linkedin.com/in/jimporzak
On Tue, Jun 24, 2008 at 9:33 PM, Esmail Bonakdarian [EMAIL PROTECTED]
wrote:
Jim Porzak wrote:
Roger,
Following on to
97 matches
Mail list logo
