Thanks Gabor,
Gabor Grothendieck wrote:
Try:
chooseCRANmirror()
It works the call to the pop-up window, but it fails now whatever mirror
I use with the following message:
> install.packages("zoo",repos=chooseCRANmirror())
/bin/sh: tar: command not found
Error in sprintf(gettext(fmt, domain
On 10/2/08, Michael Just <[EMAIL PROTECTED]> wrote:
> Dieter and Thierry:
>
> Per you suggestions I have tried:
>
> ggplot2 from Thierry:
>
> > p <- ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram()
>
> > p + facet_grid(. ~ sc_recov %in% c(21,31,41))
>
> But get the followinng error
> idx = c(1,3,5)
> as.formula(paste("Surv(time, status)~", paste("x", idx, sep = "", collapse =
> "+")))
Surv(time, status) ~ x1 + x3 + x5
Regards,
Yihui
--
Yihui Xie <[EMAIL PROTECTED]>
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.yihui.name
Scho
Dear R user,
I got a question when using the coxph function. I have 5 covariates x1, x2,
x3, x4, x5 and I want to write a function so that when given an indicator,
e.g., c(1,3,5), I can fit the model as
model=coxph(Surv(time, status)~x1+x3+x5). Any idea to play around the form
of formula?
Thank y
Dieter and Thierry:
Per you suggestions I have tried:
ggplot2 from Thierry:
> p <- ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram()
> p + facet_grid(. ~ sc_recov %in% c(21,31,41))
But get the followinng error:
Error in check_formula(formula, varnames) :
Formula contains variables
Hi Duncan, Mark and all,
Thanks for the suggestion. I ve tried the below suggestion.
I got this error now
Error in `[<-.data.frame`(`*tmp*`, , x, value = NULL) :
new columns would leave holes after existing columns
Its weird as I try to assign NULL manually and its not complaining anything.
W
On 02/10/2008 10:07 PM, Jason Lee wrote:
Hi,
I came across the below error when I try to do ifelse condition:-
Error in "[<-"(`*tmp*`, test, value = rep(yes, length = length(ans))[test])
:
incompatible types
What I am trying to accomplish is :-
for(x in 1:ncol(filterpred)){
sumno<-su
Hi,
Let say I have a normal density X~n(0,1) and I have a line y=0.01x+0.07. the
following code generate the plots.
x=seq(-10,10,length=100)
plot(x,p1,type='n',ylab="Density",main="Overlap Measure",xaxt="n",yaxt="n")
pi=dnorm(x,0,1)
points(x,p1,type='l')
abline(0.07,0.01)
you can see that the
You can also use the 'reshape' package:
> d.m <- melt(dat, measure.var='x')
> d.m
Group1 Group2 variable value
1 A Yx 1
2 B Nx 1
3 A Yx 1
4 B Nx 420164904
5 A Nx 3
> cas
Here is one way of doing it:
> dat <- read.table(textConnection("Group1 Group2 x
+ AY1
+ BN1
+ AY1
+ BN 420164904
+ AN 3"), header=TRUE, as.is=TRUE)
> closeAllConnections()
> d <- a
Hi,
I came across the below error when I try to do ifelse condition:-
Error in "[<-"(`*tmp*`, test, value = rep(yes, length = length(ans))[test])
:
incompatible types
What I am trying to accomplish is :-
for(x in 1:ncol(filterpred)){
sumno<-sum(filterpred[no,x])
sumyes<-sum(filterpred[
My dear R buddies,
I've run into a problem when doing numerical computation recently. In
a program that I've been working on, I usually get a vector of real
values which are theoretically (and it's correct) supposed to decrease
until reaching zero after a given value. However, most of the value
ju
I try to use RBloomberg to get the dividend for IBM. However,
blpGetData(conn, "IBM EQUITY", field="EQY_DVD_HIST_ALL",
start=as.chron("1980-01-01"))
doesn't work. It returns
EQY_DVD_HIST_ALL
(10/02/08 14:46:36) NA
I have to used
blpGetData(conn, "IBM EQUITY",
Try:
chooseCRANmirror()
On Thu, Oct 2, 2008 at 7:47 PM, [Ricardo Rodriguez] Your XEN ICT Team
<[EMAIL PROTECTED]> wrote:
> OK. Please, how do I force the select directory window to pop-up after
> selecting a given repos and without restarting R?
>
> When I submit install.packages("zoo", dep = TRU
Not yet. I'll keep trying.
Please, what format string must I pass to strptime to get the date as a
number? Thanks.
Cheers,
Ricardo
[EMAIL PROTECTED] wrote:
did you get zoo to work and try gabor's solution because i wasn't
really following what you wanted.
On Thu, Oct 2, 2008 at 7:25 P
OK. Please, how do I force the select directory window to pop-up after
selecting a given repos and without restarting R?
When I submit install.packages("zoo", dep = TRUE) for the first time
after R starts, I am required to select a mirror. But this doesn't
happen in the next install command at
Select a different mirror until you find one that works for you.
On Thu, Oct 2, 2008 at 7:38 PM, [Ricardo Rodriguez] Your XEN ICT Team
<[EMAIL PROTECTED]> wrote:
> Hi,
>
> Gabor Grothendieck wrote:
>>
>> Try this:
>>
>> install.packages("zoo", dep = TRUE)
>
> Here the result. I've tried with a cou
Hi,
Gabor Grothendieck wrote:
Try this:
install.packages("zoo", dep = TRUE)
Here the result. I've tried with a couple of mirrors (Bristol and Spain):
> install.packages("zoo", dep = TRUE)
--- Please select a CRAN mirror for use in this session ---
Warning: unable to access index for reposit
Hi Mark,
[EMAIL PROTECTED] wrote:
hi: I don't know if it's the best way but one way is to use an
POSIXct object. ( you can convert using as.POSIXct(
strptime(datestring)) because then each unit of digit is a second
and you can proceed from there. if you explain what you want to do
exactly,
Try this:
install.packages("zoo", dep = TRUE)
On Thu, Oct 2, 2008 at 7:15 PM, [Ricardo Rodriguez] Your XEN ICT Team
<[EMAIL PROTECTED]> wrote:
> Thanks Gabor,
>
> I am afraid I have some problems installing zoo package. Please, could you
> take a look at this?
>
>> install.packages("packagename"
Thanks Gabor,
I am afraid I have some problems installing zoo package. Please, could
you take a look at this?
> install.packages("packagename",repos="http://R-Forge.R-project.org";)
Warning message:
package 'packagename' is not available
> install.packages("zoo",repos="http://R-Forge.R-project
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")
endP
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!
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]), fi
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] "CON
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
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 co
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
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 11
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,
Willia
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 = title(main="Velocity\n(m/
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 usin
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
https://stat.ethz.ch/mailman/l
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 usi
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 th
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()
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
example(
example(matplot)
or
?points
On Thu, Oct 2, 2008 at 11:43 AM, Michel PETITJEAN
<[EMAIL PROTECTED]> 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 datas
Megh Dal wrote:
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 wa
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 the two first datasets).
Thank you ver
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
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 depe
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.
>
>
> ---
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.
Regards, Adai
Ben Bolker wrote:
Adaikalava
Adaikalavan Ramasamy 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 height
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
> Rec
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
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",
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, me
> 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)
> sum
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.o
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 <-
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
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
return
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: [EMAI
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
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
readj
Wonderful, great! That solves my problem nicely (on to the next one!).
Thanks a lot,
Keith Jewell
"Katharine Mullen" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
> a) you should put values for Ca, Cb, Cc directly into the data list as
> data=list(Ca=1,
>
> b) you can simplify
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 unders
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
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 like
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:
libra
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 from
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 s
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))
n
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 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
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 3.4 and
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:
> charter.net> writes:
>
> >
> > ff <- complex(length(fs))
> > ff[9] <- fs[9]
> > ff[5] <- fs[5]
>
i am new to R and using the random forest package. is there a way to specify
the x-axis scale range for the variable importance plot? many thanks.
-alison
--
View this message in context:
http://www.nabble.com/specifying-x-axis-scale-on-random-forest-variable-importance-plot-tp19780560p197805
a) you should put values for Ca, Cb, Cc directly into the data list as
data=list(Ca=1,
b) you can simplify the call to
# idealised data set
aDF <- data.frame( x= c(1.80, 9.27, 6.48, 2.61, 9.86, 5.93, 6.76, 5.52,
6.06, 8.62),
y= c(24.77, 2775.07, 895.15, 60.73, 3373.57, 677.82, 1021.92, 542.8
> 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'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-2
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
07-07-2008-13:39:17
07-07-2008-
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 data
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, Joan-
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
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"
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 t
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 t
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
?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, O
Uwe Ligges 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
[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 s
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 output:
>
> NA
>
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 it
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.deli
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("
-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 h
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)
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
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' i
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 n
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 tha
Joan-Josep Vallbé 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 variable na
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 Horizo
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("
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
1 - 100 of 127 matches
Mail list logo
