Hi,
I am using RStudio as my R editor.
After someday I accidentally hit something, the whole program is
repeated in the Console whenever I compile it. How can I fix it so that the
whole program won't be repeated in the future?
Thanks,
miao
[[alternative HTML version deleted]]
Dear R community
I am using the gee package to run logistics regression on paired cases from
a panel sample.
We are getting request from a reviewer to use sample weights to account for
different level of attrition.
After searching the documentation I am unable to find a way to incorporate
Hi,
This looks like an rstudio question that's not specific to R itself
There are support forums at rstudio.org which you should ask your question
on, ani think you'll get better help there
On Monday, April 30, 2012, jpm miao wrote:
Hi,
I am using RStudio as my R editor.
After
On 29.04.2012 12:01, ya wrote:
Hi Uwe and Indrajit,
Thank you very much for the response.
Since we are talking about this package, can I ask more about how to use
it to deal with statistical issues please?
I have this data with categorical missing values in it. Now I am trying
to impute
Hi Uwe,
Thank you very much for the suggestions.
ya
ya
From: Uwe Ligges
Date: 2012-04-30 15:18
To: ya
CC: Indrajit Sengupta; r-help@r-project.org
Subject: Re: [R] parallel package
On 29.04.2012 12:01, ya wrote:
Hi Uwe and Indrajit,
Thank you very much for the response.
Since we are
On 29.04.2012 16:28, Kyriakos Antoniou wrote:
Hi,
i am trying to run an ANCOVA and a bootstrapped ANCOVA analysis on a specific
data set. I am using the ancova and ancboot functions as in the following code:
setwd(C:/Users/User/Desktop/Rdatabilingualstudy2012)
Dear useRs,
Just to let you know, the results using mgcv 1.7-13 are OK. @Simon:
could you perhaps still let me know if the procedure I apply w.r.t.
assigning the labels to random slopes and intercepts is correct?
With kind regards,
Martijn
--
***
Martijn
Well, I would suggest using the code already in place in the survival
package. Here is my code for your problem.
I'm using a copy of the larynx data as found from the web resources for
the Klein and Moeschberger book.
larynx - read.table(larynx.dat, skip=12,
Hi
I have a distance matrix which is computed by user defined method. I
would like to plot the dendrogram. I would like to use different color
and want the leaves laying down bottom.
The script like this. I am not familiar with R. I followed the example
shown in
Hi,
I have created 2 models in logistic regression and got the predictor values
are significance on output. Here are 2 summaries of 2 model. HOw can we
compare 2 models by what factor or coefficient and say which model is best
Please help
Hello,
i have a question calculating the shaped area between the two curves (see
image).
http://r.789695.n4.nabble.com/file/n4597813/test.png
I try to use a Simpson-Integral but it doesn't work. R doesn't know the
command!
It is possible, that i need a library?
Curve 1: y
Curve 2: z
Hi!
I'm trying to do a lm() test on three objects. My problem is that R protests
and says that the variable lengths differ for one of the objects
(Sweden.GDP.gap). But I have double checked that the number of observations
are the same. All three objects should contain 9 observations but R only
Goodmorning everybody,
i'm an italian statistician and i'm using R for research.
Could someone tell me some indices to see the goodness of fit in multilevel
modelling?
I'm using the lmer function, and I want to know if my model fit well my
data.
I actually want to justify the use of
See the function geeglm() in package:geepack
Maurice
O. Maurice Haynes, Ph.D.
Statistician, Child and Family Research Section
Laboratory for Comparative Ethology
Eunice Kennedy Shriver
National Institute of Child Health and Human Development
NIH, DHHS
6705 Rockledge Drive, Suite 8030
Bethesda,
Hi Mark,
I'm sorry I forgot about this. I found that Gls was out of date with
respect to gls. It's now fixed and you can
source('http://biostat.mc.vanderbilt.edu/tmp/Gls.s') to override the old
version with the fixed version until the next release of rms. You won't
need the workaround below.
Sorry!! Data now attached. Note that this is a .dat file while in the R code
i import data directly from an SPSS file (.sav).
Best,
Kyriakos
Date: Mon, 30 Apr 2012 15:20:04 +0200
From: lig...@statistik.tu-dortmund.de
To: antonio...@hotmail.com
CC: r-help@r-project.org
Subject: Re: [R]
Dear Uwe.
Many thanks for your reply. I attach the data set that i use for the analyses.
Hopefully this will help you identify the problem. I'm really new in using r
so please bear with me.
All the best,
Kyriakos
Date: Mon, 30 Apr 2012 15:20:04 +0200
From:
Hi!
I'm trying to do a lm() test on three objects. My problem is that R
protests
and says that the variable lengths differ for one of the objects
(Sweden.GDP.gap). But I have double checked that the number of
observations
You shall check it again. BTW How did you checked it?
are the
On 30.04.2012 15:52, Kyriakos Antoniou wrote:
Sorry!! Data now attached. Note that this is a .dat file while in the R code
i import data directly from an SPSS file (.sav).
Oh, come on, we do not have endless time to import your data, so
inlcuding code please.
Moreover, you have not
Hi again,
Sorry for connecting you again(I am still getting this error:subscripts out
of bounds). the problem is that I am trying to apply R script for one file
to many files. It succeeded to do that for the first file but not for the
others.I think the problem in my code is that ( d ) (in the
Dear Santos,
Since your models are nested you can apply a likelihood ratio test.
M0 - glm(formula = Spend_bucket ~ Freq + Address_is_res +
last_update_days_ago, data = qdataset, family = binomial)
M1 - glm(formula = Spend_bucket ~ Freq + Address_is_res, data = qdataset,
family = binomial)
On Sun, 29 Apr 2012, Abigail Clifton wrote:
Hi,
This is my code (my data is attached):
library(languageR)
library(rms)
library(party)
OLDDATA - read.csv(/Users/Abigail/Documents/OldData250412.csv)
OLDDATA$YD - factor(OLDDATA$YD, label=c(Yes, No))
OLDDATA$ND - factor(OLDDATA$ND, label=c(Yes,
Dear Uwe,
really sorry about this.
The code i use is the following (data attached again):
setwd(C:/Users/User/Desktop/Rdatabilingualstudy2012)
bilingualdataforconferences2012-read.delim(bilingualdataforconferences2012.dat,
header=TRUE)
As far as I understand you want:
lmo - lm(Simon ~ TypeBiling + FAS, data = bilingualdataforconferences2012)
anova(lmo)
summary(lmo)
for the ordinary ancova.
Not sure if there is a proper implementation for that kind of
bootstrapped robust ancova you are using. You may have to fix the code
Missing comma between -1.180620213 and -0.525964648
Saint bahman.za...@comhem.se wrote in message
news:1335779012548-4597768.p...@n4.nabble.com...
Hi!
I'm trying to do a lm() test on three objects. My problem is that R
protests
and says that the variable lengths differ for one of the
Yes, But what is your opinion about d in
sprintf(C:\\Users\\aalyaari\\Desktop\\New folder (6)\\Yar_%00d.bin, d),
because will go from one to the number of rows(some files have 248 and
others have 224) and my results will change accordingly. for example i will
get from the first file 31 files
As you have not provided a clue of what your models are, one can only guess.
But if you mean using lots of fixed effects vs a random effect, the answer is
that there is no such animal. They are two different non-nested models, and
should be chosen based on subject matter considerations.
Hi Chiara,
If you just want to compare model fit, you could use a LRT between
models where you do and do not estimate the variance/covariance matrix
of random effects.
R^2 in mixed models do not have the same nice properties they do in
fixed effects models.
Cheers,
Josh
On Mon, Apr 30, 2012
Dear useRs,
I would like to be able to smoothly use the log base graphics argument as an
argument to pairs(). Currently in my hands, pairs() accepts the argument,
makes transformed plots accordingly, but then stumbles over some annotation
commands. Example:
begin session
Dear colleagues,
I have a question regarding controlling the jitter when plotting
predictions in the rms package. Below I've simulated some data that
reflect what I'm working with. The model predicts a continuous variable
with an ordinal score, a two-level group, and a continuous covariate. Of
You can't just decide a command exists and try to use it by your
arbitrarily chosen name did you look at the function integrate()
which does exist? It uses a smarter algorithm than Simpson's rule.
Alternatively, Simpson's rule is very easy to code with vectorization
-- you should be able to
On Mon, 30 Apr 2012, Santiago Guallar wrote:
Hello,
I'm having problems working with date values in POSIXct format.
Indeed you are.
Here is what I got (eg.lig attached):
x - read.table(eg.txt, sep = ',', col.names=c(ok,time,secs,lig)) #
it gives time as factor
z -
A quick Google suggests that sintegral is in the Bolstad2 package:
http://cran.r-project.org/web/packages/Bolstad2/
R. Michael Weylandt michael.weyla...@gmail.com wrote in message
news:CAAmySGOCtNxWNHJDka=68criphxnppwap6fm7ufkhtx7xsz...@mail.gmail.com...
You can't just decide a command exists
New version RExcel 3.2.7-1
Download wrapper for new version of RExcel.
___
R-packages mailing list
r-packa...@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-packages
__
R-help@r-project.org mailing
If you want to know if your model fit will your data, then looking at
residual-type plots is useful, as can be plotting the model-predictions and
observed data together. You might also find interesting R2 statistics for
mixed models by Matthew Kramer. Beware--as others have indicated, there
is
hc.obj is not a dendrogram.
plot(dend1, xlab='',ylab='',sub='', nodePar=list(pch = c(1,NA), cex=0.8,
lab.cex = 0.8))
See ?plot.hclust and ?plot.dendrogram. They are not using the same
arguments.
Regards,
Eloi
On 12-04-30 06:36 AM, farscape2012 wrote:
Hi
I have a distance matrix which
I am trying to make a nice 2x1 plot and add a timestamp with comment.
The pstamp function from Hmisc works nicely when mfcol=c(1,1),
but when mfcol=c(2,1), the stamp winds up in the wrong place:
require('Hmisc')
opar - par(mfcol=c(2,1))
plot(1:10)
title(main=MAIN Title)
plot(1:20)
Thank.
I made a stupid mistake. Now it is working.
On 4/30/2012 8:30 PM, Mercier Eloi wrote:
hc.obj is not a dendrogram.
plot(dend1, xlab='',ylab='',sub='', nodePar=list(pch = c(1,NA),
cex=0.8, lab.cex = 0.8))
See ?plot.hclust and ?plot.dendrogram. They are not using the same
arguments.
This question might better be asked on r-sig-db or r-sig-mac (I am not
sure which one would better, but please not both).
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 4/29/12 7:08 AM, julia.jacob...@arcor.de
A couple of things that I did not see mentioned by the others:
Generally statistics plots work better in .png files than in .jpg
files due to the type of compression each uses (detailed image plots
and some surface plots may be the exception), though that may be out
of your control (some journals
length(Sweden.GDP.gap.adjust)
[1] 8
Sweden.GDP.gap.adjust
[1] 0.673792123 1.196706756 1.196131539 0.646944002 -0.312886525
[6] -1.706584861 -0.369401194 -0.003280389
R is reading Sweden.GDP.gapp.adjust as having 8 elements. Why I don't know.
John Kane
Kingston ON Canada
-Original
On Mon, Apr 30, 2012 at 10:05 PM, Stefán Hrafn Jónsson shj...@gmail.com wrote:
Dear R community
I am using the gee package to run logistics regression on paired cases from
a panel sample.
We are getting request from a reviewer to use sample weights to account for
different level of
Because of the missing comma.
On Mon, Apr 30, 2012 at 4:56 PM, John Kane jrkrid...@inbox.com wrote:
length(Sweden.GDP.gap.adjust)
[1] 8
Sweden.GDP.gap.adjust
[1] 0.673792123 1.196706756 1.196131539 0.646944002 -0.312886525
[6] -1.706584861 -0.369401194 -0.003280389
R is reading
The sunrise and sunset are calculated for each time value on the input vector,
but all times are treated in the TZ timezone. I can see that the Canary Islands
are in the GMT timezone, so the example TZ is right.
I can see that you entered the latitude and longitude backward (longitude is
x),
Hi all,
I was wondering if you can help me with the following situation:
I have a data frame that includes weather station data for 30 years in
the form:
YEAR, MONTH, DAY, TEMP
1970, 01, 01, -15
...
1999, 12, 31, -21
I would like to add another variable JULIAN that assigns the integers
1 to
Perhaps this is easier: convert your data to Date objects and use the
yday (year day) function from the lubridate package:
i.e., something like this:
jd - apply(DAT, 1, function(x) as.Date(paste(x, collapse = ),
format = %Y %m %d)
# Just a guess you might have to debug it a little
yday(jd)
It
I don't really work with dates but thought I'd pass a solution on. I think
that there some great packages for handling dates though (lubridate) and you
may want to convert your data to a true date instead of separate columns.
# FUNCTION TO INDEX DATES
date.int - function(month, year, day){
Hey Michael and Tyler!
Thanks for the input and the effort! The lubridate package is really awesome!
What I did now is just this here:
dat$DATE-as.Date(paste(dat$YEAR,dat$MONTH,dat$DAY,sep='-'),format='%Y-%m-%d')
require(lubridate)
dat$JULIAN-yday(dat$DATE)
That's it. Nice and easy :)
Hello, again.
Try to use the values in 'listfile' directly:
# Wrong!!!
for(h in length(listfile)) # Just last value ?!?
{b=file.info(listfile[h])$size/67420/4
wind - file(listfile[h]), rb) # Parenthesis do not match.
# Right
for(h in listfile)
{b=file.info(h)$size/67420/4
wind - file(h, rb)
Dear All
I'm having trouble working out what exactly loess means by its Standard Error
of the Residuals denoted s
and in particular when the weights argument is invoked.
For example, if the weights are weights are all =1, then s^2 is nearly sum sq
res/ (n -1 - 'equiv num paras')
If the
farscape2012 farscape2...@live.com writes:
Hi
I have a distance matrix which is computed by user defined method. I
would like to plot the dendrogram. I would like to use different color
and want the leaves laying down bottom.
The script like this. I am not familiar with R. I followed the
This posting is for the help-file record in that the problem has already
been resolved via a misplaced discussion in the R-bug forum (Bug 14894). I
place it here so make it easier for others to find. The problem involved an
error in calling Rscript from Matlab on Mac OS X 10.6 Snow Leopard.
The
Hey,
i am trying to do the MLE for Garch and have a problem with the optim
function.
Initally i tried optim with Method=BFGS. Reading trhough the forum i found
out
i would neet bounds. So i went on with Method=L-BFGS-B.
But now my parameters equal the lower bounds.
out - optim(par=initial,
Hi,
I need to perform a robust regression with several dummy independent
variables (i.e., binary), which are paired (clustered). The outcome
(dependent variable) is continuous.
I have tried using the robcov function with the next syntaxis:
out1 - ols(y ~ x1 + x2 + x3 , x=T)
robustclust -
Hello,
I have a list of dates, such as
dates- as.Date(c(1996-03-29,1996-05-30,1996-09-28,1996-05-09))
from which I would like to extract the week number for each date, with week
n°1 being the week going from Dec 30th 1995 to Jan 6th 1996 (1995-30-12 to
1996-06-01). Any suggestion for a simple
On 04/28/2012 02:32 PM, Jim Lemon [via R] wrote:
On 04/28/2012 04:12 AM, bodiless wrote:
Dear R-users
I hope someone could help me on this problem.
I want to create a multiple kiteChart showing the real values with a
scalebar on each indicating the scale .
Here are some
Hello to all,
I'm new to R so I have a lot of problems with it, but I'll only ask the main
one.
I have clustered an environmental matrix with 2 different methods, and I'd
like to plot them in a PCA and a db-RDA. I mean, I want see these clusters
in the plots like points of differents colours,
Dear all,
I'm pleased to announce the release on CRAN of the mmand package
(for Mathematical Morphology in Any Number of Dimensions). It provides
functions for performing mathematical morphology (erode, dilate,
etc.), smoothing, and other kernel-based operations on array-like
objects of any
Thanks a lot! One gets blind sometimes.
--
View this message in context:
http://r.789695.n4.nabble.com/Different-varable-lengths-tp4597768p4598762.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Trying to do a regression for four variables. I get information only for two.
The rest is marked NA. Why? And what does NA mean?
Thanks!
Sweden.infl.dev - c(0.2, 0.6, -1, -0.2, -1, -1.5, -0.7, 0.7, -0.1, 0.3)
Sweden.GDP.gap - c(-0.024662769, -0.01519859, -0.027251109, -0.027129556,
Hi everyone.
I want to transpose a data frame. Lets say the following DF:
df = data.frame(matrix(ncol=4, nrow = 10))
df[,1] = c(rep(1,5),rep(2,5))
df[,2] = c(rep('a',4),rep('b',3),rep('c',3))
df[,3] = c(letters[5:14])
df[,4] = runif(10)
I would like to form a data frame with each line
Hello,
I was wondering whether there had been any progress on this bug yet?
I have checked for updates to the package and it doesn't seem that there
have been any, but I might be mistaken!
Similarly to Mark, I am getting the following error message when running a
Gls with two (factor) covariates
Hello,
Try
#
# Adapted from several posts to R-Help,
# mostly by David Winsemius, also by Gabor Grothendiek.
#
nweek - function(x, format=%Y-%m-%d, origin){
if(missing(origin)){
as.integer(format(strptime(x, format=format), %W))
}else{
x -
Hello,
Filoche wrote
Hi everyone.
I want to transpose a data frame. Lets say the following DF:
df = data.frame(matrix(ncol=4, nrow = 10))
df[,1] = c(rep(1,5),rep(2,4), 3)
df[,2] = c(rep('a',4),rep('b',3),rep('c',3))
df[,3] = c(letters[c(5:13,13)])
df[,4] = runif(10)
I would
I'm interested in this question, too, not so I can hide code, but so I can
run R code faster. It's been my experience that compiled code always runs
faster than interpreted code. Can you explain further how to compile a
front-end?
--
View this message in context:
Thank you sire, this is exactly what I was looking for.
Many thanks,
Phil
--
View this message in context:
http://r.789695.n4.nabble.com/Matrix-transposition-tp4599219p4599310.html
Sent from the R help mailing list archive at Nabble.com.
__
Hello,
I am doing a simple linear regression analysis that includes few variables.
I am using a bootstrap analysis to obtain the variation of my variables to
replacement.
I am trying to obtain the coefficients 95% confidence interval from the
bootstrap procedure.
Here is my script for the
Hello fellow R users,
I am trying to extract the coefficient values during a bootstrap operation.
Here is the list of my variables that I would like to extract the
coefficient values from:
(Intercept)
LogRds_25k
GeoRockbimodal
Hi Rui,
Thanks for your help!!
It works perfectly, but when I call stor.confint, I obtain the list of
confidence interval (10,000 times). Is there an easy way to summarize the
results and getting only one 2.5% and 97.5% values for each variable?
Thanks again for your help
Steeve
--
View this
Hello,
I am having difficulty putting a frequency table on my clipboard in R
studio. I do not know what the error message is. I was thinking it might
have to do with changing a setting on my clipboard but I do not know how to
do this. I will copy and paste the exact command and error message:
Thanks a lot, Jeff
Yes, I applied colsplit in package reshape. Sorry, I forgot to state it in the
piece of code I attached.
Actually, these are archival data from a geolocator attached to a seabird. The
origin lies on Gran Canaria and data were gathered all along its migratory
journey for one
Hi,
Is it possible to read, say 100 dat files in one time? I want to combine these
100 data set into one single dat file.
Any suggestions?
Thank you very much.
ya
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
You could try
legend(bottomright, legend = Sys.time())
to get the timestamp in the right place (or alternatively, using text()). The
documentation for pstamp mentions that it has a problem under Windows when
mfrow is used.
The mar argument allows you to adjust the plot margins. See ?par.
I am creating a package, and a rather long usage line for one function becomes
truncated in the generated .pdf manual. Is there a way that I can prevent that
truncation?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Hello,
I have a set of 100 variables with 1560 observations. I did an O.L.S
regression of three of these variables on a fourth. But there are problems
of endogeneity... So I look in my dataset for instruments to do an IV. I
can't find a good instrument because their correlation with my
Hello,
baconbeach wrote
Hello,
I am doing a simple linear regression analysis that includes few
variables. I am using a bootstrap analysis to obtain the variation of my
variables to replacement.
I am trying to obtain the coefficients 95% confidence interval from the
bootstrap
Hello,
baconbeach wrote
Hi Rui,
Thanks for your help!!
It works perfectly, but when I call stor.confint, I obtain the list of
confidence interval (10,000 times). Is there an easy way to summarize the
results and getting only one 2.5% and 97.5% values for each variable?
Thanks
First let me say thank you for all the help everyone gave me.
Below is the code I came up with to look at Rainfall pH reading for
Kentucky and Tennessee.
Also there is a more *generic* ver of the code on my web page if anybody
ever needs it in the future.
Thanks again for your swift response!!
With your last line, I get
rowMeans(sapply(stor.confint, colMeans))
2.5 %97.5 %
0.3256882 0.4604677
I need the values (2.5% and 97.5%) for each variable of my model. I don't
think this what I am getting.
This is what my script looks like now,
On Mon, Apr 30, 2012 at 12:34 PM, Zeeberg, Barry (NIH/NCI) [E]
zeebe...@mail.nih.gov wrote:
I am creating a package, and a rather long usage line for one function
becomes truncated in the generated .pdf manual. Is there a way that I can
prevent that truncation?
How long is rather long? In
Whatever your post means, it's about statistics, not R. Post to a
statistics list like stats.stackexchange.com instead.
-- Bert
On Mon, Apr 30, 2012 at 3:24 PM, kebrab67 selamgetac...@gmail.com wrote:
Hello,
I have a set of 100 variables with 1560 observations. I did an O.L.S
regression of
81 matches
Mail list logo