On 5 January 2015 at 21:08, Ben Bolker wrote:
> Roger Coppock cox.net> writes:
>
>>
>> When will "R" implement the "se.fit" option to the
>> predict.nls() function? Is there some schedule?
>>
>
> I think this is unlikely to happen, ever (sorry). The exact method
> for finding confidence inte
Because data.frame converts it. I don't know what rationale was originally used
to justify including that transformation automatically, but I think it is much
better this way than having POSIXlt in the data frame. POSIXlt is itself like a
little data frame, and having data frames inside of data
On 2015/1/5 20:51, David Winsemius wrote:
On Jan 5, 2015, at 5:47 PM, Jinsong Zhao wrote:
On 2015/1/5 17:28, Ben Bolker wrote:
Jinsong Zhao yeah.net> writes:
In the following code snippet,
#
a <- strptime("121114 0510", "%m%d%y %H%M")
b <- data.frame(date = a, res = 1:5)
class(a)
class(b[
On Jan 5, 2015, at 5:47 PM, Jinsong Zhao wrote:
> On 2015/1/5 17:28, Ben Bolker wrote:
>> Jinsong Zhao yeah.net> writes:
>>
>>
>>> In the following code snippet,
>>> #
>>> a <- strptime("121114 0510", "%m%d%y %H%M")
>>> b <- data.frame(date = a, res = 1:5)
>>> class(a)
>>> class(b[1,1])
>>> #
On 2015/1/5 17:28, Ben Bolker wrote:
Jinsong Zhao yeah.net> writes:
In the following code snippet,
#
a <- strptime("121114 0510", "%m%d%y %H%M")
b <- data.frame(date = a, res = 1:5)
class(a)
class(b[1,1])
#
I am wondering why the class of a and b[1,1] are not the same.
How to make the class
Jinsong Zhao yeah.net> writes:
> In the following code snippet,
> #
> a <- strptime("121114 0510", "%m%d%y %H%M")
> b <- data.frame(date = a, res = 1:5)
> class(a)
> class(b[1,1])
> #
> I am wondering why the class of a and b[1,1] are not the same.
>
> How to make the class of a and b[1,1] to b
Perfect. Thanks.
-Roy
On Jan 5, 2015, at 5:20 PM, Ista Zahn wrote:
> You were close, but a) if you want to map shape to something it needs
> to be inside an aes() call, and b) shape 1 is also an open circle. I
> think this does what you want:
>
> ggplot(temp,aes(x=lon,y=lat)) +
> geom_point(
You were close, but a) if you want to map shape to something it needs
to be inside an aes() call, and b) shape 1 is also an open circle. I
think this does what you want:
ggplot(temp,aes(x=lon,y=lat)) +
geom_point(aes(colour=chlcol, shape=factor(missing)),size=3) +
scale_colour_gradient(low="li
Hi there,
In the following code snippet,
#
a <- strptime("121114 0510", "%m%d%y %H%M")
b <- data.frame(date = a, res = 1:5)
class(a)
class(b[1,1])
#
I am wondering why the class of a and b[1,1] are not the same.
How to make the class of a and b[1,1] to be same?
I would really appreciate you for
QUESTION: Can anyone give me advice on how to resolve an error message
regarding a non positive definite correlation matrix so that the computed
negative cross-variogram can be used to perform co-kriging?
I am performing co-kriging on negatively correlated parameters.
data = dtrd
x y
Hi All:
I have a data frame the first part of which is the following:
date lonlat mean chlcol missing
1 2003-04-23 203.899 19.664 0.0730 0.15176637 0
2 2003-04-24 204.151 19.821NaNNaN 1
3 2003-04-30 203.919 20.351 0.0740 0.15515780
Hello,
I am looking at a data set predicting number of plants (FERO) based on fixed
effects sowing method (Method) and sowing rate (Rate). I also have a slope
random effect for the site and intercept random effects for the nested
variables Block and Strip.ID (to account for experimental design). I
Roger Coppock cox.net> writes:
>
> When will "R" implement the "se.fit" option to the
> predict.nls() function? Is there some schedule?
>
I think this is unlikely to happen, ever (sorry). The exact method
for finding confidence intervals on nonlinear fits would be
to compute likelihood p
Hello,
How can I export the sample data frame below -USvoteS with the associated
probabilities and rules for each row?
library("CHAID")
example("chaid", package = "CHAID")
To be more specific, for each of 1000 rows of the USvoteS data frame I want to
see in which node is that row in (what ar
On Jan 4, 2015, at 10:15 PM, കുഞ്ഞായി kunjaai wrote:
> Dear all,
>
> I am trying to control the font size of correlation axis and gamma label
> (0.6, 0.4 , 0.2 labeling of attached plot) of Taylor diagram in R.
> I am using **plotrix** package for making Taylor diagram.
Nothing attached in w
On Jan 5, 2015, at 2:46 AM, PIKAL Petr wrote:
> Hi
> cc to list, others can have further comments
>
> substitute is superfluous
>
> see results of
>
> substitute(paste(z,sigma)
> paste(z,sigma)
> substitute(paste(z,"sigma")
> paste(z,"sigma")
Perhaps this would be more instructive (after addi
Dear all R members
I have a problem in the R when i wanted to call winbugs as following
#Input data set for WinBUGS
thd18 <-read.table("C:/Users/hp/Desktop/thd18.txt")
#Input data set for WinBUGS
data<-list(N1=2000,N2=2000,P=9,R=Ro,z1=yo1,z2=yo2,thd18)
#Call WinBUGS
model<-bugs (
You can create your own:
http://topepo.github.io/caret/custom_models.html
I put a prototype together. Source this file:
https://github.com/topepo/caret/blob/master/models/files/chaid.R
then try this:
library("CHAID")
### fit tree to subsample
set.seed(290875)
USvoteS <- USvote[sample(1:
Dear Pushpa Methekar,
You can use which.max(), as in
> library(car)
> mod <- lm(prestige~ income + education, data=Duncan)
> which.max(abs(rstudent(mod)))
minister
6
In fact, this is exactly what outlierTest() does, as you can see by looking at
car:::outlierTest.lm
Beyond that, I'd en
When will "R" implement the "se.fit" option to the predict.nls() function? Is
there some schedule?
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://
Hi, this is my function to find rstudent of model which will give me outlier
But I wonder if I could find out the exact no. of outlier in data set.
Like outlierTest() does from car package .
rm.outliers = function(dataset,model){
dataset$rstudent = rstudent(model)
for(i in 1:length(dataset$rstu
Hello,
Is there an option of cross validation for CHAID decision tree? An example of
CHAID is below:
library("CHAID")
example("chaid", package = "CHAID")
How can I use a 10 fold cross-validation for CHAID?
I've read that caret package is to cross-validate on many times of models, but
model CHA
Dear all,
I am trying to control the font size of correlation axis and gamma label
(0.6, 0.4 , 0.2 labeling of attached plot) of Taylor diagram in R.
I am using **plotrix** package for making Taylor diagram.
Thank you all in advance
--
DILEEPKUMAR. R
J R F, IIT DELHI
___
Dear R-help list,
I have a problem regarding text manipulation in R, where my basic knowledge
doesn't suffice anymore. It might be a bigger problem, but any help would be
greatly appreciated and acknowledged.
As input, I have a character string representing some Boolean function, such
as "aB+Bc+D
Hi
cc to list, others can have further comments
substitute is superfluous
see results of
substitute(paste(z,sigma)
paste(z,sigma)
substitute(paste(z,"sigma")
paste(z,"sigma")
Cheers
Petr
> -Original Message-
> From: Levent TERLEMEZ [mailto:lterle...@anadolu.edu.tr]
> Sent: Monday, Janu
On 05 Jan 2015, at 00:21 , Pete Brecknock wrote:
> n <- c(1,2,3,4,5)
> lambda <- c(0.1,0.8,1.2,2.2,4.2)
>
> mapply(function(x,y) rpois(x,y), n, lambda)
Yes. I'd throw in a SIMPLIFY=FALSE to avoid getting results in a different
format if n is constant (then again, sapply() in the original q
Hi
I am not sure what do you want to achieve.
I would change your text function
{text(min(gecici$Tarih)+0.1,y,paste(z,"sigma"),cex=.85,pos=3)}
and instead of copy paste you can open and close appropriate device (pdf, jpg
or emf - from devEMF package)
Cheers
Petr
> -Original Message-
Hi Dear Users,
I abonden ggplot2 adventure for the moment because of urgent results (i will of
course will do this with ggplot2 because I would like to use ggplot2 much more
than standart plots). But this time I came across with a different problem.
Problem is different graphical outputs from sc
28 matches
Mail list logo