Please read ?lm! -- where it says:
method:
the method to be used; for fitting, currently only method = "qr" is
supported; method = "model.frame" returns the model frame (the same as
with model = TRUE, see below).
More to the point, your request for a "spearman" method for lm() makes
little or no
Hi,
A following function was kindly provided by GGally’s maintainer, Barret
Schloerke.
function(data, mapping, ...) {
p <- ggplot(data = data, mapping = mapping) +
geom_point(color = I("blue")) +
geom_smooth(method = "lm", color = I("black"), ...) +
theme_blank() +
> On Apr 28, 2016, at 9:32 AM, BISWAJIT KAR wrote:
>
> Respected all,
>I am working on a socio-economic survey (named as
> National Sample Survey in India provided by National Sample Survey
> Organization, Govt. of India) data of individual as well as households.
> This is a
> On Apr 28, 2016, at 4:51 AM, Gordon, Fabiana
> wrote:
>
>
> Maybe I wasn't clear about my query.
>
> I'm very familiar with pre-allocation and vectorization and I had already
> wrote an R code for this problem in this way. My question wasn't about the
> most efficient way to solve the
Respected all,
I am working on a socio-economic survey (named as
National Sample Survey in India provided by National Sample Survey
Organization, Govt. of India) data of individual as well as households.
This is a sample survey where stratified random sapling method has been
use
What's not possible, as far as I know, is to create a variable using an
expression that refers to elements within the variable.
For example, suppose that a some point during my R session, there is no
variable named tmpx:
> exists('tmpx')
[1] FALSE
If I try to reference tmpx, in this case by usin
We are pleased to announce the following book:
Title: Beginner's Guide to Zero-Inflated Models with R
Authors: Zuur, Ieno
Book website: http://www.highstat.com/BGZIM.htm
Paperback or EBook can be order (exclusively) from:
http://www.highstat.com/bookorder.htm
TOC: http://www.highstat.com/BGS/ZI
Hello,
I'm trying to run a code that uses the snowfall package. Here is the structure
of my code.
sfInit(parallel=T, cpus = 5, slaveOutfile="ErrorMessage.txt")
sfExportAll()
sfLibrary(rgdal)
sfLibrary(raster)
sfLibrary(sp)
sfLibrary(rgeos)
sfLibrary(snowfall)
system.time( sfClusterA
Thanks, adding asp = 1 did the trick!
-Original Message-
From: David L Carlson [mailto:dcarl...@tamu.edu]
Sent: Thursday, April 28, 2016 12:17 PM
To: Greg Snow <538...@gmail.com>; Simon, Heather
Cc: r-help@r-project.org
Subject: RE: [R] polygon angle option perpendicular to axis
The ang
The angle is not based on the polygon edges, but it can seem that way if you do
not use the asp=1 argument in your plot. Try this example,
> plot(1:10, 1:10, type='n')
> x <- c(1, 3, 5, 3)
> y <- c(3, 5, 3, 1)
> polygon(x, y, angle=0, density=10)
> polygon(x, y + 4.5, angle=45, density=10)
> poly
I think you are seeing the Moire (interference) pattern arising from
the interaction of lty=4 (dotted/dashed lines) and closely spaced
parallel lines. Use lty=1 (solid lines) to see this.
[BTW, your script did not work because you changed the capitalization
of top and bottom halfway through it.]
Hi All,
please find enclosed the missing attachment.
Kind regards
Georg
-- cut --
#-[ Header ]
--
# Program : Framework for R scripts
# Author: Georg Maubach
# Date : 2016-03-03
# Update: 2016-04-27
> On Apr 27, 2016, at 11:37 PM, WU Qingwei wrote:
>
> hi all,
>
>
> We have developed a software called HCBA (Highly Consumable Business
> Analytics) using R as modelling engine.
> We now want to donate the whole software (R serves only as modelling engine
> in this project. There are other
Maybe I wasn't clear about my query.
I'm very familiar with pre-allocation and vectorization and I had already wrote
an R code for this problem in this way. My question wasn't about the most
efficient way to solve the problem. It was about whether in R it was possible
to use the same index
If you want to include random intercepts in a model fit by bam/gam, then
include
s(g,bs="re")
in the model formula, where g is a factor variable with one level for
each group requiring a random intercept.
Now suppose that for each level of group g you want a random slope
w.r.t. x. You should
R Users,
The d/p/q/r functions for the bridge distribution are now available in
bridgedist.
When a random intercept follows the bridge distribution, as detailed in
Wang and Louis (2003) http://dx.doi.org/10.1093/biomet/90.4.765>>, a marginalized
random-intercept logistic regression will still be
Hi all,
2 days ago a new R packages appeared on CRAN
https://cran.r-project.org/web/packages/RZabbix/index.html
RZabbix is an R interface to the 'Zabbix API' data <
https://www.zabbix.com/documentation/3.0/manual/api/reference>. Enables
easy and direct communication with 'Zabbix API' from 'R'.
Y
Thanks for the suggestions, all of you!
I first began to think about using somehow permutations of the
gtools-package. I will continue utilizing Peter's solution.
My purpose is to divide basic musical rhythm units (whole note, half
note, quarter note etc durations) to meaningful entities in al
hi all,
We have developed a software called HCBA (Highly Consumable Business Analytics)
using R as modelling engine.
We now want to donate the whole software (R serves only as modelling engine in
this project. There are other core functions.) to the R community and
open-source it. But we don't
Thanks for the question. Here is a sample of the code for my plot:
Top = c(34, 39, 42, 45, 46, 41, 41, 40, 43, 38, 33, 33)
Bottom = c(24, 29, 32, 36, 32, 34, 32,41, 40, 39, 29, 24)
plot(1,1, col = "white", xlim = c(1.3,11.7), ylim = c(0,80), axis =
FALSE, xaxt = "n")
Greetings, R users!
It is with pleasure that I am announcing the release of msmtools package on
CRAN:
https://cran.r-project.org/web/packages/msmtools/index.html
msmtools provides a fast and general method for restructuring classical
longitudinal data into augmented ones. The reason for this i
Filling polygons with lines is a throwback to the time when the height
of quality graphics was the mechanical pen plotter (a device that used
a pen in a mechanical arm to draw the plot on a piece of paper).
Computing and printing technology has advanced quite a bit from that
day, so you may want to
Dear all,
I’ll need your help with obtaining robust clustered errors. I use polr command
in MASS package m<–porl(y~x1+x2,data=mydata, method=probit). In the rms
package, this is as simple as: clusterSE<–robcov(m, mydata$id). Is it possible
to do something similar for polr object as well? Thank
Give us reproducible data using dput(). There are easier ways to accomplish
what you want in R. The overwriting problem stems from the fact that you have a
double loop so each iteration of the outer loop overwrites the data produced in
the inner loop.
Also your if/else specification should be
Hi
?factor
help page says rather cryptically
The encoding of the vector happens as follows. First all the values in exclude
are removed from levels. If x[i] equals levels[j], then the i-th element of the
result is j. If no match is found for x[i] in levels (which will happen for
excluded value
Thx a lot Gabor!
Aljosa Aleksandrovic, FRM, CAIA
Quantitative Analyst - Convertibles
aljosa.aleksandro...@man.com
Tel +41 55 417 76 03
Man Investments (CH) AG
Huobstrasse 3 | 8808 Pfäffikon SZ | Switzerland
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Se
Hi Gabor,
Thanks a lot for your help!
I tried to implement your nonlinear least squares solver on my data set. I was
just wondering about the argument start. If I would like to force all my
coefficients to be inside an interval, let’s say, between 0 and 1, what kind of
starting values are norm
The nls2 package can be used to get starting values.
On Thu, Apr 28, 2016 at 8:42 AM, Aleksandrovic, Aljosa (Pfaeffikon)
wrote:
> Hi Gabor,
>
> Thanks a lot for your help!
>
> I tried to implement your nonlinear least squares solver on my data set. I
> was just wondering about the argument start
Hi
your initial ds
> str(ds)
'data.frame': 2 obs. of 3 variables:
$ var1: num 1 1
$ var2: logi TRUE FALSE
$ var3: logi NA NA
first result
> str(ds)
'data.frame': 2 obs. of 6 variables:
$ var1 : num 1 1
$ var2 : logi TRUE FALSE
$ var3 : logi N
Sorry
these
T&NA = T (you can decide that regardless value in NA the result must be T)
F&NA = NA (you cannot decide hence NA)
should be
T | NA = T (you can decide that regardless value in NA the result must be T)
F | NA = NA (you cannot decide hence NA)
Cheers
Petr
> -Original Message-
Hi
Your script is not reproducible.
Creating Check_U_0__Kd_1_2011 from Umsatz_2011 and Kunde01_2011
Error in ifelse(Kunden01[[Umsatz]] == 0 & Kunden01[[Kunde]] == 1, 1, 0) :
object 'Kunden01' not found
>
This is interesting
x <- c(NA, FALSE, TRUE)
names(x) <- as.character(x)
outer(x, x, "&") #
Hi All,
my script tries to do the following on factors:
> ## Check for case 3: Umsatz = 0 & Kunde = 1
> for (year in 2011:2015) {
+ Umsatz <- paste0("Umsatz_", year)
+ Kunde <- paste0("Kunde01_", year)
+ Check <- paste0("Check_U_0__Kd_1_", year)
+
+ cat('Creating', Check, 'from', Umsatz,
32 matches
Mail list logo
