Thank you Steve for the solution: As per your suggestion I spend some time
to make it work for 2 variables.
nvar = 3 # number of variables
ncol<- nvar*2
func1<- function(x) {sapply( strsplit(as.character(x), ""),
match, table= c("1","2","3","4",NA))}
mydf1 <- data.frame(t
Hi Max,
I tried your suggestion but cam up with errors:
fitControl<-trainControl(number=1)
LOGISTIC_model <- train(mdrrDescr,mdrrClass
,method='glm'
,trControl = fitControl
)
Fitting: parameter=none
Error in if (all.equal(so
That's great - thanks all for your help.--
View this message in context:
http://r.789695.n4.nabble.com/Sorting-dataframe-by-number-of-occurrences-of-factor-tp3485443p3488978.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-proj
Thanks a lot! It works!
xiagao1982
2011-05-02
åä»¶äººï¼ Henrique Dallazuanna
åéæ¶é´ï¼ 2011-05-02 08:31:02
æ¶ä»¶äººï¼ xiagao1982
æéï¼ r-help
主é¢ï¼ Re: [R] How to pass objects from local() to GlobalEnv
Try this:
local(x <<- 1)
On Sun, May 1, 2011 at 9:20 PM, xiagao1982 w
Try this:
local(x <<- 1)
On Sun, May 1, 2011 at 9:20 PM, xiagao1982 wrote:
> Hi all,
>
> I create some objects in local(), and want to pass them to GlobalEnv. How can
> I do this? Thanks!
>
>
>
>
> xiagao1982
> 2011-05-02
>
> [[alternative HTML version deleted]]
>
> _
Hi all,
I create some objects in local(), and want to pass them to GlobalEnv. How can I
do this? Thanks!
xiagao1982
2011-05-02
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/
Unfortunately, it seems that vcd doesn't return grobs but draws
directly to the device, which prevents a concise solution. You could
try the following,
library(gridExtra)
library(vcd)
data("Titanic")
p = grid.grabExpr(mosaic(Titanic))
grid.arrange(p, p, p, ncol=2)
Or, more versatile but also mor
Thanks again Max - a great time saver this is.
Now just for my sanity, if I use glm.fit to build a model where I have the
matrices, how do I then use the predict function without getting an error
message?
> LOGISTIC_model1 <- glm.fit(mdrrDescr,mdrrClass,
> family=binomial(link="logit"))
Warning
Here is one approach. It would be good to provide a reasonable sample of data:
> x <- unclass(Sys.time()) # today's date
> # create some data
> # increments by ~ 0.1 seconds
> len <- cumsum(runif(100, 0, 0.1))
> dataFile <- data.frame(time = x + len,
+flag = sample(c("Y",
Not all modeling functions have both the formula and "matrix"
interface. For example, glm() and rpart() only have formula method,
enet() has only the matrix interface and ksvm() and others have both.
This was one reason I created the package (so we don't have to
remember all this).
train() lets yo
I would like to display multiple mosaic plots from vcd (not defined by a model
but derived from different data sets)
side by side.
Neither par(mfrow=...)
nor layout seem to allow to arrange multiple mosaic plots in a grid.
Is there an easy way of arranging mosaics in a grid?
_
Hi Yao,
I can't answer that question, but I offer the following thoughts for
your consideration.
Generally it's best to approach the package maintainers directly with
questions like these. You can find their contact details in the
package documentation.
Also, you will want to make sure that yo
glm.fit - answered my own question by reading the manual!--
View this message in context:
http://r.789695.n4.nabble.com/caret-prevent-resampling-when-no-parameters-to-find-tp3488761p3488923.html
Sent from the R help mailing list archive at Nabble.com.
_
This is a nice demonstration of the formula interface to aggregate. A
less elegant alternative is to pass lists as arguments.
with(dd,
aggregate(Correct,
by = list(Subject = Subject,
Group = Group),
FUN = function(x) sum(x == 'C')))
Thanks for the clarification Max - I should have realised that.
One final question, I like caret because it lets me pass in data to all
functions in the same way. For glm I have only ever used the formula
notation and did not see a way to pass in predictors and a target
individually. How do I do t
No, the sampling is done on rows. The definition of a bootstrap
(re)sample is one which is the same size as the original data but
taken with replacement. The "Accuracy SD" and "Kappa SD" columns give
you a sense of how the model performance varied across these bootstrap
data sets (i.e. they are not
When you say "variable" do you mean predictors or responses?
In either case, they do. You can generally tell by reading the help
files and looking at the examples.
Max
On Fri, Apr 29, 2011 at 3:47 PM, ypriverol wrote:
> Hi:
> I'm starting a research of Support Vector Regression. I want to obta
Hi Max,
But in this example, it says the sample size is the same as the total number
of samples, so unless the sampling is done by columns, wouldn't you get
exactly the same model each time for logistic regression?
ps - great package btw. I'm just beginning to explore its potential now.--
View th
It isn't building the same model since each fit is created from
different data sets.
The resampling is sort of the point of the function, but if you really
want to avoid it, supply your own index in trainControl that has every
index (eg, index = seq(along = mdrrClass)). In this case, the
performan
One thing to watch out for using POSIXct is 1us is about the limit of
accuracy due to floating point (see FAQ 7.31). Notice that printing
out today's date requires about 15 digits with microsecond
granularity. Notice in the example that if the time difference
between intervals is 0.1 us, you have
I want to use caret to build a model with an algorithm that actually has no
parameters to find.
How do I stop it from repeatedly building the same model 25 times?
library(caret)
data(mdrr)
LOGISTIC_model <- train(mdrrDescr,mdrrClass
,method='glm'
Hi,
There are a couple of ways to do what you want.
I'll provide the fodder and let you finish the implementation.
On Sun, May 1, 2011 at 4:26 PM, Ram H. Sharma wrote:
> Hi R experts
>
> I have a couple of quick question:
>
> Q1
> #my data
> set.seed(12341)
> SN <- 1:100
> pool<- c(12,13,14, 2
On Apr 30, 2011, at 4:18 PM, Alice Wines wrote:
Hello all,
I have a quandry I have been scratching my head about for a
while. I've searched the manual and the web and have not been able to
find an acceptable result, so I am hoping for some help.
I have two data frames and I want to in
Can anyone point me to examples with R code where bwplot in lattice is used to
order the boxes in ascending order? I have found the following discussion and
it partly works. But, I have a conditioning variable, so my example is more like
bwplot(var1 ~ var2|condition, dat)
Th example in the disc
Apology for reposting, but the format of earlier message got
distorted; hopefully this time it will be readable:
From: wilds...@hotmail.com
To: r-help@r-project.org
Subject: Longitudinal data with non-randomized subjects
Date: Sun, 1 May 2011 00:34:08 -0700
Dear List,
I have a theoretical quest
Hi R experts
I have a couple of quick question:
Q1
#my data
set.seed(12341)
SN <- 1:100
pool<- c(12,13,14, 23, 24, 34)
CT1<- sample(pool, 100, replace= TRUE)
set.seed(1242)
CT2 <- sample(pool, 100, replace= TRUE)
set.seed(142)
CT3 <- sample(pool, 100, replace= TRUE)
# the number of variables r
Hi,
many thanks for helpful answer.
Best
Marcin M.--
View this message in context:
http://r.789695.n4.nabble.com/Kolmogorov-Smirnov-test-tp3479506p3488364.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing l
For this case B=200 should work well if using the bootstrap. For cross-val.
you can use B=10-fold cross-val and repeat the process 100 times for
adequate precision, averaging over the 100 as done in
http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/logistic.val.pdf (note
this was using the Design
Hi:
Instead of using ifelse(), you can multiply the logical statement by
the rest of the expression. When the logical statement is false, its
numerical value is zero.
HTH,
Dennis
On Sun, May 1, 2011 at 9:52 AM, Andrew Coop wrote:
> I have data vectors x and y both with 179 observations. I'm tr
Hi:
I would do something like the following:
(1) Create a vector of the file names.
(2) Use lapply() to read the files into a list.
(3) Use the reshape or reshape2 package to melt the individual files
into 'long' form.
(4) rbind together the resulting data frames.
(5) Use a summarization function
On Sun, May 1, 2011 at 15:33, Joel Reymont wrote:
>
> Does R have support for microseconds in timestamps, e.g. when reading this in
>
> "Time","Include","Kind","Duration"
> 2011-04-01 14:20:36.368324,Y,U,1.03238296509
> 2011-04-01 14:20:35.342732,Y,C,0.0252721309662
> 2011-04-01 14:20:34.337209,Y,
Hi,
thank you very much, both methods worked perfectly.
Regards
On Fri, Apr 29, 2011 at 4:17 PM, Berend Hasselman wrote:
>
> David Winsemius wrote:
> >
> > On Apr 29, 2011, at 4:27 AM, ivan wrote:
> >
> >> Hi All,
> >>
> >> I am trying to create a function which evaluates whether the values
>
I have data vectors x and y both with 179 observations. I'm trying to
fit a nonlinear model with five parameters using nls. The formula is
only defined within a range of x-values, it should be zero otherwise,
thus my attempted use of ifelse:
> df<-data.frame(x,y)
> nlsfit<-nls(y~ifelse(x>m&xhttp
Does R have the capability to perform tests for the need of clustering
analysis (e.g., in prabclus)? I am using panel data with two-way fixed
effects but am unsure about whether I should be using a cluster option as
well to estimate my model.--
View this message in context:
http://r.789695.n4.nab
I have 100+ .csv files which have the basic format:
> test
X Substance1 Substance2 Substance3 Substance4 Substance5
1 Time1 10 0 0 0 0
2 Time2 9 5 0 0 0
3 Time3 8 10 1
Does R have support for microseconds in timestamps, e.g. when reading this in
"Time","Include","Kind","Duration"
2011-04-01 14:20:36.368324,Y,U,1.03238296509
2011-04-01 14:20:35.342732,Y,C,0.0252721309662
2011-04-01 14:20:34.337209,Y,R,0.00522899627686
Thanks, Joel
--
I'm trying to run "rfe" for variable selection in the caret package, and am
getting an error. My data frame includes a dummy variable with 3 levels.
x <- chlDescr
y <- chl
#crate dummy variable
levels(x$State) <- c("AL","GA","FL")
dummy <- model.matrix(~State,x)
z <- cbind(dummy, x)
#remove Stat
On 29.04.2011 15:06, Pete Brecknock wrote:
... is the apply function what you are looking for?
A=matrix(1,2,4)
apply(A,1,sum)
Thanks for providing answers to R-help, but:
1. Please quote the original question for the mailing list readers of us.
2. Please reply also to the original poster
Hmmm, your subject says you are going to get statmod,
On 28.04.2011 04:31, Jim Silverton wrote:
Hello all,
I keep on getting the following error message when I try downloading
statmod:
install.packages("statmod")
Installing package(s) into ‘C:\Users\Isaac\Documents/R/win-library/2.12’
(as ‘li
On 29.04.2011 22:20, hck wrote:
Barth sent me a very good code and I modified it a bit. Have a look:
Error<-rnorm(1000, mean=0, sd=0.05)
estimate<-(log(1+0.10)+Error)
DCF_korrigiert<-(1/(exp(1/(exp(0.5*(-estimate)^2/(0.05^2))*sqrt(2*pi/(0.05^2
))*(1-pnorm(0,((-estimate)/(0.05^2)),sqrt(1/(
On 01.05.2011 05:39, Chee Chen wrote:
Dear All,
For function: rmvnorm{mvtnorm} in (library
No, it is a package, not a library!
mvtnorm, not splus2R), if I generate 2 bivariate normal samples as follows:
rmvnorm(2,mean=rep(0,2),sigma=diag(2))
[,1] [,2]
[1,] 2.0749459 1.
On 29.04.2011 14:39, Jim Lemon wrote:
On 04/29/2011 08:35 PM, hck wrote:
Dear all
Problem: hist()-function, scale = “percent”
I want to generate histograms for changing underlying data. In order
to make
them comparable, I want to fix the y-axis (vertical-axis) to, e.g.,
0%, 10%,
20%, 30% as
On Sat, Apr 30, 2011 at 10:26 AM, xiagao1982 wrote:
> Hi all,
>
> I am a C++/C# programmer who is new to R. I would like to use something like
> "namespace" to organize my functions without creating a package. How can I do
> this? Thanks!
>
You can arrange them in classes using reference classe
Please upgrade to a current release (or R-patched). Version 2.10.1 is
quite old.
Duncan Murdoch
On 30/04/11 3:18 PM, Oliver wrote:
Hello,
when using
Sys.getenv() during startup-phase (.First or .Rprofile)
to get the env-variables
COLUMNS as well as HOST I get empty strings.
After the sta
On 30/04/11 10:26 AM, xiagao1982 wrote:
Hi all,
I am a C++/C# programmer who is new to R. I would like to use something like
"namespace" to organize my functions without creating a package. How can I do
this? Thanks!
You could do it with the local() function or other explicit use of
environ
Thanks for the answer and for the link. I was lookin for a search trough the
forum posts
So the slope of the line is not important as long as the data is approx. on the
line?
Thanks, m
-Original Message-
From: Joshua Wiley [mailto:jwiley.ps...@gmail.com]
Sent: Saturday, April 30, 2
Dear List,
I have a theoretical question related to epidemiological data analysis:
If the treatment status (tx = 0,1) changes over time for the patients in a
non-randomized cohort, is there a way to estimate the treatment effect?
(i.e., after joining the study, some patients may have to wait fo
On Sat, Apr 30, 2011 at 11:39:43PM -0400, Chee Chen wrote:
> Dear All,
> For function: rmvnorm{mvtnorm} in (library mvtnorm, not splus2R), if I
> generate 2 bivariate normal samples as follows:
> > rmvnorm(2,mean=rep(0,2),sigma=diag(2))
>[,1] [,2]
> [1,] 2.0749459 1.4932752
> [
Dear R users:
Recently, I learn to use penalized logistic regression. Two packages
(penalized and glmnet) have the function of lasso.
So I write these code. However, I got different results of coef. Can someone
kindly explain.
# lasso using penalized
library(penalized)
pena.fit2<-penalized(HRLNM,
On Apr 30, 2011, at 6:06 PM, Thomas Lumley wrote:
On Sun, May 1, 2011 at 4:49 AM, David Winsemius > wrote:
On Apr 30, 2011, at 10:44 AM, Jabba wrote:
Dear useRs,
I was asked to produce a survival curve like this:
http://www.palug.net/Members/jabba/immaginetta.png/view
with the cardinality
On Sun, May 1, 2011 at 8:08 AM, azam jaafari wrote:
> Dear All
>
> I want to import the vector file ( .shp) to R. I could import the file by
> rgdal package before, by following:
>
> geology<-readOGR('C:/geology//saga/geo.geom','finalgeology')
>
> but now there is an error:
>
> Error in ogrInfo
Dear All
I want to import the vector file ( .shp) to R. I could import the file by
rgdal package before, by following:
geology<-readOGR('C:/geology//saga/geo.geom','finalgeology')
but now there is an error:
Error in ogrInfo(dsn = dsn, layer = layer, input_field_name_encoding =
input_fi
52 matches
Mail list logo
