Hi,
I am trying to run R summary command through JRI to get the result for
mulitvariate Linear Regression
eg. result - lm(Performance Score ~ Department+Grade,data = StudentData)
summary(result)
on running the above cmd using in R Console will fetch me below result:
Call:
lm(formula
On Jun 23, 2015, at 3:20 PM, boB Rudis wrote:
You can do something like:
aaa - function(data, w=w) {
if (class(w) %in% c(integer, numeric, double)) {
I think you will find that inherits(w, numeric) is more compact and safer.
Both integer and double do inherit from numeric (and double is
Steve Taylor steve.tay...@aut.ac.nz
on Wed, 24 Jun 2015 00:56:26 + writes:
Note that objects can have more than one class, in which case your == and
%in% might not work as expected.
Better to use inherits().
cheers,
Steve
Yes indeed, as Steve said, really do!
I don't use JRI, but the data seem to be there. If you are looking for the row
names, try ?rownames.
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics:
Dear members,
I would like to conduct nonparametric two way ANOVA (repeated
measures) and post hoc test by using ARTool and lsmeans package.
I have the following data,
x-runif(120)
A-gl(3,40,labels=c(a1,a2,a3))
B-gl(2,20,120,labels=c(b1,b2))
ID-gl(20,1,120)
data-data.frame(ID,x,A,B)
Then,
On 24 Jun 2015, at 03:28 , gianni lavaredo gianni.lavar...@gmail.com wrote:
I am doing an Repeated Measures ANOVA and the Bonferroni post hoc test for
my data using R project. The ANOVA gives a significantly difference between
the data but not the Bonferroni post hoc test.
anova(aov2)
Does the package semPlot not do what you want? I notice that you got
an error when you used library(semPlot) because you don't have all the
dependencies installed. When you have semPlot working, it should be
able to produce a graphical output of the results, including
'covariance arrows'.
David
On 24/06/2015 7:08 AM, Martin Batholdy via R-help wrote:
Hi,
I would like to define the size for tick-marks, axis-titles, legends, drawing
symbols etc. absolute,
meaning that regardless of the size of the plot device, the font-size /
character size is the same.
Thus if I output my plot
I have a large spatiotemporal database (131 spatial locations, 9 years of
daily weather data) in long table format which I would like to convert to
an ST object.
head(tempmean)
site latlon alt var year mth
day value date doy numdate
518941
I am interested in using quantile regression to fit the following model at
different quantiles of a response variable:
(1) y = b0 + b1*g1 + b2*g2 + B*Z
where b0 is an intercept, g1 and g2 are dummy variables for 2 of 3
independent groups, and Z is a matrix of covariates to be adjusted for in
Hi,
That's the default, isn't it?
I am sorry – one of my plots was actually set up with mfrow.
But the documentation actually explains the change in cex when using mfrow;
In a layout with exactly two rows and columns the base value of cex is
reduced by a factor of 0.83: if there are three or
Have you considered using the semPlot package? It works nicely with
lavaan models (among other sem packages). There is also the DiagrammeR
package.
Rick
On 06/23/2015 10:48 AM, DzR wrote:
Dear Senior users of R/R Studio,
I am very new to this environment hence am unable to plot the SEM
Another option would be to use segments instead of lines.
library(lubridate)
temp$end - temp$time + minutes(temp$duration)
library(ggplot2)
ggplot(temp, aes(x=time, xend = end, y=Typ, yend = Typ, colour=person)) +
geom_segment()
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek /
Hi Petr,
I'm not exactly sure this is what you are looking for, but try:
start_indices-which(!is.na(temp$duration))
pp_gantt_info-list(
labels=paste(as.character(temp$person[start_indices]),
temp$Akce,temp$Typ,sep=-),
starts=temp$time[start_indices],
ends=temp$time[start_indices+1],
Hi,
I would like to define the size for tick-marks, axis-titles, legends, drawing
symbols etc. absolute,
meaning that regardless of the size of the plot device, the font-size /
character size is the same.
Thus if I output my plot with pdf(width=5, height=5) or pdf(width=15,
height=15), the
Bert, can you be more specific about which article for those of us who don't
subscribe?
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics: ##.#.
Andrew Gelman's Working Through Some Issues and the two Letters to
the Editor that follow responding to the editorial decision to ban P
values from The Journal of Basic and Applied Social Psychology (BASP).
You may wish also to read ASA President's David Morgenstern's
reflexive and entirely
On 6/23/15 4:07 PM, My List wrote:
All,
I am new to the vcd package and new to R too.
Welcome to R and glad you found the vcd package.
1) I have a lickert analysis based data set.
2) I am doing a hypothesis tests on the variables ( like, is there a
relationship between the choice of a Doctor
On Jun 24, 2015, at 10:52 AM, Chel Hee Lee wrote:
Could you also kindly check the following links in the rw-FAQ manual at
http://cran.r-project.org/bin/windows/base/rw-FAQ.html?? The links list in
the below seem to be broken. I hope these links are fixed in the very near
future.
You may want to consider another way of getting your answer that takes
advantage of some of R's features:
# Make some example data
cods - LETTERS[1:10] # Ten companies
yrs - 2010:2014 # 5 years
set.seed(42) # Set random seed so we all get the same values
# Chances of revenue for a given
Giacomo,
Please include some representative data. It is not clear why your offset of 4
(z$cod[i - 4]) is going to be an accurate surrogate for complete data.
Since I do not have your data set or its true structure I am having to guess.
# make 5 copies of 200 companies
companies -
I used the dplyr do function to apply a kernel regression smoother to a 3
column data table (grouping index, x, y) with about 7 M rows and 45000
groups.
This runs quickly, about 1-2 minutes.
It creates an data table (44,326 by 2) - grouping index, kernel smoothing
output.
The kernel smoothing
Maybe something like the punch cards on github?
https://github.com/hadley/ggplot2/graphs/punch-card
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie Kwaliteitszorg / team Biometrics Quality Assurance
Kliniekstraat 25
1070
Not an answer to your question, but you should not be using dummy
variables in R. Use factors instead. Please read a R tutorial or text
-- there are many -- to learn how to fit models in R. You might also
wish to consult a local statistician or post on a statistics list like
Hi Thierry
Thanks a lot. This is the option I found later.
Still not 100% satisfactory. Maybe somebody did similar task and will have
different opinion how to visualise such data.
Jim’s gant diagrams are worth consideration so I will try to elaborate it
further.
Cheers
Petr
From: Thierry
Hi Jim
Thanks a lot.
gantt.chart is worth trying, beside Thierry's segment solution. I need to think
it over if it can be better for visualisation.
Cheers
Petr
-Original Message-
From: Jim Lemon [mailto:drjimle...@gmail.com]
Sent: Wednesday, June 24, 2015 12:50 PM
To: PIKAL Petr
Hi
I am not sure what you want
plotMatrix
, , group1
a b c d
1 1 0 0 0
2 1 0 0 0
3 1 1 0 0
4 0 1 0 0
5 0 1 1 0
, , group2
a b c d
1 0 0 1 0
2 0 0 1 1
3 0 0 0 1
4 NA NA NA NA
5 NA NA NA NA
If you do not want to show the NA's without giving them a different colour then
The easiest way is to find out what the actual projection is, and
reconstruct a sensible representation of the grid.
Usually you have to guess, but it's possible to figure out. You *can* plot
by building a proper mesh in longlat, but my preference is full rescue. You
need to know the projection
Thanks for the suggestions, Gunter.
On Wed, Jun 24, 2015 at 10:33 AM, Bert Gunter bgunter.4...@gmail.com wrote:
Not an answer to your question, but you should not be using dummy
variables in R. Use factors instead. Please read a R tutorial or text
-- there are many -- to learn how to fit
The Details section of ?par starts of with:
Each device has its own set of graphical parameters.
(So this is not Mac-specific.)
Strictly speaking, the options you set with par() are not reset when you
open a new graphics device. Rather, when a new device is opened, it is
initialized with
Hi Edwin,
If you look at the build output you will notice that the C++11 compiler
flag is not being used. I just created a small package using Rcpp11 and
your function and it worked without a problem. I can't give you a specific
reason without seeing your package but there are some
Hi all,
I've just started using Rcpp and am trying to get cpp11 support working. As
suggested I added [[Rcpp:plugins(cpp11)]] to my source file and a test
function:
// [[Rcpp::export]]
int useCpp11() {
auto x = 10;
return x;
}
This works fine when using:
sourceCpp(filename)
from R, but I
In R lattice, I am trying to label predefined groups of rows in a
matrix of data with strips. Currently, the length of the strips fail
to match the different sizes of the groups as the data representation
only allows groups with the same size.
One possibility to solve this might be to suppress
Hi, I have a dataset (728 obs) containing three variables code of a company,
year and revenue. Some companies have a complete history of 5 years, others
have not a complete history (for instance observations for three or four
years).I would like to determine the companies with a complete
i tried the semPlot but it flopped. various other packages also did not
perform. i will try the DiagrammeR package and revert.
meanwhile, i tried going to the onyx package and its neat, though i have
yet to spend some time on it to familiarise myself with the nuts and bolts
of the process.
Could you also kindly check the following links in the rw-FAQ manual at
http://cran.r-project.org/bin/windows/base/rw-FAQ.html?? The links
list in the below seem to be broken. I hope these links are fixed in
the very near future.
Under the section 2.4 Can I customize the installation?
*
Glad to help,
The SystemRequirements is for a package. I believe the example in the
gallery is intended to demonstrate a function where if you set the
CXX_FLAGS with:
Sys.setenv(PKG_CXXFLAGS=-std=c++11)
And then compiled a single *.cpp file with Rcpp::sourceCpp(test.cpp)
I believe it should
Could you kindly check if the following links are working fine in the
R-FAQ page at http://cran.r-project.org/doc/FAQ/R-FAQ.html? The links
listed in the below seem to be broken. I hope these links are fixed in
the very near future.
Under the section 2.6 Are there Unix-like binaries for R?,
I would **strongly** recommend that you speak with a local statistical
expert before proceeding further. Your obsession with statistical
significance is very dangerous. (see the current issue of SIGNIFICANCE
for some explanation).
Cheers,
Bert
Bert Gunter
Data is not information. Information is
Please repost your question in plain text rather than HTML - you can
see below that your code got rather mangled. Please also include some
sample data using dput() - made-up data of similar form is fine, but
it's very hard to answer a question based on guessing what the data
look like.
Sarah
On
Comments below
On 24/06/2015 19:26, giacomo begnis wrote:
Hi, I have a dataset (728 obs) containing three variables code of a company,
year and revenue. Some companies have a complete history of 5 years, others
have not a complete history (for instance observations for three or four
Thank you! I was missing the SystemRequirements. I guess it could be useful
to add this to the example given here:
http://gallery.rcpp.org/articles/simple-lambda-func-c++11/
Cheers, Edwin
On Wed, 24 Jun 2015 at 17:50 Charles Determan cdeterma...@gmail.com wrote:
Hi Edwin,
If you look at the
Thank you, Thierry. And yes, Bert, it turns out that it is more of a
statistical question after all, but again, since my question used specific R
functions, R experts are well placed to help me.
As pairewise.t.test was recommended in a few tutorials about repeated-measure
Anovas, I assumed it
I have spatial data from the WRF-model. I don't know what kind of projection
was used to produce these data (neither longitude nor latitude are constant
at any borders), but I have 2D matrixes for each longitude and latitude. The
area covered is North America. How can I add a map using these 2D
Hello,
I am trying to add a legend to my PCA plot so that it looks neat. I think
plotting this outside of the chart area would be good but I cannot seem to
fathom the correct code for this. I wondered if anyone could help please?
The code I am using is as follows:
grp- with(Matan,
Hi all,
I use ireadLines function to iterate large '.csv' files from '.zip' file.
When I execute the nextElem function in R console, I can only fetch the first
line of file content no matter how many times . Example code show as bellow.
library(iterators)
Hello,
I have few scripts that have been written in MATLAB. I need to
translate or convert them into R. They all deal with reading in a netcdf
file and doing some plots. I managed to read in the netcdf file with these
API calls
# Read input file
input_dir = /home/aan/aa/data/r
file =
Tena koe Samantha
You probably need to set some graphics parameters such as xpd and mar (see
?par), and then give the × and y location of the legend rather than
'bottomright' (see ?legend).
HTH
Peter Alspach
PS Please don't post in html (see the posting guide) ... P
-Original
Edwin van Leeuwen edwinvanl at gmail.com writes:
Thank you! I was missing the SystemRequirements. I guess it could be useful
to add this to the example given here:
http://gallery.rcpp.org/articles/simple-lambda-func-c++11/
No.
If you actually read the Rcpp documentation--eg the Rcpp
49 matches
Mail list logo