Thanks, Gabor! This workaround fixes my immediate needs.
--sundar
Gabor Grothendieck said the following on 10/1/2008 10:42 PM:
That may be a bug in R but I think there is another problem on top of that
as I don't think bquote descends into function bodies:
z - 2
bquote(function(x) {x^.(z)})
Hi,
I have tried this with R-2.7.2 then I got
alist[['goodbye']]
NULL
alist[['hi']]
NULL
alist[['hello']]
[1] 10
The results make sense to me thus you might want to update your R version.
chunhao
B Fox wrote:
This seems odd. When I try to look up a list element which has a space in
Hi I believe I have long time solution to the problems with the
read.spss() funtion.
In non-ascii locales the import of spss files into R has consistently
produced errors and sometimes wouldn't allow import at all.
I've made a custom version of the read.spss() function that appear to
rkevinburton at charter.net writes:
ff - complex(length(fs))
ff[9] - fs[9]
ff[5] - fs[5]
Include the DC component:
ff[1] - fs[1]
Take the inverse
fi - fft(ff, inverse=TRUE) / length(ff)
Plot
plot(fi)
Notice that the plot is the Re vs. Im on the x and y axis'
Michael Just mgjust at gmail.com writes:
I am trying to plot multiple histograms with the same scales, etc into one
plot. The commands below produce a 3 page PDF with each histogram occupying
the upper right quadrant. And use slightly different scales on the X and Y
axes.
I suggest that you
Daniela Reicheneder Daniela.Reicheneder at tz.agrar.tu-muenchen.de writes:
I was using the function pvals.fnc from package 'languageR' until April.
I do not know which version. Yesterday I updated all my packages
and tried to run my loop again. Now I get the following error message:
error
Dear R users,
Assume I have three models with the following AIC values:
model AIC df
model1 -10 2
model2 -12 5
model3 -11 2
Obviously, model2 would be preferred, but it wastes 5 df compared to the
other models.
Would it be allowed to select model3 instead, simply because
Dear Michael,
Try ggplot2. Use something like
install.packages(ggplot2)
library(ggplot2)
recov - 0:2
n - 1000
all - data.frame(bbED = rnorm(3 * n, mean = recov), recov =
factor(rep(recov, n)))
ggplot(data = all, aes(x = bbED)) + geom_histogram() + facet_grid(. ~
recov)
ggplot(data = all, aes(x =
Duncan Murdoch wrote:
Daniel Rabczenko wrote:
Hello everybody,
Two - I hope Simple questions about working with workspaces.
Is there a way to force R to start in clean workspace / avoid
previously
saved workspace restored?
Start with --no-restore. (Start with --help for the full list
Christoph Scherber wrote:
Dear R users,
Assume I have three models with the following AIC values:
modelAICdf
model1 -102
model2-125
model3-112
Obviously, model2 would be preferred, but it wastes 5 df compared to
the other models.
Would it be allowed to
ram basnet wrote:
Dear R users,
I am using RandomForest package. While using this package, i got
Error: cannot allocate vector of size 117.3 Mb .message.
I had this problem earlier too but could not manage. Is there any way to solve
this problem or to increase
I think you will find that this follows the partial
matching rules:
* exact match
* unique partial match
* take evasive action if non-unique partial match
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S User)
B Fox
The function in 'lapply' gets each component, so
write a function that does what you want given a
component of the list:
lapply(BigList, function(x) mean(x$label1))
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S
Hi there,
I'm attempting to use
color2D.matplot(x,cellcolors=cellcol, do.hex=TRUE)
from the plotrix package, but do not want to have borders around the hexagons.
I've tried various approaches such as
fg=white, col=white, lty=0, lwd=0 and border=NULL/FALSE/white, but none
of them works.
I'd
Hello R-Users!
I need a little help to build up a contingency table out of several
variables.
A-c(F,M,M,F,F,F,F,M,F,M,F,F)
B-c(0,0,0,0,0,0,1,1,1,1,0,1)
C-c(0,1,1,1,1,1,1,1,1,0,0,0)
ABC-as.data.frame(cbind(A,B,C))
ABC
A B C
1 F 0 0
2 M 0 1
3 M 0 1
4 F 0 1
5 F 0 1
6 F 0 1
7
Hello,
If I use:
p - ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram()
and then:
p + facet_grid(. ~ sc_recov)
Its a little crazy because I have 48 different values in 'sc_recov'.
Instead I want to select cases from 'sc_recov' and only use three at a time:
I tried:
p +
There are still places available on the following courses. For more
detailed information and a registration form please contact
[EMAIL PROTECTED], or visit our website at
www.mango-solutions.com http://www.mango-solutions.com/ .
R for Finance
Dieter Menne:
Thanks for the suggestion and link, it looks good. I think my trouble now is
lack of basic R knowledge.
Cheers,
Michael
For example,
I tried:
histogram(~bbContag | sc_recov, data=dat)
This would work fine if I could select cases from sc_recov. How can I select
more than one case?
Dear R Users,
Are error messages in R documented somewhere ? I am getting the following
error message and would like to track down what it means and which
specific routine triggers it:
Error: no function to return from, jumping to top level
Once I figure out which snippet of code triggers it,
I think you wrote a function, but pass only snippets of code from your script
to your R-console.
The part which return the error is the following:
return()
return can only be used to return from a function and not as a command in
the R-console.
Bart
TU wrote:
Dear R Users,
Are error
Michael Just mgjust at gmail.com writes:
sc_recov.21 - dat[dat$sc_recov=21,]
Be careful: it should be ==, not = when want a comparison.
df = data.frame(grp = letters[1:5],val=rnorm(100))
dfab = df[df$grp==a | df$grp==b,] # The simple way with or
dfcde = df[df$grp %in% c(c,d,e),] # The
Michael,
Use %in% to select multiple cases:
dat[dat$sc_recov %in% c(21, 31, 41), ]
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie
Hello everybody,
I'm a new user and I'm trying to use read.table with a tab delimited
file but the system tells me (i) that there are more columns than
column names, and (ii) that headers and col.names are of different
lengths. I have already checked my variable names and there's no blank
Hi there,
How are you? I'm having some issues with the following set up:
WINDOWS XP
R-2.7.2
RBloomberg-0.1-10
RDCOMClient-0.92-0
zoo-1.5-4
chron-2.3-24
Bloomberg is installed in the machine. The following code crashes R:
library(RBloomberg)
.bbfields - blpReadFields('C:\\path\\to\\blp\\API')
Stephen, I think this is what you want:
diag(as.matrix(d)[-1, ])
the '-1' above means that you ignore row one in the matrix. Given that
it's a square matrix, I suppose you could also also do:
diag(as.matrix(d)[-1,-19])
the '-19' above mean you ignore the 19th column of hte matrix. Both
give
First of all your construction of ABC leads to a structure with 3 factor
variables due to the way cbind processes the input variables - which is
not intended I think.
You can do sth like
ABC-data.frame(A,B,C)
aggregate(ABC[,2:3],by=list(A),sum)
hth.
Birgitle schrieb:
Hello R-Users!
I
Try this:
rapply(BigList, mean)[c(1, 3)]
On Wed, Oct 1, 2008 at 6:29 PM, eric lee [EMAIL PROTECTED] wrote:
Hi. I have a list where each object in the list has multiple parts. I'd
like to take the mean of just one part of each object. Is it possible to do
this with lapply? If not, can you
Thanks for your answer.
It is intended, that the variables are treated as class factor, because
these are binary variables with, for example, the presence or the absence of
a plant organ.
As far as I understood, I have to treat them for other calculations as
factor. Therefore I classified these
you can see the different methods of mean by: methods(mean)
Then you can type mean.default and you will see the complete code.
Bart
Alphonse Monkamg wrote:
Dear ALL,
Does anyone know how to get the complete code program for any build-in
function
in R, e.g. when I tape mean in
On Thu, Oct 2, 2008 at 1:34 PM, Alphonse Monkamg [EMAIL PROTECTED] wrote:
Dear ALL,
Does anyone know how to get the complete code program for any build-in
function
in R, e.g. when I tape mean in the R-console, I get the following:
mean
function (x, ...)
UseMethod(mean)
I have drawn a 3D scatter plot :
library(mnormt)
library(scatterplot3d)
dat = cbind(rmnorm(3, rep(0,2), diag(2)), 1:3)
scatterplot3d(dat)
Now I want to do 2 things :
1 : In the Z-axis (i.e. height), I want to see only numbers 1,2,3, etc NOT,
1,1.5,2,2.5.
2. I want to add two
Joan-Josep Vallbé pep.vallbe at uab.cat writes:
I'm a new user and I'm trying to use read.table with a tab delimited
file but the system tells me (i) that there are more columns than
column names, and (ii) that headers and col.names are of different
lengths. I have already checked my
Dear R gurus,
As part of finding initial values for a much more complicated fit I want to
fit a function of the form y ~ a + bx + cx^d to fairly noisy data and have
hit some problems.
To demonstrate the specific R-related problem, here is an idealised data
set, smaller and better fitting than
Ok, then treat them as factors - but if they are really binary and coded
0 and 1, which kind of calculation would lead to different results for a
factor instead of a numeric variable?
Anyway,
ABC-as.data.frame(cbind(A,B,C))
aggregate(ABC[,2:3],by=list(A),FUN=function(x)sum(x=='1')) # '1' is
I think it makes a difference if I want to use a classification method like
rpart () or if I use a modelling approach like glm().
Many thanks for the kind and fast help. I am still very untrained and it is
difficult for me to create such codes.
B.
Eik Vettorazzi wrote:
Ok, then treat them
Dear ALL,
Does anyone know how to get the complete code program for any build-in function
in R, e.g. when I tape mean in the R-console, I get the following:
 mean
function (x, ...)
UseMethod(mean)
environment: namespace:base
but I need the full mean function.
Thank in advance,
Hi, I've tried it with read.delim:
ayuso - read.delim(/Users/user/Rwork/ayusotriada.txt)
and now the error message is Error with make.names(col.names, unique
= TRUE) :
string multibyte 187 invalid
thanks!
On 02/10/2008, at 12:21, Joan-Josep Vallbé wrote:
Hello everybody,
I'm a new
Look at ?count.fields to diagnose this.
On Thu, Oct 2, 2008 at 8:05 AM, Joan-Josep Vallbé [EMAIL PROTECTED] wrote:
Hi, I've tried it with read.delim:
ayuso - read.delim(/Users/user/Rwork/ayusotriada.txt)
and now the error message is Error with make.names(col.names, unique =
TRUE) :
string
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Trying to work out a model formula that will do what I want ... suppose
I want to model
y = b_i x + epsilon
(i.e. a linear model with zero intercept and with slopes differing
by groups), and I want to parameterize the slopes
in the usual way of
On Thu, Oct 2, 2008 at 6:34 AM, Alphonse Monkamg [EMAIL PROTECTED] wrote:
Dear ALL,
Does anyone know how to get the complete code program for any build-in
function
in R, e.g. when I tape mean in the R-console, I get the following:
mean
function (x, ...)
UseMethod(mean)
Hi, I've done ?count.fields and it helped me identify a blank space
that I had not noticed (thanks!), but when I try to upload the table
again with
ayuso - read.table(/Users/user/Rwork/ayusotriada.txt, header=TRUE,
sep=\t, na.strings=NA, dec=., strip.white=TRUE)
or with read.delim, the
In answer to your subject heading question: to suggest a patch, send an
email to the package maintainer (as listed in the DESCRIPTION file).
Your suggestion will get a better reception if you offer sample
code/data to illustrate the problem, and a nicely written patch to the
sources to fix
klaus.steenberg.larsen at risoe.dk writes:
Dear R-Help,
I have used R2.6.0 until I recently installed also R2.7.2 (see details below)
In R 2.6.0, the following script using the function 'var' (cor(stats)):
x.test - c(NA, NA, NA, NA)
var(x.test, na.rm=T)
gives the
[EMAIL PROTECTED] wrote:
Dear R Users,
Are error messages in R documented somewhere ? I am getting the following
error message and would like to track down what it means and which
specific routine triggers it:
Error: no function to return from, jumping to top level
Once I figure out which
Uwe Ligges ligges at statistik.tu-dortmund.de writes:
Christoph Scherber wrote:
Dear R users,
Assume I have three models with the following AIC values:
modelAICdf
model1 -102
model2-125
model3-112
Obviously, model2 would be preferred,
?data.frame
Lists and data.frames allow mixed components, matrices do not
--- On Thu, 10/2/08, ZHU, Justin [EMAIL PROTECTED] wrote:
From: ZHU, Justin [EMAIL PROTECTED]
Subject: [R] matrix with different type of column [SEC=UNCLASSIFIED]
To: r-help@r-project.org
Received: Thursday, October
I don't have R on this machine to check but are you sure that sc_recov still
exists in s.all?
Otherwise do str(s.all) and see what the object looks like.
--- On Wed, 10/1/08, Michael Just [EMAIL PROTECTED] wrote:
From: Michael Just [EMAIL PROTECTED]
Subject: [R] Aggregate issues with
Did you try skipping the header as one respondent mentioned?
read.table(..., skip = 1, header = FALSE)
On Thu, Oct 2, 2008 at 9:05 AM, Joan-Josep Vallbé [EMAIL PROTECTED] wrote:
Hi, I've done ?count.fields and it helped me identify a blank space that I
had not noticed (thanks!), but when I try
Duncan Murdoch wrote:
[EMAIL PROTECTED] wrote:
Dear R Users,
Are error messages in R documented somewhere ? I am getting the
following error message and would like to track down what it means
and which specific routine triggers it:
Error: no function to return from, jumping to top level
Many thanks, that is very helpful Peter.
Regards,
Tolga
Peter Dalgaard [EMAIL PROTECTED]
02/10/2008 14:24
To
Duncan Murdoch [EMAIL PROTECTED]
cc
[EMAIL PROTECTED], r-help@r-project.org
Subject
Re: [R] error message documentation: Error: no function to return from,
jumping to top level
Hi All,
I seem to be having a few troubles with aggregating data back onto the the
dataframe,
I want to take the max value of a user, and then apply this max value back
against all id's that match (i.e a one to many matching)
Can anyone offer any advice? is there a better way of doing this?
Dummy
Sorry, I didn't, and removing the first row it effectively works.
Thank you very much for your help.
On 02/10/2008, at 15:20, Gabor Grothendieck wrote:
Did you try skipping the header as one respondent mentioned?
read.table(..., skip = 1, header = FALSE)
On Thu, Oct 2, 2008 at 9:05 AM,
See ?ave (in package stats)
/hans
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Michael Pearmain
Sent: den 2 oktober 2008 15:28
To: r-help@r-project.org
Subject: [R] Help with aggregation
Hi All,
I seem to be having a few troubles with aggregating
I'm trying to convert a variable that is imported from CSV into a datetime,I'm
trying to use the strptime function but with no joy, can anyone offer any
advise?
i have a vector
timestamp
07-07-2008-21:03:55
07-07-2008-19:52:55
07-07-2008-19:54:14
07-07-2008-15:24:01
07-07-2008-13:39:17
Hello,
I am using R version 2.7.2 on Windows. I would like to create a script
that draws many plots and pauses before drawing a new plot, so that I
could copy each plot consecutively to clipboard (Copy as metafile from
context menu) or save it to a file (Save as metafile from context
menu).
I
I'm trying to convert a variable that is imported from CSV into a
datetime,I'm
trying to use the strptime function but with no joy, can anyone offer
any
advise?
i have a vector
timestamp
07-07-2008-21:03:55
07-07-2008-19:52:55
07-07-2008-19:54:14
07-07-2008-15:24:01
My question is how does 'plot' know to implicilty call the plot.ts (in the case
of the full exact spectrum being fed back into the inverse?
Kevin
Dieter Menne [EMAIL PROTECTED] wrote:
rkevinburton at charter.net writes:
ff - complex(length(fs))
ff[9] - fs[9]
ff[5] - fs[5]
rkevinburton at charter.net writes:
My question is how does 'plot' know to implicilty call the plot.ts (in the
case of the full exact spectrum
being fed back into the inverse?
So the title should be How does the specific incarnation of object orientation
in R work? Try, for example, section
Thanks a lot! It seems to work fine now with what you suggested. I
also tried another approach by running the following on the terminal
R -f ~dir1/prog.R
and it read in some lines in prog.R, but then the execution halted at
one point for some reason I could not figure out. Any way I could get
R
I'd like to control my namespace thoroughly, separated by task. Is
there a way, in R session, to introduce namespaces for tasks
dynamically and switch them as needed? Or, is there a combination of
load/save workspace steps which can facilitate this?
Cheers,
Alexy
Thanks to all.
I summarized (in order to thank the list) the solutions to help future workers
searching subjects like this at R help.
# Number of rows
nr = 10
# Data set
dataf =
as.data.frame(matrix(c(rnorm(nr),rnorm(nr)*2,runif(nr),sort(c(1,1,2,2,3,3,sample(1:3,nr-6,replace=TRUE,ncol=4))
Hello,
Our company has been looking at SAS enterprise guide 4 (Insightful Miner is a
similar product from Insightful, I believe) which seems to provide a nice
graphical way of displaying/managing a process or project. It is also
accessible to non-programmers.
Does any of the R GUIs provide
Yes, I could prefix everything with taskN$, but that's boring! Can I
now attach and detach tasks to have one prefix-less, since they look
like dataframes? What does proto buy us?
To respond to Uwe here as well, I want dynamic task switching and
*also* an ability to load and save task
You could have an environment for each task and place your objects for
each task in its environment.
Note that this is getting close to object oriented ideas where each
enviroment/task
is an object containing your R variables and methods and the proto package
can be used to facilitate that:
What you mean is probably that you can have different workspaces.
A Workspace is associated to the current R session and is saved/loaded
in the directory from which R is started.
Note that Namespaces (for packages) in R have a different meaning.
Best,
Uwe Ligges
Alexy Khrabrov wrote:
I'd
On 02/10/2008 11:03 AM, Alexy Khrabrov wrote:
I'd like to control my namespace thoroughly, separated by task. Is
there a way, in R session, to introduce namespaces for tasks
dynamically and switch them as needed? Or, is there a combination of
load/save workspace steps which can facilitate
On Thu, Oct 2, 2008 at 11:26 AM, Alexy Khrabrov [EMAIL PROTECTED] wrote:
Yes, I could prefix everything with taskN$, but that's boring! Can I now
attach and detach tasks to have one prefix-less, since they look like
dataframes?
Read the proto vignette and other info to get a better
Hey list,
I have a barplot with a line plot overlayed (on the 2nd y axis) and it
looks fine except that the origin (0) on the 2nd y axis (side=4) is not
lining up with the origin on the 1st y-axis (side=2, from the barplot).
Is there a way to either shift the 2nd y-axis (line plot) down or
Here is what I do:
I have a different folder/directory for each of my different projects, then I
start a new R session for each project that I am working on. If you start R
with the folder/directory as the current working directory then the workspace
for that project is loaded (and the script
The updateusr function in the TeachingDemos package can be used to change the
user coordinates of the current plot and may do what you want.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
801.408.8111
-Original Message-
From: [EMAIL
Hopefully this will not seem too ignorant of a question. I am having a hard
time picking out the sources of the differences between:
abs(fft(x))^2/length(x)
and
spec.pgram(x, taper=0, log=no, plot=FALSE)
Also from the limited testing that I have done since the DC frequency is not
returned
The trick about setting the windows title is quite nifty and I've just
added that to my profile.
Just wanted to add two additional ways to start up R in the current directory
on Windows. They are provided in the batchfiles distribution. See home
page at:
http://batchfiles.googlecode.com
1.
I was wondering if I could get some advice on how best to diagnose a problem
that I am having with StructTS. To illlustrate the problem I have generated a
rigged time series:
genseq - function()
{
x - numeric(4*365)
s - seq(as.Date(2005-01-01), as.Date(2008-12-31), by=month)
ob -
help(%in%) may be of use. However I don't see any overlap in your example.
--- On Wed, 10/1/08, [Ricardo Rodriguez] Your XEN ICT Team [EMAIL PROTECTED]
wrote:
From: [Ricardo Rodriguez] Your XEN ICT Team [EMAIL PROTECTED]
Subject: [R] time segments intersection
To: r-help@r-project.org
Cc:
y-c(131.79, 131.79, 135.02, 135.55, 136.46, 136.83, 137.82, 138.00,
138.06, 138.04, 140.04, 142.44, 145.47, 144.34, 146.30, 147.54, 147.80)
x-c(194.5, 194.3, 197.9, 198.4, 199.4, 199.9, 200.9, 201.1, 201.4, 201.3,
203.6, 204.6, 209.5,208.6, 210.7, 211.9, 212.2)
fitted.results-lm(y~x)
Dear all,
I used the glm.nb with the default values from the MASS package to run a
negative binomial regression. Here is a simple example:
set.seed(123)
y - c( rep(0, 30), rpois(70, lambda=2) )
smoke - factor( sample( c(NO, YES), 100, replace=T ) )
height - c( rnorm(30,
Sorry, I must be looking at a different section but when I look at 3.4 in
r-intro.pdf I see:
3.4 The class of an object
All objects in R have a class, reported by the function class. For simple
vectors this is just the
mode, for example numeric, logical, character or list, but matrix,
array,
Hi there
I have 386 list of vectors.In that i want to delete 23 vector in the 53 list
element.
[[1]]
abc dfg dft err
[[2]]
elm erg trr rtt
similarly i have [[386]]
now i want to delete or avoid trr in the [[2]].
i just tried something like this
list[-list[[53]][3]]
--
View this message in
What version are you running? Under what OS. I don't seem to see the problem
with 2.7.2 under Windows XP.
--- On Thu, 10/2/08, Jordan Wong [EMAIL PROTECTED] wrote:
From: Jordan Wong [EMAIL PROTECTED]
Subject: [R] missing output in summary() and anova()
To: r-help@r-project.org
Received:
Adaikalavan Ramasamy a.ramasamy at imperial.ac.uk writes:
Dear all,
I used the glm.nb with the default values from the MASS package to run a
negative binomial regression. Here is a simple example:
[snip -- thanks for the example!]
The question now is how do I report the results, say, for
It is 2.7.2 under Windows XP SP3. I have tried 2.6.1,2.6.2 in another
computer under Windows XP SP3 also and the problem still occur.
2008/10/3 John Kane [EMAIL PROTECTED]
What version are you running? Under what OS. I don't seem to see the
problem with 2.7.2 under Windows XP.
--- On Thu,
Adaikalavan Ramasamy wrote:
Ben, fantastic. Thank you for confirming it.
One more question. What do you call the back transformed variable?
In my domain, people use something called the ratio of mean but I am not
sure if it is the same. I am not what the ratio is between.
It depends.
Try this:
library(zoo)
# create two time series to test with
z1 - zooreg(0:10, Sys.Date())
z2 - lag(z1, 5)
z1
2008-10-02 2008-10-03 2008-10-04 2008-10-05 2008-10-06 2008-10-07
2008-10-08 2008-10-09 2008-10-10 2008-10-11 2008-10-12
0 1 2 3 4
imho better: ?xyplot using the groups argument.
-- Bert
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Kingsford Jones
Sent: Thursday, October 02, 2008 11:17 AM
To: Michel PETITJEAN
Cc: r-help@r-project.org
Subject: Re: [R] 3 curves / 1 plot
Hello R-users
Kindly assist me with the following plot problem in R-
My data looks like this (shortened because I have about 5000 observations).
dat - read.table(textConnection(Id Time Y
1 0 194
1 5.22 179
1 5.97 190
2 1.61 265
2 2.1 234
2 16.4 300
2 29.5 345
3 0 212
3 0.36
Michel PETITJEAN-2 wrote:
I am a new user of R.
Please does somebody knows how to plot 3 datasets
(x1,a1),...,(xn,an), (x1,b1),...,(xn,bn), and (x1,c1),...,(xn,cn)
on a single x,y plot, each of the three datasets being plotted with
its own character pch() ?
(three calls to plot() erase
I have a list, myList, with each of its 9 components being a 15X15
matrix. I want to run a t-test across the list for each component in
the matrix. For example, the first t-test is on myList[[1]][1, 1],
myList[[2]][1, 1], ..., myList[[9]][1, 1]; and there are totally 15X15
t-tests. How can I run
I got following error :
library(mnormt)
library(scatterplot3d)
dat = cbind(rmnorm(3, rep(0,2), diag(2)), 1:3)
s3d - scatterplot3d(dat, lab.z=2, scale.y=0.7, angle=20)
s3d$plane3d(2,0,0, solid, col=grey)
Error in object$coefficients : $ operator is invalid for atomic vectors
I am using
I want to run a R program, prog.R, interactively. My question is, is
there a way I can start prog.R on the shell terminal when invoking R,
instead of using source() inside R?
TIA,
Gang
__
R-help@r-project.org mailing list
The default range is the range of the variable in the data. It wouldn't
make much sense to extend beyond that. If you want to restrict the
range, the varImpPlot() function does return the x-y pairs, so you can
construct your own plot from that.
Andy
From: A Van Dyke
i am new to R and using
Hello there,
I have a fairly simple request (I hope!)
I have produced a filled contour plot like this:
library(grDevices)
library(gplots)
library(plotrix)
filled.contour(contour, axes=F, frame.plot=TRUE, color=terrain.colors, ylab=
Length Along Flume (m), key.title =
Hi Shiva,
Maybe you are interested in the following paper:
Learning when Training Data are Costly: The Effect of Class Distribution
on Tree Induction. G. Weiss and F. Provost. Journal of Artificial
Intelligence Research 19 (2003) 315-354.
For validating the models in those enviroments,
I am sure you will get helpful answers. I am almost as sure that you
shouldn't be doing this. I suggest you consult with your local statistician.
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Gang Chen
Sent: Thursday, October 02, 2008
Is this what you want to do - it was not clear from your example:
x
[[1]]
[1] 1
[[2]]
[1] 2
[[3]]
[1] 3
[[4]]
[1] 4
x - x[-3]
x
[[1]]
[1] 1
[[2]]
[1] 2
[[3]]
[1] 4
Or maybe this:
x - list(c( 1,2 ,3,4),c(4,5,6,7),c(8,9))
x
[[1]]
[1] 1 2 3 4
[[2]]
[1] 4 5 6 7
[[3]]
[1] 8 9
x[[2]]
I appreciate your suggestion. The example I provided was fabricated
because I was only focusing on the programming perspective on how to
deal with such a data structure, not real statistical issues. Do you
mind elaborating a little more why that would not be appreciate? I
know I can do it with a
Hi,
To pretty print aggregates by various dimensions I needed to add a
empty row in output of aggregate.
For example.
d-(aggregate(data[,cbind(x)], by=list(data$group1,data$group2),
sum))
Group.1 Group.2 x
1 A N 3
2 A Y
Dear all
When I used the method, L-BFGS-B, in OPTIM, I've got the following message.
-
$par
[1] 0.176166426835580
$value
[1] 1322.17600079332
$counts
function gradient
88
$convergence
[1] 0
$message
[1]
Hi!
Yes, the example in the help of pvals.fnc is working.
Is something wrong with my lmer-object? It´s created with
m-lmer(file2[,11] ~ file2[, 1] + file2[, 2] + file2[, 3]+ file2[, 4] +
file2[, 5] + file2[, 6]+ file2[, 7]+ file2[, 8] + file2[, 9] + file2[,
10] + (file2[, 5] | file2[, 2]),
Yes, the estimates are reliable when you get that messsage. If you want a
more rigorous convergence, you can decrease factr, whose default is 1e07.
However, I am not sure why you need this.
Read the help for `optim'
?optim
Ravi.
Hi John,
John Kane wrote:
help(%in%) may be of use. However I don't see any overlap in your example.
You are right! I've choosen an example where the intersection is 0. But we
could consider this other case...
endPeriod-as.POSIXlt(2008-09-30)
startPeriod-as.POSIXlt(2007-10-01)
1 - 100 of 115 matches
Mail list logo
