On Tue, 25 Mar 2014, Alicia Ellis wrote:
I am running glmulti with the lme function and selection using BIC. When I
do this, I get the following error:
In general, in such situations, it's easiest to use AIC() rather than
BIC() and supply the correct log-number of observations:
On 03/26/2014 05:35 AM, Si Qi L. wrote:
Hi, it's so urgent that I really need your help on R. This following is
part of my data, do u know how to create a dummy coding about it in R? My
control group is Owner occupier. Many thanks for your big help!!!:)
*AppHomeStatus*
Owner occupier
Living
Hello John
As per your suggestion, I have now used the dput function. So my Rdata
dataset bodysize looks like this.
dput(bodysize)
structure(list(neck = c(36.2, 38.5, 34, 37.4, 34.4, 39, 36.4,
37.8, 38.1, 42.1, 38.5, 39.4, 38.4, 39.4, 40.5, 36.4, 38.9, 42.1,
38, 40, 39.1, 41.3, 33.9, 35.5,
Dear r-help mailing list readers,
I am facing a problem using the c.trellis function from the
latticeExtra package, to merge 2 trellis objects (1 levelplot and 1
xyplot).
Using the following example, it works well without customizing the
y-axis of levObj1. But when the y-axis is customized
Hi friends
I would like to rbind two matrices with different column number. How can I
do that? (I dont want to use Plyr package). Is there any way except plyr?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Look at the results of summary(body.pc). You will see that the
first component accounts for 70% of the variability. That is
very large and suggests that the variables are highly correlated
with one another. You could check with cor(bodysize). None of
the correlations is negative and the smallest
This problem is poorly defined, since the meanings of row vectors of different
lengths do not naturally agree. Perhaps you are thinking of binding data frames?
In any event, if you think the behaviour provided by the plyr package is what
you want, why DON'T you want to use it? Further, if plyr
I wrote a simple script to build graphs:
require(RSQLite)
require('igraph')
drv - dbDriver(SQLite)
con - dbConnect(drv, case.db)
dataset - dbGetQuery(con, select * from stream)
table - aggregate(rep(1, nrow(dataset)), by = list(identifier =
dataset$actor_id, name = dataset$actor, page =
Out of curiosity, what's the reason for not using plyr?
If it has to do with installing the package, here's how the rbind.fill
function is implemented in plyr. It may help show you what you'd like to
do:
plyr::rbind.fill
function (...)
{
dfs - list(...)
if (length(dfs) == 0)
Hi,
I have data in a dataframe in following structure
var1 - c(a,b,c,a,b,c,a,b,c)
var2 - c(X,X,X,Y,Y,Y,Z,Z,Z)
var3 - c(1,2,2,5,2,6,7,4,4)
df - data.frame(var1,var2,var3)
Now I'd like to calculate relative values of var3. This values
should be relative to the base value (where var1=c) which is
Hello,
Hopefully there is an answer for this, but I need an ifelse statement that
replaces and returns a value based on a different dataframe. For example:
CurrentDataFrame-ifelse(CurrentDataFrame$Name==NA,match(CurrentDataFrame$Code
with PastDataFrame$Code),replace(CurrentDataFrame$Name) with
I believe this will generalize. But check carefully!
Using your example (Excellent!), use ave():
with(df,ave(seq_along(var1),var2,FUN=function(i)
var3[i]/var3[i][var1[i]==c]))
[1] 0.500 1.000 1.000 0.833 0.333 1.000 1.750
[8] 1.000 1.000
This is kind of a
On 26-03-2014, at 17:09, Johannes Radinger johannesradin...@gmail.com wrote:
Hi,
I have data in a dataframe in following structure
var1 - c(a,b,c,a,b,c,a,b,c)
var2 - c(X,X,X,Y,Y,Y,Z,Z,Z)
var3 - c(1,2,2,5,2,6,7,4,4)
df - data.frame(var1,var2,var3)
Now I'd like to calculate relative
Hello!
I've run SEM using lavaan and after I used summary(myfit) I saw the
following fit indices:
Model Chi Squared
CFI
TLI
RMSEA
SRMR
I was wondering if these are the only fit indices lavaan produces, e.g.:
GFI
AGFI
RMR
Also - does lavaan automatically estimate Chronbach's Alphas for
Please read and in the future follow the Posting Guide, which requests that you
provide a reproducible example... that is, a series of R statements that we can
run to get us to your problem point with a small sample data set that resembles
yours. Forging on anyway...
The ifelse function
On 26 March 2014 10:13, Dimitri Liakhovitski dimitri.liakhovit...@gmail.com
wrote:
Hello!
I've run SEM using lavaan and after I used summary(myfit) I saw the
following fit indices:
Model Chi Squared
CFI
TLI
RMSEA
SRMR
I was wondering if these are the only fit indices lavaan produces,
Great, thanks.
Sure, it's easy to calculate Cronbach's alpha separately...
On Wed, Mar 26, 2014 at 1:36 PM, Jeremy Miles jeremy.mi...@gmail.comwrote:
On 26 March 2014 10:13, Dimitri Liakhovitski
dimitri.liakhovit...@gmail.com wrote:
Hello!
I've run SEM using lavaan and after I used
dplyr's group_by and mutate can create those columns for you:
var1 - c(a,b,c,a,b,c,a,b,c)
var2 - c(X,X,X,Y,Y,Y,Z,Z,Z)
var3 - c(1,2,2,5,2,6,7,4,4)
df - data.frame(var1,var2,var3)
dt - tbl_df(df)
dt %.%
group_by(var2) %.%
mutate(
div = var3[var1 == c],
result_calc = var3/div
)
Hi, I have got a problem with dummy coding and I really couldn't figure it
out. Would you please help me out? this is my codes:
idx-sort(unique(Employment.time$V1));
dummy - matrix(NA, nrow=nrow(Employment.time), ncol= length(idx))
for (i in 1:nrow(Employment.time)) {
for (j in 1:length(idx)) {
Hello
I need to ask questions about the k-means clustering function. Mainly I would
like to know why, with the use of nstart=enough number of times, kmeans always
finds the same clustering arrangements; and this happens even when the input
dataset is sorted in different ways or I take out few
Hello,
I am a student at the University College of London in the Civil Engineering
Department and I need to use the R program for my final year project even if I
do not have any programming background.
I have downloaded the program and I was following the instructions given to me
by my
Hello,
Try, at an R prompt,
install.packages(lmtest, dependencies = TRUE)
install.packages(betareg, dependencies = TRUE)
Hope this helps,
Rui Barradas
Em 26-03-2014 17:44, Campaner, Fabio escreveu:
Hello,
I am a student at the University College of London in the Civil Engineering
No you don't. You have a problem with not learning how R works. Dummy
coding is not needed in R. R uses factors instead. Read An
Introduction to R (ships with R) or a web tutorial of your choice to
learn how to work with R.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650)
On Wed, 26 Mar 2014 18:35:34 + Tomassini, Letizia
tomass...@vetmed.wsu.edu wrote:
Hello
I need to ask questions about the k-means clustering function. Mainly I would
like to know why, with the use of nstart=enough number of times, kmeans
always finds the same clustering arrangements;
On Mar 25, 2014, at 11:02 AM, Michal Kvasnička wrote:
Hallo.
What are the residuals returned by resid() in this code?
library(AER)
m - tobit(y ~ x1 + x2, data=dat)
rr - resid(m2, type = response)
Can I use them to test the normality of the random variable distribution
this way?
I
To add to Ranjan's reply, k-means can potentially find different
results with large nstart= numbers in a large data set. But you
are correct, with a large enough value, the results will be the
same unless there are two solutions that have exactly the same
between sum of squares (unlikely but not
On Mar 25, 2014, at 5:31 PM, Rolf Turner wrote:
On 26/03/14 12:51, David Winsemius wrote:
On Mar 25, 2014, at 9:52 AM, Luo Weijun wrote:
Dear Robert and R project team,
I notice that the Google search function on the R mail list archives page
has stopped working for quite a while,
On Mar 26, 2014, at 4:14 PM, David Winsemius dwinsem...@comcast.net wrote:
On Mar 25, 2014, at 5:31 PM, Rolf Turner wrote:
On 26/03/14 12:51, David Winsemius wrote:
On Mar 25, 2014, at 9:52 AM, Luo Weijun wrote:
Dear Robert and R project team,
I notice that the Google search
On 03/27/2014 05:34 AM, Si Qi L. wrote:
Hi, I have got a problem with dummy coding and I really couldn't figure it
out. Would you please help me out? this is my codes:
idx-sort(unique(Employment.time$V1));
dummy- matrix(NA, nrow=nrow(Employment.time), ncol= length(idx))
for (i in
On Mar 26, 2014, at 11:34 AM, Si Qi L. wrote:
Hi, I have got a problem with dummy coding and I really couldn't figure it
out. Would you please help me out? this is my codes:
idx-sort(unique(Employment.time$V1));
dummy - matrix(NA, nrow=nrow(Employment.time), ncol= length(idx))
for (i in
I have 2 time series (class 'ts') and I want to plot the acf for each and
the joint ccf using acf(ts.union(ts1, ts2)).
According to the ts.union help, ‘ts.union’ pads with ‘NA’s to the total
time coverage but the problem seems to be that one series has a single NA
while the other has 2 NAs.
I would like to understand why the fastclus procedure in SAS is affected by the
initial order of the data. So, with the same dataset, but sorted in a different
way, I get different clusters rearrangements. I find this really disturbing. R
seems to find the stable solution with the use of
Dear useRs,
A similar question has previously been asked by another user
(https://stat.ethz.ch/pipermail/r-sig-geo/2011-September/012791.html) but i'll
try to discuss it from another angle. Its about data reading. I am trying to
read to read a data-set APHRO_MA_050deg_V1101R2.1961.gz from
Hi,
I need to process a data stream using ARIMA model in real-time, that is : I
need to build a ARIMA model and do prediction in every 1 minute using a slide
window with 200 samples .
To improve the performance of prediction, what I plan to do is: 1) to build
ARIMA model using some history
I think rugarch package is doing what you need . Also it uses cluster option to
increase performance :
http://www.unstarched.net/r-examples/rugarch/a-short-introduction-to-the-rugarch-package/
-Original Message-
From: cyl123 [505186...@qq.com]
Date: 03/26/2014 09:43 PM
To: r-help
Please keep the mailing list included by using reply-all... I am not
doing this as a private consultation.
Your sample data is a step forward, but it is still not reproducible. You
could Google R reproducible example and find
This question is about a specific data set you found on the internet... it is
not in any way a question about R. You should be asking the suppliers of this
data, not us.
---
Jeff NewmillerThe
Hello,
Some hints:
- for the year 1961, the total number of values is 27594000,
- there are 180 longitudes and 140 latitudes,
- there are 365 days,
- there are 3 variables,
Compare the total number of values and the result of (180 x 140 x 365 x 3).
The order is precip, rstn, flag,
On Sat, Jul 16, 2011 at 10:26 AM, Roberto Molinari
roberto.molin...@hotmail.it wrote:
Dear R users,
I am currently using the locpoly function from the KernSmooth package to
estimate densities. However, I have some trouble understanding how this
estimation technique is implemented in R. My
39 matches
Mail list logo
