You want to use list:
a = list()
a[[1]] = c(1,2,3)
a[[2]] = c(1,2,3)
Note the [[..]] operator - check the "An Introduction to R" manual for
more details
Schalk Heunis
On Wed, Sep 9, 2009 at 6:00 AM, wrote:
> I just learned that vectors can't contain vectors, which frankly simply
> confuses m
I just learned that vectors can't contain vectors, which frankly simply
confuses me (can you imagine arrays or lists in any other language not being
able to contain arrays or lists?). At any rate, I need to create a data
structure (size to be determined at runtime) which I will instantiate and t
As you see in the data.. each section has table header followed by a row of
table column headers and then rows of data. I felt read.table may not be the
optimal function to extract both the texts (headers) and numerics
(rows&columns of data). As of now, I use "scan" to get the text boundaries
and u
Try
?Rscript
?system
?shell
On Tue, Sep 8, 2009 at 11:24 PM, Daniel Bradley wrote:
> Hi!
>
> I have a need to do the ridiculous and run a batch file if certain
> conditions are met within an R script. is there a simple way to do so? or
> to run dos commands?
>
> Thanks!
> Daniel
>
> [[alt
You would normally have R or SAS run on the server, not on the client machine so
whether it opens a window or not makes no difference. The window is not on the
client's machine anyways. The client only has a browser.
>From within R see ?Rscript and ?system
On Tue, Sep 8, 2009 at 11:19 PM, Carlo
Try this:
Lines <- "Date ABC D E
1978-10-22 18 20.64 0.0 0.176 -1.76
1978-10-23 15 17.06 0.4 0.147 2.52
1978-10-243 7.588 0.0 0.068 -6.86
1978-10-259 11.491 0.0 0.102-1.01
1978-10-26 13 14.98 1.4
Carlos,
The last time I checked the SAS license for web use was $25,000 per
year. R costs $0 per century. Still interested in SAS?
Not sure about your question about 'DOS window'. R on a web server
interacts with your browser and doesn't open other windows. There are
several online demos
Hi!
I have a need to do the ridiculous and run a batch file if certain
conditions are met within an R script. is there a simple way to do so? or
to run dos commands?
Thanks!
Daniel
[[alternative HTML version deleted]]
__
R-help@r-project.org
Hello Everyone,I think this is a very simple problem, I have been struggling
with it for a few days now.
I have a 10-year daily data in the following format.
Date ABC D E
1978-10-22 18 20.64 0.0 0.176 -1.76
1978-10-23 15 17.06 0.4 0.14
Good evening,
I have been asked to investigate the pros and cons of using SAS vs. R in a web
application. Either SAS or R would be the engine used to make some very simple
calculations and to produce graphs, preferably in png format.
The advantages of R are pretty obvious as there would be no l
Ben Ridenhour wrote:
>
> Hello all,
> I'm trying to establish some confidence intervals on predictions I am
> making using
>
>>predict(nls(...))
>
> and predict.nls (unfortunately) does not utilize the se.fit option. A
> little more background is that I am trying to match the output with old
Or try the following:
require(rms)
data(lung)
survfit(Surv (time, status) ~ ph.ecog, data=lung)
lung <- transform(lung, ph.ecog=factor(ph.ecog))
f <- cph(Surv(time, status) ~ strat(ph.ecog), data=lung, surv=TRUE)
## doesn't assume PH
quant <- Quantile(f) # print quant to see how all quantiles a
Szumiloski, John wrote:
>
> Dear useRs,
>
> I have a complicated function to be optimized with optim(), and whose
> parameters are passed to another function within its evaluation. This
> function allows for the parameters to enter as arguments to various
> probability distribution functions.
Hi,
I have a strong recollection of a relatively recent communication on
this issue with Terry, but cannot seem to locate it at the moment and
I don't recall if it was on-list or off-list. I am either on drugs or
Terry indicated that he was going to modify either the returned
survfit() or
How do you want to extract the data? You can use 'readLines' to read
in the data and then 'grep' to find the header lines and delete them.
On the other hand, do you want to separate each section into a
differnet object/list? You can again use readLines and determine
where the breaks are and then
Hi
Out of curiosity, how can i determine the minimum required sample size of this
data in a combined form say 100 entries which r functions can i use?
R helpers
Please help me combine the simulated data to a form of table where: Hypermarket
have 10 rows, supermarket have 15 rows,...,
Suggest you try writing this in C.
On Tue, Sep 8, 2009 at 8:00 PM, R_help Help wrote:
> Hi Gabor,
>
> I appreciate your suggestion. I'm wondering if there's any faster
> implementation that one could achieve. The dataset I have is about
> 400K rows and I have many of them. Just wondering if you ha
First, dependent on how highly dimensional your data is, I would recommend
inspecting it visually first. Does the step-function seem a reasonable
assumption?
If I understand you correctly when you say, "piecewise linear which has the
gradient constrained to zero," you mean dummy-variable coding f
Try this:
matplot(matrix(1:length(d), NCOL(d), byrow = TRUE), t(d), type = 'h', lty =
1, lwd = 2)
On Tue, Sep 8, 2009 at 11:44 AM, Rafael Moral
wrote:
> Dear useRs,
>
> I want to plot the following barplot with lines instead of bars. Is there a
> way?
>
> data <- data.frame(cbind(k = 0:3, fk = c
Michael Chajewski fordham.edu> writes:
>
> Dear All,
>
> Has anyone tried to use this plug-in? Since I am running R-2.9.1 it will not
even let me install it. Further,
> since I am running Windows I cannot use the R provided R-2.7.0 Linux
installation file from the archive
> (tried to install
R helpers
Please help me combine the simulated data to a form of table where: Hypermarket
have 10 rows, supermarket have 15 rows,..., spazashops with 35 rows.
Hypermarket <- rnorm(10, mean=2, sd=7000)
Supermarket <- rnorm(15, mean=12000, sd=4000)
Minimarket <- rnorm(20, mean=1, s
I don't understand the cbind(bi) sintax, but you can do this with the
folowing:
(Using iris data from R)
form <- formula(paste(paste(names(iris), collapse = " + "), "~ 1"))
2009/9/8 "Biedermann, Jürgen"
> Hi there,
>
> I have the following problem:
>
> I have a package called "polLCA" which has
Hi Gabor,
I appreciate your suggestion. I'm wondering if there's any faster
implementation that one could achieve. The dataset I have is about
400K rows and I have many of them. Just wondering if you have any
thoughts. Thanks.
Sakda
On Tue, Sep 8, 2009 at 6:57 AM, R_help Help wrote:
> --
Hi
Thanks for the further analysis on this Ted. I think the problem is
that, with such a "wide" plot, you are running into the default paper
size. If you look at the EPS produced by ghostscript, you will see a
line like this ...
612 792 /letter setpagesize
... and notice that the value 61
The survfit.object help page says:
"The print.survfit method does more computation than is typical for a
print method and is documented on a separate page."
It takes a bit of digging, but after first trying:
getAnywhere(print.survfit) # and then following code and trying
getAnywhere(survme
njhuang86 wrote:
Hi all,
Does anyone know how to plot overlapping qqnorm plots on the same
window?
Suppose I have data in the vector x and y:
qqnorm(x)
lines(qqnorm(y))
I though these two lines will do the job... However, lines doesn't
seem to
work. Anyways, thanks in advance!
E.g.:
You can do this in lattice:
library(lattice)
df = data.frame(val = rnorm(100), group = rep(c('x', 'y'), each = 50))
qqmath(~val, groups = group, data = df)
Hope that helps,
Greg
njhuang86 wrote:
Hi all,
Does anyone know how to plot overlapping qqnorm plots on the same window?
Suppose I have
Hi,
I'm sure this should be simple but I can't figure it out! I want to get the
median survival calculated by the survfit function and use the value rather
than just be able to print it. Something like this:
library(survival)
data(lung)
lung.byPS = survfit(Surv (time, status) ~ ph.ecog, data=
Have a look at ?segments and/or ?arrows.
--- On Tue, 9/8/09, rafamoral wrote:
> From: rafamoral
> Subject: Re: [R] barplot with lines instead of bars
> To: r-help@r-project.org
> Received: Tuesday, September 8, 2009, 2:12 PM
>
> How can I draw thin bars in a barplot?
> Rafael
>
>
> hadley w
--- On Tue, 9/8/09, hadley wickham wrote:
> From: hadley wickham
> Subject: Re: [R] barplot with lines instead of bars
> To: "rafamoral"
> What's the difference between a line
> and a thin bar?
> Hadley
Diet?
>
> On Tue, Sep 8, 2009 at 12:17 PM, rafamoral
> wrote:
> >
> > I'm sorry, but
Hello,
Maybe you want something like this:
desD2 <- optFederov(~(Intro+Duration+GOTO+Fee+Color)^2,dat, nTrials =
30)
In any case, both SAS's proc optex and R's optFederov implement a
non-exhaustive search algorithm and nothing guarantees that the final
design will be the same.
However, you can
Dec (6/5) 2006 Newsletter mentioned that EBImage might include image
registration/alignment as "of foreseeable interest".
My present Newsletter and Help Digest search finds no mention of this. Does
anyone know of someone who might have implemented an image alignment or
image registration algorit
On Tue, 8 Sep 2009, Welma Pereira wrote:
Hello, I have the following makefile. The problem is that the bibliography
doesn t work. Any help would be appreciated! I really don t don t what to
do..:-(
# The sources of the report (tex, Rnw and other files (e.g. bib, idx))
TEX_CMPS = Report proble
Petar Milin wrote:
I would like to have a line on this plot, instead of two points:
x1 = rnorm(100, 10, 2.5)
x2 = rnorm(100, 26, 3.2)
x1 = as.data.frame(x1)
x2 = as.data.frame(x2)
colnames(x1) = 'anxiety'
colnames(x2) = 'anxiety'
x1$gender = 'male'
x2$gender = 'female'
dat = rbind(x1, x2)
requi
pietro bulian wrote:
I have two alternative Cox models with C-statistics 0.72 and 0.78. My question
is if 0.78 is significantly greater than 0.72. I'm using rcorrp.cens. I cannot
find the U statistics in the output of the function. This is the output of the
help example:
x1 <- rnorm(400)
x2
I have two alternative Cox models with C-statistics 0.72 and 0.78. My question
is if 0.78 is significantly greater than 0.72. I'm using rcorrp.cens. I cannot
find the U statistics in the output of the function. This is the output of the
help example:
> x1 <- rnorm(400)
> x2 <- x1 + rnorm(400)
>
On Tue, 08-Sep-2009 at 06:40AM -0400, Duncan Murdoch wrote:
>
> I don't see the underscores in that posting, but I do see this as the
> last line in the headers:
That's interesting. I see them using Mutt as well as with Thunderbird.
I bounced the message to a Yahoo account, and the underscore
Hi all,
Does anyone know how to plot overlapping qqnorm plots on the same window?
Suppose I have data in the vector x and y:
qqnorm(x)
lines(qqnorm(y))
I though these two lines will do the job... However, lines doesn't seem to
work. Anyways, thanks in advance!
--
View this message in context:
Hi all--
I'm new to R, statistics and programming, so sorry if this is a really basic
question!
I have plotted a directional variogram, and I want to
a. overlay the omni-directional line over each directional panel
b. display the directional variograms in a single panel with a legend that
ass
On 8/09/2009, at 9:07 PM, FMH wrote:
Dear All,
I'm looking for a way on computing the derivative of first and
second order of a smoothing curve produced by a nonprametric
regression. For instance, if we run the R script below, a smooth
nonparametric regression curve is produced.
provid
On Tue, Sep 8, 2009 at 7:05 PM, Lauri Nikkinen wrote:
> Well, yeah, Henrique's solutions works fine with this data. Thanks for
> that, although this is not so generic solutions which I was looking
> after. As I originally posted, I was looking for solution which uses
> the field-width information,
Genentech has a full time job opening for a statistician in its nonclinical
biostatistics group in South San Francisco, CA. For details, please go to
the R-sig-jobs mailing list here:
https://stat.ethz.ch/pipermail/r-sig-jobs/2009q3/70.html
(you may need to subscribe to view this).
Bert Gu
Dear All,
I am attempting to use R to perform an ANOVA with three factors:
feature (3 levels), group (5 levels), and patient (246 levels), where
patient is nested within group.
Currently I am using the "lm" function to fit the model, with the
following form:
fit <- lm(intensity ~ featu
Hello all,
I'm trying to establish some confidence intervals on predictions I am making
using
>predict(nls(...))
and predict.nls (unfortunately) does not utilize the se.fit option. A little
more background is that I am trying to match the output with older SAS routines
to maintain consistency
How can I draw thin bars in a barplot?
Rafael
hadley wrote:
>
> What's the difference between a line and a thin bar?
> Hadley
>
> On Tue, Sep 8, 2009 at 12:17 PM, rafamoral
> wrote:
>>
>> I'm sorry, but I think I was misunderstood. What I need is something like
>> this:
>>
>> http://img525.ima
Well, yeah, Henrique's solutions works fine with this data. Thanks for
that, although this is not so generic solutions which I was looking
after. As I originally posted, I was looking for solution which uses
the field-width information, as Ted pointed out. But as I already
mentioned, it seems that
Hi,
I am using ada to predict a data set with 36 variables
ada(x~.,data=train,iter=Iter,
control=rpart.control(maxdepth=4,cp=-1,minsplit=0,xval=0))
can any one tell me in in laymans terms
maxdepth- how do you set this, how do you change this to improve
predictions success
cp- same question fo
What's the difference between a line and a thin bar?
Hadley
On Tue, Sep 8, 2009 at 12:17 PM, rafamoral wrote:
>
> I'm sorry, but I think I was misunderstood. What I need is something like
> this:
>
> http://img525.imageshack.us/img525/2818/imagemyu.jpg
>
> Lines instead of bars
>
> Thanks!
>
> Raf
It's all simply a matter of definitions, and there are many who disagree. See
?quantile , specifically the "type" argument. Since IQR does not appear to
have a type argument, you could easily write your own versions of these that do
what SAS does (assuming that is your goal).
With x defined
Hi James,
Take a look at the "recode" function in the "car" package. It might be
useful in this case.
HTH,
Jorge
On Tue, Sep 8, 2009 at 12:01 PM, wrote:
> Hello,
>
> I am trying to map a factor variable within a data frame to a new variable
> whose entries are derived from the content of the o
use 'ifelse'
# not tested; you supply data for the '%in%'
map <- function(x){
ifelse(x %in% c('a','b'), "North",
ifelse(x %in% c('c','d'), "South",
ifelse(x %in% c('e', 'f'), "East",
ifelse(x %in% c('g', 'h'), "West", NA
}
On Tue, Sep 8, 2009 at 12:01 PM, wrot
If your data.frame was called "test" below,
nrow(unique(test))
would do what you want, I believe.
Erik
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of "Biedermann, Jürgen"
Sent: Tuesday, September 08, 2009 9:24 AM
To: r-help@
Dear Jürgen,
If x is your data, here are two suggestions:
# Suggestion 1
length(unique( apply(x, 1, paste, sep="", collapse="") ) )
# [1] 2
# Suggestion 2
res <- as.data.frame( xtabs( ~. , as.data.frame( x ) ) )
dim(res[res$Freq > 0,])[1]
# [1] 2
HTH,
Jorge
2009/9/8 "Biedermann, Jürgen"
> Hi
try this:
> x
V1 V2 V3
1 0 3 1
2 1 2 0
3 1 2 0
> unique(x)
V1 V2 V3
1 0 3 1
2 1 2 0
> nrow(unique(x))
[1] 2
>
2009/9/8 "Biedermann, Jürgen" :
> Hi there,
>
> Does anyone know a method to calculate the number of different patterns in a
> given data frame. The variables are of p
try this:
DF <- data.frame(V1 = c(0,1,1), V2 = c(3,2,2), V2 = c(1,0,0))
DF
nrow(unique(DF))
I hope it helps.
Best,
Dimitris
Biedermann, Jürgen wrote:
Hi there,
Does anyone know a method to calculate the number of different patterns
in a given data frame. The variables are of polytomous ty
Hello, I am using "B" as a vector to store all the t.tests. since i am a
newbie to both R and statistics I am not sure if "B" is the list. Also I see
"c" used in the do.call formula and do not know what it is being used for. I
used aggregate but getting this error" Error in aggregate.data.frame(eo
Vishal gmail.com> writes:
>
> Jari, thanks for the quick answer.
>
> > sqrt(eigen(cov.trob(mydataforellipse)$cov)$values)
>
> what will this return?
>
> For my data, I get:
>
> > sqrt(eigen(cov.trob(r)$cov)$values)
> [1] 1.857733e-05 4.953181e-06
>
> Is this Left hand value the major or the
Hi Folks,
After seeing some tips on the web, I managed to get RSPython-0.7-1
installed with R 2.9.0 (on Ubuntu 8.04). I had to get rid of the extra
comma in the PythonPath function in Python.S:
.PythonPath <-
function(path=NULL, merge = TRUE)
{
dirSep <- ifelse(R.version$os == "Win32", ";", "
okay fixed it by putting "c" in quote marks.
1Rnwb wrote:
>
> Hello, I am using "B" as a vector to store all the t.tests. since i am a
> newbie to both R and statistics I am not sure if "B" is the list. Also I
> see "c" used in the do.call formula and do not know what it is being used
> for. I
Hello,
I am trying to map a factor variable within a data frame to a new variable
whose entries are derived from the content of the original variable and there
are fewer factors in the new variable. That is, I'm trying to set up a
surjection.
After first thinking that this would be a common op
Dear R users,
Something is strange in summary and IQR. Suppose, I have a data set and I
would like to find the Q1, Q2, Q3 and IQR.
x<-c(2,4,11,12,13,15,31,31,37,47)
> summary(x)
Min. 1st Qu. MedianMean 3rd Qu.Max.
2.00 11.25 14.00 20.30 31.00 47.00
> IQR(x)
[1] 19.75
SOLVED.
Thanks to a reply off-list it appears that the 'space' in "published
11" is actually some kind of multibyte character. If I physically
delete the 'space' and replace it by using the spacebar on my
keyboard, then strsplit() behaves as expected.
I had got the text from a hyperlink and copy
Hi there,
Does anyone know a method to calculate the number of different patterns
in a given data frame. The variables are of polytomous type and not
binary (for the latter i found a package called "countpattern" which
unfortunately only functions for binary variables).
V1 V2 V3
0 3
Hi there,
I have the following problem:
I have a package called "polLCA" which has the following syntax:
poLCA(formula, data)
and needs the following formula definition:
formula <- cbind(V1,V2,V3,...)
So far so good.
What I tried now was the following:
#Get "data" with the "read.table" fuct
swertie voila.fr> writes:
>
> Hello, I performed a Mantel test and plotted communitiy similarities. I
> would like to add a least square line. I thought about using abline taking
> as slope the r-statistic of the Mantel test and calculating the y-intercept
> analytically. Is this method correct?
I'm sorry, but I think I was misunderstood. What I need is something like
this:
http://img525.imageshack.us/img525/2818/imagemyu.jpg
Lines instead of bars
Thanks!
Rafael.
ONKELINX, Thierry wrote:
>
> Here is a solutions using ggplot2 and reshape
>
> library(reshape)
> library(ggplot2)
> d
Consider the following simple example using R-2.9.0 and 'perm' 2.9.1:
> require('perm')
> p<- c(15,21,26,32,39,45,52,60,70,82)
> g<- c('y','n','y','y', rep('n',6)) #Patients ranked 1,3,4 receive treatment
> permTS(p ~ g, alternative = 'two.sided', method='exact.ce') #find
p-value by complete enu
Dear useRs,
I have a complicated function to be optimized with optim(), and whose
parameters are passed to another function within its evaluation. This
function allows for the parameters to enter as arguments to various
probability distribution functions.
However, I am violating some scoping con
On 08-Sep-09 16:17:00, David Winsemius wrote:
> On Sep 8, 2009, at 12:00 PM, Lauri Nikkinen wrote:
>> Ok, I think that I have to give up and try to get this data separated
>> by some char. It seem pretty much impossible to separate those fields.
>> Thanks for your help and efforts.
>
> The solutio
Hi Robin --
Robin Hankin wrote:
> Hi
>
> I deal with long vectors almost all of whose elements are zero.
> Typically, the length will be ~5e7 with ~100 nonzero elements.
>
> I want to deal with these objects using a sort of sparse
> vector.
>
> The problem is that I want to be able to 'add' two
Hello, I have the following makefile. The problem is that the bibliography
doesn t work. Any help would be appreciated! I really don t don t what to
do..:-(
# The sources of the report (tex, Rnw and other files (e.g. bib, idx))
TEX_CMPS = Report problem
RNW_CMPS = prop1 prop2 ExeExps
OTHER = Re
Hi
Actually there is no more trouble with R 2.9.2 version. I was using 2.5.0 and
really needed to update!
Thanks anyway
Edwige
De : William Dunlap
Envoyé le : Mardi, 8 Septembre 2009, 18h12mn 59s
Objet : RE: [R] R-crash when loading workspace - Windows
Coul
On Sep 8, 2009, at 12:00 PM, Lauri Nikkinen wrote:
Ok, I think that I have to give up and try to get this data separated
by some char. It seem pretty much impossible to separate those fields.
Thanks for your help and efforts.
The solution that Henrique offered seems to be a complete one:
rea
Lauri,
Having looked at your example file, and examined its byte-by-byte
content, it is a plain ASCII file which gives exactly the same layout
as you originally posted. This, along with the field-width information
you originally supplied, is not sufficient to determine a
unique dcomposition into fi
Ok, I think that I have to give up and try to get this data separated
by some char. It seem pretty much impossible to separate those fields.
Thanks for your help and efforts.
-L
2009/9/8 Lauri Nikkinen :
> This is the file (see the attachment) that represents the problem I'm
> facing with the ori
This is the file (see the attachment) that represents the problem I'm
facing with the original file. I'm looking for some generic way to
solve this problem. Thank you for your time.
-L
2009/9/8 Barry Rowlingson :
> On Tue, Sep 8, 2009 at 1:52 PM, Lauri Nikkinen wrote:
>
>> But this is not the sol
Hi list,
I have a data frame with a "Date" column and a "Price" column - for example:
Date Price
01/01/2009 5.45
01/03/2009 6.53
01/04/2009 7.55
01/06/2009 6.76
01/08/2009 4.12
01/18/2009 5.87
...
As you can see, there are days for which I don't have any data
Win7, R-2.9.1
Dear list,
I have a large (nxm) matrix which is the output of an analysis. Since it
is so large, I would like to use output formatting to make it easier to
find particular values and patterns in the matrix. In particular, I want
to print the matrix as follows:
- cells with a va
Alex,
It's mainly speculation, as I cannot check the Excel add-in nor Vassar, but
I'll give it a try.
For the Friedman-test: Results of R coincide with those reported by Hollander &
Wolfe, which I'd take as a point in favor of R. In any case, my guess is that
ties are handled differently (aver
Thanks to Justin, Baptiste, and Sebed for your answers.
The solutions work well. I have been putting them to good use today: the
code now works wonderfully and I learnt some useful tricks!
thanks, Peter
<-Original Message->
>From: justin bem [justin_...@yahoo.fr]
>Sent: 9/8/2009 9:06:2
Here is a solutions using ggplot2 and reshape
library(reshape)
library(ggplot2)
data <- data.frame(k = 0:3, fk = c(11, 20,7,2), f0k = c(13.72, 17.64, 7.56,
1.08), fkest = c(11.85, 17.78, 8.89, 1.48))
Molten <- melt(data, id.vars = "k")
ggplot(Molten, aes(x = k, y = value, colour = variable)) + g
dear Dan,
As far as I know, the strucchange package can be helpful for you..
On the other hand, if your regression function is continuous at the
unknown break points to be estimated, you could try the segmented package.
Hope this helps you,
vito
Daniel Brewer ha scritto:
Hello,
I would lik
Dear useRs,
I want to plot the following barplot with lines instead of bars. Is there a way?
data <- data.frame(cbind(k = 0:3, fk = c(11, 20,7,2), f0k = c(13.72, 17.64,
7.56, 1.08), fkest = c(11.85, 17.78, 8.89, 1.48)))
d <- t(data[,2:4])
barplot(d, beside=TRUE)
Regards,
Rafael.
__
Not sure what you mean by 'store', but you can use a list:
> a <- matrix(1, nrow = 5, ncol = 1)
> b <- matrix(2, nrow = 10, ncol = 1)
> c <- matrix(3, nrow = 15, ncol = 1)
> myList <- list(a, b, c)
> str(myList)
List of 3
$ : num [1:5, 1] 1 1 1 1 1
$ : num [1:10, 1] 2 2 2 2 2 2 2 2 2 2
$ : num
> "Robin" == Robin Hankin
> on Tue, 08 Sep 2009 14:58:49 +0100 writes:
Robin> Hi guys
Robin> thanks for this, it works fine, but I'm not sure the Matrix package
does
Robin> what I want:
>> a = sparseMatrix(i=c(20, 30, 10), j=rep(1, 3), x=c(2.2, 3.3,
4.4))
Pharmacokinetic and pharmacodynamic modeling and simulation
By Dr. Jan Freijer
September 24, 2009
Amsterdam, The Netherlands
http://www.can.nl/events/details.php?id=57
This course is aimed at users of R or S-PLUS in the bio-pharmaceutical
sciences who would like to use R for clini
Gavin Simpson ucl.ac.uk> writes:
>
> On Fri, 2009-09-04 at 17:15 +0200, Kim Vanselow wrote:
> > Dear all,
> > I would like to calculate a cca (package vegan) with species and
> > environmental data. One of these environmental variables is
> > cos(EXPOSURE).
> > The problem: for flat releves ther
Hello,
I would like to test some data to see whether it has the shape of a step
function (i.e. y1 up until x_th and then y2 where x_th is the
threshold). The threshold x_th is unknown and the x values can only
take discrete values (0,1,2,3,4).
An example would be:
data<- data.frame(x=1:20,y=c(rn
Hi,
On Sep 8, 2009, at 9:56 AM, Abbas R. Ali wrote:
my dimentions of trining set dim(trainingset) = 7 x 96 and dim
(validation) = 3 x 96
other thing if i want to predicit trainingset accuracy it is also
giviing me same error. Now its no issue of dimentions from my side.
1. Please keep rep
Hi guys
thanks for this, it works fine, but I'm not sure the Matrix package does
what I want:
> a = sparseMatrix(i=c(20, 30, 10), j=rep(1, 3), x=c(2.2, 3.3,
4.4))
Error in asMethod(object) :
Cholmod error 'out of memory' at file:../Core/cholmod_memory.c, line 148
Surely an efficien
library(Matrix)
a = sparseMatrix(i=c(20, 30, 1), j=rep(1, 3), x=c(2.2, 3.3,
4.4))
b = sparseMatrix(i=c(3, 30), j=rep(1, 2), x=c(0.1, 0.1), dims=dim(a))
theSum = a+b
summary(theSum)
hth,
b
On Sep 8, 2009, at 10:19 AM, Henrique Dallazuanna wrote:
Try this:
abMerge <- merge(a, b, by
Hi,
On Sep 8, 2009, at 9:28 AM, Abbas R. Ali wrote:
by using predict(model.ksvm, validationset, ...) I am facing
following error:
"Error: '...' used in an incorrect context"
No. I didn't mean to explicitly use "...", sorry.
Actually, I had a two line fix here, but I'm looking back at your
Thanks a lot, Mohamed.
Kind regards,
Ezhil
--- On Mon, 9/7/09, Mohamed Lajnef wrote:
> From: Mohamed Lajnef
> Subject: Re: [R] Averaging rows if a condition is true.
> To: "A Ezhil"
> Cc: r-help@r-project.org
> Date: Monday, September 7, 2009, 9:22 PM
> Hi,
>
> Try to use aggregate function
On 08-Sep-09 13:06:28, Robin Hankin wrote:
> Hi
> I deal with long vectors almost all of whose elements are zero.
> Typically, the length will be ~5e7 with ~100 nonzero elements.
>
> I want to deal with these objects using a sort of sparse
> vector.
>
> The problem is that I want to be able to 'a
Dear All,
Has anyone tried to use this plug-in? Since I am running R-2.9.1 it will not
even let me install it. Further, since I am running Windows I cannot use the R
provided R-2.7.0 Linux installation file from the archive (tried to install it
through cygwin and it was a mess). Suggestions? Id
On Tue, Sep 8, 2009 at 1:52 PM, Lauri Nikkinen wrote:
> But this is not the solution I was looking for. Thanks.
I think the only way you'll get the solution you are looking for is
if you can let us have a copy of the original input file, or at least
the first few lines - and not pasted into an e
Hi,
On Sep 8, 2009, at 9:06 AM, Robin Hankin wrote:
Hi
I deal with long vectors almost all of whose elements are zero.
Typically, the length will be ~5e7 with ~100 nonzero elements.
I want to deal with these objects using a sort of sparse
vector.
Would using sparse matrices (from the Matrix
by using predict(model.ksvm, validationset, ...) I am facing following error:
"Error: '...' used in an incorrect context"
And I want to use regression in SVM what can I do for it.
Thanks
Abbas
--- On Tue, 9/8/09, Abbas R. Ali wrote:
From: Abbas R. Ali
Subject: [R] Regarding SVM using
Try this:
abMerge <- merge(a, b, by = 'index', all = TRUE)
list(index = abMerge$index, val = rowSums(abMerge[,2:3], na.rm = TRUE))
On Tue, Sep 8, 2009 at 10:06 AM, Robin Hankin wrote:
> Hi
>
> I deal with long vectors almost all of whose elements are zero.
> Typically, the length will be ~5e7 w
I have to start by saying that I am new to R, so I might miss something crucial
here. It seems to me that the results of friedman.test and ks.test are "wrong".
Now, obviously, the first thing which crossed my mind was "it can't be, this is
a package used by so many, someone should have observed"
1 - 100 of 152 matches
Mail list logo