OK, for the last step I have tried this (among other things):
library(data.table)
DT = data.table(read.table(textConnection("A B C
1 1 a 1999
2 1 b 1999
3 1 c 1999
4 1 d 1999
5 2 c 2001
6 2 d 2001
7 3 a 2004
8 3 b 2004
9 3 d 2004"),head=TRUE,stringsAsFactors=FALSE))
Dear Ramya
On 18 February 2011 05:08, Ramya wrote:
>
> http://r.789695.n4.nabble.com/file/n3312061/x_and_y_values.txt
> x_and_y_values.txt
>
>
> I have the absorbance values form HPLC Chromatogram. I need to find the
> peaks in datapoints and area under each peak.
>
> I am not sure how how to fin
On 02/18/2011 04:31 AM, Alaios wrote:
Dear all
when I call color.scale like this:
require('plotrix')
colcolor<-color.scale(c(range_sr,sr),extremes=c("red","blue"))
Error in if (min(reds)< 0 || max(reds)> 1) reds<- rescale(reds, c(0, :
missing value where TRUE/FALSE needed
range_sr
[1]
It has been a while back, but I believe I had to install libgtk2.0-dev
(that was on Ubuntu)
You could also try to install the r-cran-rgtk2 debian-package using
dpkg, aptitude, or whatever you use as package manager. This makes
rgtk available for all users.
HTH,
Jan
Quoting Lorenzo I
Hi John,
seems there is no easy way. I'll just precondition it with AWK as described
here http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg53401.html
There are some remarks in the thread that R is not supposed to read too large
files for "political" reasons. Maybe that's it.
Many thank
It has been a while back, but I believe I had to install libgtk2.0-dev
(that was on Ubuntu)
You could also try to install the r-cran-rgtk2 debian-package using
dpkg, aptitude, or whatever you use as package manager. This makes
rgtk available for all users.
HTH,
Jan
Quoting Lorenzo Isel
say, 'Dat' is your data frame, then one way to do it is:
with(Dat, tapply(freq, list(hair, eye, sex), c))
I hope it helps.
Best,
Dimitris
On 2/18/2011 8:50 AM, Lao Meng wrote:
The data is in the attachment.
What I wanna get is:
, , Sex = Male
Eye
HairBrown Blue Hazel Green
Thanks Jan, I am CCing to the list just to let everybody know that your
suggestion pretty much fixes my problem. Installing RGtk2 was the last
step to get rattle up and running. Somehow I had problems in running
rattle with the RGtk2 package from the standard debian testing
repositories (I supp
Damian Betebenner nciea.org> writes:
>
> In building a package, is it possible to make the package depend upon another
package that is only available
> on R-Forge (not CRAN). For example, by doing something in the DESCRIPTION file
> I'd like to add a dependency to my package such that when the u
On 02/18/2011 07:56 AM, juliemcwh wrote:
Sorry, I have read other posts on tick marks but nothing is working.
Here is my R code. I was able to turn off the x axis at one point (I
removed that code), but never got the tick marks on by 20s. It always goes
to the default of 50, 100, 150, 200.
T
Dear Jim,
this is just to confirm that it your solution worked :)
Regards
Alex
--- On Fri, 2/18/11, Jim Lemon wrote:
> From: Jim Lemon
> Subject: Re: [R] color.scale error
> To: "Alaios"
> Cc: R-help@r-project.org
> Date: Friday, February 18, 2011, 9:10 AM
> On 02/18/2011 04:31 AM, Alaios
> w
Hi:
Here's another approach using the reshape package:
library(reshape)
# id defines the grouping variables; the others are stacked, creating
# new variables 'variable' (containing the variable names as factor levels)
and
# 'value' (which holds the corresponding values)
dfm <- melt(df, id = c('co
On Fri, Feb 18, 2011 at 8:35 PM, Dennis Murphy wrote:
> Hi:
>
> Here's another approach using the reshape package:
>
> library(reshape)
> # id defines the grouping variables; the others are stacked, creating
> # new variables 'variable' (containing the variable names as factor levels)
> and
> # '
This is a built-in dataset in R - see ?HairEyeColor and str() it. I smell
homework...
Dennis
On Thu, Feb 17, 2011 at 11:50 PM, Lao Meng wrote:
> The data is in the attachment.
>
> What I wanna get is:
> , , Sex = Male
> Eye
> HairBrown Blue Hazel Green
> Black32 1110 3
Hi,
I am not a statistics expert, so I have this question. A linear model
gives me the following summary:
Call:
lm(formula = N ~ N_alt)
Residuals:
Min 1Q Median 3Q Max
-110.30 -35.80 -22.77 38.07 122.76
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(
Hello all,
In package "ade4", does anyone know whether there is a simple way of
setting the line type for the ellipses in an s.class plot? The
argument "lty" does not work.
Thanks,
Duarte
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailma
On Fri, 18 Feb 2011, Jan wrote:
Hi,
I am not a statistics expert, so I have this question. A linear model
gives me the following summary:
Call:
lm(formula = N ~ N_alt)
Residuals:
Min 1Q Median 3Q Max
-110.30 -35.80 -22.77 38.07 122.76
Coefficients:
Estima
Try this:
xtabs(Freq ~ Hair + Eye + Sex, Dat)
Using Dimitri's Dat example.
On Fri, Feb 18, 2011 at 5:50 AM, Lao Meng wrote:
> The data is in the attachment.
>
> What I wanna get is:
> , , Sex = Male
> Eye
> HairBrown Blue Hazel Green
> Black32 1110 3
> Brown53
--begin included message
Hi, when ploting Kaplan-Meier estimate curves as below, the censoring
symbols
(crosses) to not change thickness along the lines
plot(survfit(surv ~ I(x>=cut.off) ),lty=c(1,2), lwd=2)
is there any strightforward way to make it happen? thanks
-- End inclusion ---
Th
Hi:
On Fri, Feb 18, 2011 at 2:49 AM, Jan wrote:
> Hi,
>
> I am not a statistics expert, so I have this question. A linear model
> gives me the following summary:
>
> Call:
> lm(formula = N ~ N_alt)
>
> Residuals:
>Min 1Q Median 3Q Max
> -110.30 -35.80 -22.77 38.07 122.76
Hello
I have a multiple linear regression with two cofactors, I would like to
represent a plane but I could not find any help which worked out.
Any suggestions.
Regards and thanks in advance.
Rosario
__
R-help@r-project.org mailing list
https://stat.e
Hello,
is there an implementation available for a shared nearest neighbor
(SNN) clustering algorithm?
//Jay
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/pos
Dear R-people,
I am trying to read a matlab file, unfortunately without success. I enclosed
the zipped file.
Do you have an idea how to read this file by ReadMat package into the R system?
Thanks for your time!
Wim Krijnen
__
R-help@r-project.org
On Thu, 17 Feb 2011, Andrew Ziem wrote:
After ctree builds a tree, how would I determine the direction missing values
follow by examining the BinaryTree-class object? For instance in the example
below Bare.nuclei has 16 missing values and is used for the first split, but
the missing values
An update for the benefit of the list/posterity: I resolved this issue by
switching over to using the lrm() function of package rms. It seems to pick
better starts, or something; in any case, it has been able to converge on a
solution for every model I've tried, although for the most complex
Hello Everyone,
I'm pretty new to R and I'm trying to get some idea of the capabilities of
the language. I work with some pretty large data sets and the consensus
seems to be that R does not work well with big data. I've started talking to
the guys at Revolution, but I need to get some outside opi
Date: Fri, 18 Feb 2011 12:25:36 +0100
From: achim.zeil...@uibk.ac.at
To: jrheinlaen...@gmx.de
CC: r-help@r-project.org
Subject: Re: [R] lm without intercept
On Fri, 18 Feb 2011, Jan wrote:
> Hi,
>
> I am not a statistics expert, so I have this ques
Hi,
I am modeling a gen linear mixed model with binomial repeated measures
responses (y=cbind(correct/not correct)) of some individuals (id).
I have a 2X2 design (effort and costs) and another continuous independent
variable (scr). This will lead to something like this:
mod.1<-lmer(y~effort+cost
I'm attempting to sort a 3 dimensional array that looks like this
> x
, , 1
[,1] [,2]
[1,]99
[2,]79
, , 2
[,1] [,2]
[1,]65
[2,]46
, , 3
[,1] [,2]
[1,]21
[2,]32
Such that it ends up like this
> y
, , 1
[,1] [,2]
[1,]21
No, this is a cute problem, though: the definition of R^2 changes
without the intercept, because the
"empty" model used for calculating the total sums of squares is always
predicting 0 (so the total sums
of squares are sums of squares of the observations themselves, without
centering around the sam
Krijnen WP, Wim pl.hanze.nl> writes:
>
> Dear R-people,
>
> I am trying to read a matlab file, unfortunately without success.
I enclosed the zipped file.
>
> Do you have an idea how to read this file by ReadMat package
into the R system?
>
> Thanks for your time!
>
> Wim Krijnen
>
>
I
sprocket gmail.com> writes:
>
>
> Hi,
>
> I am modeling a gen linear mixed model with binomial
repeated measures
> responses (y=cbind(correct/not correct)) of some
individuals (id).
> I have a 2X2 design (effort and costs) and another
continuous independent
> variable (scr). This will lead
Dear Jim,
This works for your small example. I was not exactly sure what you
were sorting on. In my example, I sum all values at each level of the
third dimension and sort by that sum. Of course, if there are many
levels or the 2 dimensions being summed are quite large, this will be
quite compu
Hello Achim,
> Not quite. Consult your statistics textbook for the correct interpretation
> of p-values. Under the null hypothesis of a true intercept of zero, it is
> very likely to observe an intercept as large as 13.52 or larger.
thank you for that help. I suppose the net doesn't have a detai
Hi,
thanks for your help. I'm beginning to understand things better.
> If you plotted your data, you would realize that whether you fit the
> 'best' least squares model or one with a zero intercept, the fit is
> not going to be very good
> Do the data cluster tightly around the dashed line?
No, a
On 18/02/2011 5:44 AM, Michael Holt wrote:
Hello Everyone,
I'm pretty new to R and I'm trying to get some idea of the capabilities of
the language. I work with some pretty large data sets and the consensus
seems to be that R does not work well with big data. I've started talking to
the guys at R
The following is my R script which I am struggling with to assess ICESat
data..perhaps it is the ID_min or ID_max that is wrong? I don't know, any
help would be greatly appreciated :(
# OPTIONS - CHANGE THESE VARIABLES IF NEEDED\par
###
> Date: Fri, 18 Feb 2011 08:39:05 -0500
> From: murdoch.dun...@gmail.com
> To: mich...@aers.ca
> CC: r-help@r-project.org
> Subject: Re: [R] Boundaries of R
>
> On 18/02/2011 5:44 AM, Michael Holt wrote:
> > Hello Everyone,
> >
> > I'm pretty new
It looks to me like when you assign ID_min and ID_max that you want them
to be vectors. In R, this is accomplished by ?c.
Example:
ID_min <- c(976834501, 290018099.3145, 1000)
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Ke
danielepippo wrote:
>
>I'm building a matrix in R with a cycle for like this:
> pp_ris2=matrix(NA,6,6)
> for(i in 0:6){
> .
>
R is not like c, indexing starts with 1.
Dieter
--
View this message in context:
http://r.789695.n4.nabble.com/Matrix-in-R-tp3312748p3312764.html
Sent from
Dear James,
this is what I understood your sorting along the third dimension to be:
> x <- array(c(9, 9, 7, 9, 6, 5, 4, 6, 2, 1, 3, 2), dim = list(2, 2, 3))
> y <- apply (x, 1:2, sort)
> y
, , 1
[,1] [,2]
[1,]21
[2,]65
[3,]99
, , 2
[,1] [,2]
[1,]32
[2
Mike Marchywka hotmail.com> writes:
> Duncan Murdoch:
> > Vectors are limited to about 2 billion entries (2^31 - 1). Matrices are
> > vectors, so that limit applies to the total count of entries.
> > Dataframes are lists of vectors, so that limit applies separately to the
> > numbers of rows and
Hi Rex,
Thanks, this is exactly what I want but have to do it with many big arrays ...
thus if there were a way to do it with a vectorized function would it not be a
lot more efficient?
Much appreciated!
J
>Subject: RE: sort a 3 dimensional array across third dimension ?
>
>Although I suggest
but if in my function
pp_ris2[i,j]=myfunction}
must be the indexes 0-0,0-1,0-2,0-3, ?
--
View this message in context:
http://r.789695.n4.nabble.com/Matrix-in-R-tp3312748p3312780.html
Sent from the R help mailing list archive at Nabble.com.
__
R
Hi everyone,
I'm building a matrix in R with a cycle for like this:
pp_ris2=matrix(NA,6,6)
for(i in 0:6){
for(j in 0:6){
if(i>j){
pp_ris2[i,j]=myfunction}
else if(i==j){
print(c(i,j))
pp_ris2[i,j]=myfunction}
Dear Uwe,
Thanks for your time.
Definitely something was bugging my script (my mistake).
Since your replay I have been trying to find the bug, to a simple script :( ,
without success.
Preparing a reply to post with data sample and the script, everything start to
work normally.
Sadl
Hi All,
I would like to know where to find technical papers to how calculate ANOVA
using ERGM. Which parameters and so on. Thanks.
--
View this message in context:
http://r.789695.n4.nabble.com/ERGM-Anova-tp3312792p3312792.html
Sent from the R help mailing list archive at Nabble.com.
_
hi everybody,
i know it a quite complicate subject but someone might have the
solution.
I am doing a delta model coupling a binomial glm and a lognormal one.
Using the Laurent correction I can predict mean values and I would
like to know if you know how to predict the variance?
Do you know the t
Hello
I've created a 3d scatterplot, and had no problems labeling the points.
However, I've been really struggling to change the color of the points based
on a factor (see 'group' below). Is such a thing possible?
My data look like this:
food group x y z
apple fruit 0.216 -0.11
Has anyone any experience of applying observational weights in bagging? I am
performing regression trees (continuous data on bird abundance) and need to
account for sampling intensity. In a single tree, i.e. a call of rpart, I can
specify weights either by having a separate vector called weights
Yes, remove the call to intersect, and rely on the results of match to tell you
whether there is an overlap. If there are any matches, all(is.na(index)) will
be false. Read help for match.
?match
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.
I apologize if you have already seen my question. I am new to the mailing list
and I did not get any responses the first time I posted my question. However,
I am not sure my post went through.
Here is my situation:
I have a spreadsheet with columns of fish species (text) and length (numb
Although I suggested to someone else that for-loops be avoided, they are not in
the inner loop in this code, and it's probably easier to understand than some
sort of apply:
a = array(round(100*runif(60)),dim=c(3,4,5))
a
for (i in 1:dim(a)[1])
for (j in 1:dim(a)[2])
a[i,j,] = sort(a[i,j,])
a
On 02/18/2011 04:11 PM, Maas James Dr (MED) wrote:
Hi Claudia,
It does help a lot, but not quite there yet ... I'm sure you are correct and
is much appreciated, I need some sort of generalized form, actual arrays in
my case are 3x3x1000. Do you suspect it could be done in one step with
sapply?
OK thanks Jonathon, I will have a go.
Is the plotdir file technically my output file? Or what is this otherwise?
Also, having never used R script before, is it not very good at reading
data, may I have to alter my input file in excel first?
Thanks
--
View this message in context:
http://r.789
Have you read an Intro to R or, indeed, any R documentation? Perhaps
you received no reply because you do not appear to have made an honest
effort to answer the question for yourself. Maybe you should consider
accessing R thru R Commander or one of the other R GUI's -- search
CRAN for details.
--
Dear Users,
>From previous analysis (semi-variograms using package gstat), I found
spatial autocorrelation in my dataset.
The best fitted model to this spatial correlation structure is the Gaussian
model (Spherical, Exponential, Linear tested and comparison done by Sum of
Square errors).
So I used
On Fri, Feb 18, 2011 at 06:32:01AM -0800, danielepippo wrote:
>
> but if in my function
> pp_ris2[i,j]=myfunction}
> must be the indexes 0-0,0-1,0-2,0-3, ?
You'll have to take care of that yourself with a bit of index
arithmetics. It's the same you encounter in C, if you are
modelling somet
What does this have to do with R?
-- Bert
On Fri, Feb 18, 2011 at 6:37 AM, Roger Gomes wrote:
>
> Hi All,
>
> I would like to know where to find technical papers to how calculate ANOVA
> using ERGM. Which parameters and so on. Thanks.
> --
> View this message in context:
> http://r.789695.n4.nab
On Feb 18, 2011, at 9:32 AM, danielepippo wrote:
but if in my function
pp_ris2[i,j]=myfunction}
must be the indexes 0-0,0-1,0-2,0-3, ?
From a search of RSiteSearch() that started with with terms: zero
matrix indexing
http://finzi.psych.upenn.edu/R/Rhelp02/archive/39031.html
http://
On Feb 18, 2011, at 9:32 AM, danielepippo wrote:
but if in my function
pp_ris2[i,j]=myfunction}
must be the indexes 0-0,0-1,0-2,0-3, ?
I came across a posting in r-help that called this package "blasphemy":
http://cran.r-project.org/web/packages/Oarray/index.html
--
David Winsemius
Dears UseRs,
In Zuur's book (Mixed effects models and extensions in ecology with R), p
167, it is mentionned that the Variogram function (package nlme) scale the
sill to 1.
Here is the plot of the semi-variogram given by the Variogram function on my
"lme" object. (http://imm.io/3OLe)
Note that the
Well, this script only shows the option settings, so I assumed that the
string being used as your infile was formatted for R. As for how it will
have to look for R to read it in, there are many ways to do this - I would
start by reading "R Data Import/Export" that ships with every copy of R.
As
I echo Prof. Ripley's and Jari Oksanen's concerns. I see two
options in case the R-Forge package maintainer does not submit the
package to CRAN:
(1) Copy the function(s) you want from that package into
your own. In such cases, I recommend you change the name of the
functi
Hi,
I am running a simulation study with the survival::coxph. Some of the
simulations result in problematic fits due to flat partial likelihood. So, you
get the warning message:
Warning message:
In fitter(X, Y, strats, offset, init, control, weights = weights, ... :
Loglik converged before
May be I could do:
for (i in 1:nsim) {
last.warning <- NULL
# do model fitting
if(!is.null(last.warning)) # discard simulation result
}
I think this might work. Any other ideas?
Ravi.
Ravi Varadhan, Ph.D.
Assistant Profe
Hi all,
I am new to heatmap function in stats package, so any advice will be
greatly appreciated. So I have a dataset that I want to draw a heatmap,
for example:
> dput(a)
structure(c(-8.256023, -14.764429, -10.145779, -5.123104, -11.974729,
-16.329029, -16.329029, -16.329029, -2.764114, -4.9
Nope - that does not work. The value of last.warning is not reset after the
initial NULL.
Ravi.
Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins Univers
On Feb 18, 2011, at 11:35 AM, Ravi Varadhan wrote:
May be I could do:
for (i in 1:nsim) {
last.warning <- NULL
# do model fitting
if(!is.null(last.warning)) # discard simulation result
You can raise warnings to behave as errors with options:
options(warn=2)
And you can define where R go
If you have a function that may produce warnings you
can use withCallingHandlers() to make a new function
that does what the original does but adds some
contextual information to the usual warning report.
E.g.,
addWarningReporting <- function(FUN) {
force(FUN)
function(args) {
withCall
On Feb 18, 2011, at 9:45 AM, Mike Gibson wrote:
I apologize if you have already seen my question. I am new to the
mailing list and I did not get any responses the first time I posted
my question. However, I am not sure my post went through.
It didn't.
Here is my situation:
I have a s
hello list,
I have this array of distribution of frequencies:
freq<-c(45,18,10,8,13,5,9,3,4,1,2,2,2,3,4,4,1,1,2,1,2,1,1,1,1,1,2,3,1,1,1,1,1,1,1,1,2,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
for this array of values:
values<-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,16,18,19,20,21,
Thanks Claudia (and all others!0
This works beautifully! The instructions for aperm are a bit vauge!
J
--
View this message in context:
http://r.789695.n4.nabble.com/sort-a-3-dimensional-array-across-third-dimension-tp3312612p3313126.html
Sent from the R help mailing list archive at Nabble.co
Dear all
I'm trying to install the hdf5 library into R. The HDF5 package is
installed in a special directory, distributed accross my cluster:
/share/apps/HDF5
So i put the args option to the configure command as i read in previous
post in the list:
jerome]# R CMD INSTALL --configure-vars='
I was going to say:
The problem with for-loops (as best I understand it) is that the R code gets
interpreted over and over; what you normally want to do is design the
computation so that you jump into the internals of R and stay there. But the
inner loop is in the R internals of the sort in th
Hi Bert, thanks about your answer.
Well, I would like to calculate anova.ergm in R. Wich values (edges,
vertices) the library use to calculate F values or p-value with graphos.
The calculus is like a Pandora Box... I put there a network (ergm) and then,
values show me like a magic.
Do you know
i have got a svm MODEL from e1071
by
svmModel<-svm(x, y, kernel = "linear", cost = gamma, scaled=FALSE)
who can tell me how can i get the alpha of decision functions
i just remembered rho is the bias of decision functions
i think the manual of e1071 is a piece of usely garbage
thank u
gao shan
Hello
I would like to know if the Diebold Mariano Test in the forecast Package is
adjusted to small samples (as Harvey, Leybourne, Newbold suggest) If not, how
can I do that manually?
Paka
__
R-help@r-project.org mailing list
https://stat.ethz.ch
Kari-4 wrote:
>
>
> I've created a 3d scatterplot, and had no problems labeling the points.
> However, I've been really struggling to change the color of the points
> based
> on a factor (see 'group' below). Is such a thing possible?
>
>
Its always polite to make the example self-running by
Dear R-users,
Is there any way of calculating the perimeter of a polygon in a shapefile
object? Furthermore, how to calculate the length of the common border of two
polygons?
I've searched the code of spded, but could not find a hint on how to do it.
Thank you very much,
Leo.
[[alternativ
Hi All,
I am using the package, cmprisk, to plot competing risks. In my case, I have
four lines showing risk of going on dialysis (by a lab test [fgf-23] in
quartiles), where the 4 lower lines are for the competing risk of death.
I am trying to edit the function to just plot the upper curves (
If do:
> library("e1071")
> example(svm)
I get:
svm> data(iris)
svm> attach(iris)
svm> ## classification mode
svm> # default with factor response:
svm> model <- svm(Species ~ ., data = iris)
svm> # alternatively the traditional interface:
svm> x <- subset(iris, select = -Species)
svm> y <- S
Hi everyone,
if I have the bootsrap values of the parameters of my regression model,
how can I calculate the standard error of them to refer to the parameters?
Thank you
--
View this message in context:
http://r.789695.n4.nabble.com/Bootstraps-standard-error-tp3313322p3313322.html
Sent from
I'm trying to use lattice graphics to produce some small plots for inclusion in
a LaTeX file. I want the LaTeX fonts to be used in the plots, but to be scaled
down to match the size of the plot. I have written the following code to apply
a scaling factor to all the "cex" and "padding" entries in
All,
I'd like to automatically output text from R to HTML. In doing this I've
run into trouble with non-ascii characters, as my browser (and
presumably others) does not render such characters correctly. For
example, the 'fancy' single quotes associated with summary.lm are
multi-byte character
Hi,
Can anyone help me to plot a chart graph? I have a data set like this. I
would like a bar chart graph which mouse1's two treatments are clustered
together, so on.
I tried with barplot, but couldn't get it right.
Thanks
Johnny
Mouse1
Mouse2
Mouse3
Tal,
One interactive capability that I have repeatedly wished for (but
never taken the time to develop with the existing R tools) is the
ability to interactively zoom in on and out of a data set, and to
interactively create "call-outs of sections of the data. Much of the
data that I deal with take
On Feb 18, 2011, at 1:36 PM, danielepippo wrote:
Hi everyone,
if I have the bootsrap values of the parameters of my regression
model,
how can I calculate the standard error of them to refer to the
parameters?
The same way you would any other vector?
Thank you
--
View this message i
On 2011-02-18 11:16, Juergen Rose wrote:
If do:
library("e1071")
example(svm)
I get:
svm> data(iris)
svm> attach(iris)
svm> ## classification mode
svm> # default with factor response:
svm> model<- svm(Species ~ ., data = iris)
svm> # alternatively the traditional interface:
svm> x<
On Feb 18, 2011, at 2:10 PM, Gerard Smits wrote:
Hi All,
I am using the package, cmprisk, to plot competing risks. In my
case, I have four lines showing risk of going on dialysis (by a lab
test [fgf-23] in quartiles), where the 4 lower lines are for the
competing risk of death.
I am t
Am Freitag, den 18.02.2011, 11:53 -0800 schrieb Peter Ehlers:
> On 2011-02-18 11:16, Juergen Rose wrote:
> > If do:
> >> library("e1071")
> >> example(svm)
> >
> > I get:
> >
> >
> > svm> data(iris)
> >
> > svm> attach(iris)
> >
> > svm> ## classification mode
> > svm> # default with factor res
Dear Jürgen,
did you update.packages (checkBuilt = TRUE) ?
I recently had segfaults, too on 64bit linux (with rgl, though) and they
disappeared only after updating with checkBuilt (including also the packages
originally installed via Dirk's .deb packages.
HTH,
Claudia
On 02/18/2011 09:32
On 2011-02-18 12:32, Juergen Rose wrote:
Am Freitag, den 18.02.2011, 11:53 -0800 schrieb Peter Ehlers:
On 2011-02-18 11:16, Juergen Rose wrote:
If do:
library("e1071")
example(svm)
I get:
svm> data(iris)
svm> attach(iris)
svm> ## classification mode
svm> # default with factor res
Am Freitag, den 18.02.2011, 21:42 +0100 schrieb Claudia Beleites:
> update.packages (checkBuilt = TRUE)
Dear Claudia,
I just tried update.packages (checkBuilt = TRUE), it did not solved the
issue.
Nevertheless thanks,
Juergen
__
R-help@r-project.org
Listers,
Is there a command/function to get the population standard deviation (N) and
the sample standard deviation (n-1)
Thanks
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-
On 18-Feb-11 20:51:40, Dmitry Berman wrote:
> Listers,
>
> Is there a command/function to get the population standard deviation
> (N) and
> the sample standard deviation (n-1)
>
> Thanks
Using data (1:10), with N=10, as an example:
c(SampSD=1, PopSD=sqrt(9/10))*sd((1:10))
# SampSDPopS
Dear Jürgen,
I just tried update.packages (checkBuilt = TRUE), it did not solved the
issue.
Was just a guess.
Nevertheless thanks,
Bitte, bitte.
In case that helps tracking down the problem:
it works with my normal R as well as with r-devel. In case it helps you, I also
put the installati
Thank you Joshua and Jim--I got the tick marks on finally!
One last question: I am supposed to make bins in 1, 2 and 5 mm increments for
fish up to 170 mm. I assume these are the breaks? So breaks = 170 for 1 mm
and 340 for 2 mm and 1550 for 5 mm? I couldn't find any specific examples in
the
Sorry the data were all messed up. Let me try again.
Can anyone help me to plot a chart graph? I have a data set like this. I
would like a bar chart graph which mouse1's two treatments are clustered
together, so on.
I tried with barplot, but couldn't get it right.
Mouse1 Mouse2 Mouse3
Hi all,
I'm currently trying to get accustomed to R and some of its packages, thus
need some help for a package called R2wd. As I am trying to use the wdGET()
function I'm getting the following message:
Error in if (!(tmp[["ActiveDocument"]][["Name"]] == filename))
tmp$Open(paste("path", :
1 - 100 of 124 matches
Mail list logo
