Hello,
I am working with zero inflated models for a current project and I am
getting wildly different results from R's zeroinfl(y ~ x, dist=negbin)
command and Stata's zinb command. Does anyone know why this may be? I find
it odd considering that zeroinfl(y ~ x, dist=poisson) gives identical to
Hi,
source(file,echo=T) will not echo the last comment in file if it is the last
line in the file.
For instance, when sourcing a file containing the following lines
#comment 1
a-1
#comment 2
R will echo
#comment 1
a-1
What is the solution to have R echo all of the comment lines?
Specific
To be fair I don't really understand what the error is telling me so I don't
know how to correct it.
I'm trying to reshape my data so that I do repeated measures ANOVA on it.
The warning message: Error in d[, timevar] - times[1L] : subscript out of
bounds
The program:
library(QRMlib)
Dear R People:
When I have a time series, I can get the frequency of the series via tsp(x)[3].
Is there a similar function for xts and zoo objects, please?
Thanks,
Erin
--
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
Note that the preferred way of doing this with ts series is to use the
frequency() function rather than using tsp. That function also works
for zoo objects.
On Mon, Apr 12, 2010 at 11:20 PM, Erin Hodgess erinm.hodg...@gmail.com wrote:
Dear R People:
When I have a time series, I can get the
Thank you!
Here is something that I just ran into. I have a business day series,
spc, which I obtained using get.hist.quote.
When I tried frequency(spc), I got 1.
str(spc)
‘zoo’ series from 1998-01-02 to 2010-04-09
Data: num [1:3086, 1] 975 977 967 964 956 ...
- attr(*, dimnames)=List of
1 is the correct answer. The difference between successive time
values, deltat, is 1 or a multiple of 1 and frequency is the
reciprocal of deltat.
On Mon, Apr 12, 2010 at 11:46 PM, Erin Hodgess erinm.hodg...@gmail.com wrote:
Thank you!
Here is something that I just ran into. I have a
On Apr 12, 2010, at 9:07 PM, Noah Silverman wrote:
Hi,
I'm trying to setup some complicated binning with statistics and could
use a little help.
I've found the bin2 function from the ash package, but it doesn't do
everything I need. My intention is to copy some of their code and
then
I have produced a series of graphs with the png command, however when I have
finally printed these out the black text appears to have a colour shadow with
blue or red on either side of the letter.
I tried increasing my res to 1200 and this improved it somewhat, but the text
is still not sharp
David,
That helps me a lot. Thanks!!!
-N
On 4/12/10 9:06 PM, David Winsemius wrote:
dat - as.data.frame(matrix( rnorm(200), 100 , 2)) # bivariate normal
n=100
ab - matrix( c(-5,-5,5,5), 2, 2) # interval [-5,5) x [-5,5)
nbin - c( 20, 20) # 400 bins
bins - bin2(dat, ab, nbin) # bin
Hi,
I have data with negative spatial autocorrelation at short lags. It is well
described by a spherical model, but flipped upside from what is normally
seen. It seems that gnls will work with this correlated residual type IF
one parameter in the spherical model can be negative (it is usually
Dear R-helpers
I want to express my respect and thankfulness to all of you.
When I tried to follow a draft book on lme4 authored by Prof. Doug.
Bates (http://lme4.r-forge.r-project.org/book/) to learn how to fit
mixed-effects model, I found I can not execute two functions which
listed in the
Shige Song wrote:
What is the preferred way to get Bayesian analysis results (such as
those from MCMCpacki, MCMCglmm, and DPpackage) into LaTeX table
automatically?
It depends on what you want as output. Let's assume the summary of an
MCMCpack function:
First direct it to a variable,
Dear R-helpers
I want to express my respect and thankfulness to all of you.
I just switch from windows xp to Ubuntu linux platform. I am facing
two problems when I run R in Ubuntu.
I need your helps and advice, thanks in advance.
I almost wasted 3 hours to make lattice functions work. :(. Phew!
I installed package lattice, lattice extra and other supportive packages to
use functions
1.make.groups 2.forplot 3.xyplot.
I also tried local installing. They are installing very well but not
working? Any one have any idea why R
David.Epstein wrote:
I have a data frame with two variables that are factors. One is actually a
TRUE/FALSE factor, and I have coded it as 1/0, a continuous variable, but
I could turn it back into a factor. The second is an ordered factor and
consists of five timepoints. There are several
Dear Gabor,
Thanks for your reply. however:
tail(DJd)
^DJI.Close
2010-04-01 10927.07
2010-04-05 10973.55
2010-04-06 10969.99
2010-04-07 10897.52
2010-04-08 10927.07
*2010-04-09 10997.35*
tail(ag)
2009-11-30 10344.84
2009-12-31 10428.05
2010-01-31 10067.33
2010-02-28
On Apr 12, 2010, at 1:35 AM, newbie_2010 wrote:
I almost wasted 3 hours to make lattice functions work. :(. Phew!
I installed package lattice, lattice extra and other supportive
packages to
use functions
I thought lattice was now part of the default installation? It's not
lattice
Bernd Dittmann wrote:
(Sample boot code see below)
...
So far so good. But how can I bring it into table-form whereby the
combination of all 10 different weights (setting steps as 1) and their
respective bootstrapped means are printed:
weights mean.boot
1
Always provide sessionInfo() and traceback() when reporting errors. /Henrik
On Mon, Apr 12, 2010 at 8:34 AM, Mao Jianfeng jianfeng@gmail.com wrote:
Dear R-helpers
I want to express my respect and thankfulness to all of you.
I just switch from windows xp to Ubuntu linux platform. I am
require(lattice)
Loading required package: lattice
Error in library(package, lib.loc = lib.loc, character.only = TRUE,
logical.return = TRUE, :
'lattice' is not a valid installed package
library(lattice)
Error in library(lattice) : 'lattice' is not a valid installed package
I got
Dear Iurie,
I want to run Shapiro-Wilk test for each variable in my dataset, each
grouped by variable groupFactor.
Note that, at least on a single dependent variable with a grouping variable, a
possible simplification may arise when homogeneity of variances is assumed and
reasonable. You may
Hi
r-help-boun...@r-project.org napsal dne 12.04.2010 07:14:14:
David Winsemius wrote:
I am guessing that the first time through when i= 5200 that i+1 is
indexing an entry that does not exist. What does str( Price[[1]]
[5200+1] ) return? What about str(Ca)?
So what is
On Apr 12, 2010, at 2:20 AM, newbie_2010 wrote:
require(lattice)
Loading required package: lattice
Error in library(package, lib.loc = lib.loc, character.only = TRUE,
logical.return = TRUE, :
'lattice' is not a valid installed package
library(lattice)
Error in library(lattice) :
hello,
i ran a glmer for binomial data, with one (ordered) between and one within
subjects factor. now i was thinking of how to plot my dependent.
i plotted the probablities (p=X(incidents)/n(observations)) for each group
of the between factor in a line chart, with its corresponding standard
Dear R-helpers
I want to express my respect and thankfulness to all of you. I have
ever ask this question about one hour ago. Following Mr. Henrik
Bengtsson's advice I re-submit this query here.
I just switch from windows xp to Ubuntu linux platform. I am facingtwo
problems when I run R in
Dear all,
I am a relative novice with R, so please forgive any terrible errors...
I am working with a GLM that describes a response variable as a function of
a categorical variable with three levels and a continuous variable. These
two predictor variables are believed to interact.
An example of
Dear Users:
I am able to read the file precip.nc (CMAP monthly precip) using the Rclim
package. How ever, when I try to write the data in a netcdf file, I get the
following error about the missing value attribute:
pcpmon-netcdfread(precip.nc,lon,lat,time,precip,unpack=T)
I'm fitting a regression tree with rpart and I want to write the values
for every leaf in a dataset. As an example take the variable turnover.
Let's suppose my tree for turnover has 30 leaves and I want to have 30
datasets with dataset 1 containing the turnover values of the units in
leaf 1,
Dear Matthew,
The easiest way the get the estimates (and their standard error) for the
different slopes it to reparametrise your model. Use resp ~ var1 : cat +
0 instead of resp ~ var1 * cat
HTH,
Thierry
ir. Thierry
On 2010-04-12 3:15, Matthew Carroll wrote:
Dear all,
I am a relative novice with R, so please forgive any terrible errors...
I am working with a GLM that describes a response variable as a function of
a categorical variable with three levels and a continuous variable. These
two predictor
Petr Pikal wrote:
You did not tell much more about your data and procedures. Each object
type has some distinct way of indexing and you can not mix them up.
x-1:10
x[5]
[1] 5
x[[5]]
[1] 5
x[5,]
Error in x[5, ] : incorrect number of dimensions
x-list(1:10)
x[1]
[[1]]
[1]
I'd like to use a string to refer to an R object with the end
objective of going through a loop and saving various files of the same
name with different contents using a numbered suffix.
# This will be the loop counter and file suffix.
master.i - 1
# This is the generic file name.
unislopes -
On 2010-04-11 14:27, Ajay Shah wrote:
Folks,
I have a strange situation where:
library(MatchIt)
f- function(d) {
m- matchit(treatment ~ lsales + major.industry,
data=d, method=nearest, discard=hull.treat)
treatmentfirms- match.data(m, group=treat)
save(list=paste(unislopes, master.i, sep=),
file=paste(unislopes,master.i,.Rdata,sep=))
On Mon, Apr 12, 2010 at 12:06 PM, David Young dyo...@telefonica.net wrote:
I'd like to use a string to refer to an R object with the end
objective of going through a loop and saving various files of the same
?get
On 2010-04-12 5:06, David Young wrote:
I'd like to use a string to refer to an R object with the end
objective of going through a loop and saving various files of the same
name with different contents using a numbered suffix.
# This will be the loop counter and file suffix.
master.i- 1
#
Dear Dieter,
That sounds like a good solution, thanks.
Shige
On Mon, Apr 12, 2010 at 2:27 AM, Dieter Menne
dieter.me...@menne-biomed.de wrote:
Shige Song wrote:
What is the preferred way to get Bayesian analysis results (such as
those from MCMCpacki, MCMCglmm, and DPpackage) into LaTeX
Dear all,
After I run package lme4 and function glmer, I got the error message as the
following.
Matrix is not symmetric under lme4
Is it possiable to avoid this error when my data is huge? Or use the other
package? Thank you.
[[alternative HTML version deleted]]
Hello List,
I am new to R and statistics in general.
My two colleges and I would like to get some training in the proper
usage of R.
A quick internet search provided no results for the region of Basel in
Switzerland.
Does some one on this list knows of any trainers or schools in the
region
I assume what you want is not the last data value in each month
labelled with the nth day of the month (or the last day of the month
if there are fewer than n days) but the last day in the data that is
on or prior to the nth of the month. If that is what you want then
try this.
nth.of.month -
Dear all,
Using rGraphviz, I would like to make my edge labels a smaller size than
the node labels, but I am unable to change the font size of edge labels
independently of node labels. The edge labels always seem to use the
font size defined for nodes.
I have searched the archives and found a
Hi
r-help-boun...@r-project.org napsal dne 12.04.2010 12:51:24:
Petr Pikal wrote:
snip
Nobody except you has your data available, so without providing more
clues
you can not expect mor relevant answers.
Try str(your.objects) and maybe you could use debug to see how
Dear Dominik,
Mango Solutions is a long-established R training and consulting company based
in the UK. In February we started Mango Solutions AG based right here in
Basel, Switzerland.
We have an Introduction to R course in Basel on May 17th to 18th. Here's a
list of the other courses
-- begin inclusion --
I'm fitting a regression tree with rpart and I want to write the values
for every leaf in a dataset. As an example take the variable turnover.
Let's suppose my tree for turnover has 30 leaves and I want to have 30
datasets with dataset 1 containing the turnover values of
Dear R users,
I have a question on multilevel modeling with R when there is endogenous
regressor(s) involved. As far as my economics background concerned, I
understand that Hausman-Taylor estimator (via Generalized Least Square)
deals with this situation and the package plm does the trick.
Ben, Ravi, Chuck, and Haneef,
Note that the standard normal density can be written (in pseudo-TeX) as
1/sqrt(2*pi) * exp[- ( x/sqrt(2) )^2].
The exponential on the right is a special case of the stretched
exponential function
exp[- (x/sqrt(tau))^beta].
The stretched exponential function has
Hello,
I'm using dmnorm from the package {mnormt} and getting strange results.
First, according to the documentation, dmnorm should return a vector of
densities, and I'm only getting one value returned (which is what I would
expect). I've been interpreting this as the joint density of all values
Dear all
i want to fit the self start model in nls. i have two question. i have a
function,
(asfr ~ I(((a*b)/c))+ ((c/age)^3/2)+ exp((-b^2)*(c/age)+(age/c)-2)
i am wondering how to build the selfstart model. there is lost of example,
(i.e. SSgompertz, SSmicman, SSweibull, etc). my question is, how
Hi,
I am trying to get figures from multiple source files in Sweave. My test
file is as follows
\documentclass{article}
\usepackage{Sweave}
\begin{document}
\begin{section}{notitle}
This is a simple Sweave test
gethypergraphs_rhea,fig=TRUE,echo=F=
source(./testfig.r)
@
\\
End of the simple
Quoth Martin Morgan on Boomtime, the 29th of Discord:
It seems like it would be hard to think about the tasks that are
being executed, how many processes there are, how load balancing
works, etc.
A flat representation of our nested data is, alas, too large to
contain in memory; and dynamically
Hi David,
Thanks! I just tried it with new sets version (1.0-6) and it works great!
Thanks again for a very nice package.
Best regards,
Ryszard
--
Confidentiality Notice: This message is private and may contain
On 12/04/2010 12:08 PM, michele donato wrote:
Hi,
I am trying to get figures from multiple source files in Sweave. My test
file is as follows
\documentclass{article}
\usepackage{Sweave}
\begin{document}
\begin{section}{notitle}
This is a simple Sweave test
gethypergraphs_rhea,fig=TRUE,echo=F=
Dear List members,
its me again, fighting with lattice graphics. I am trying to plot a world map,
add some points on different locations with different colors and add a legend,
but did not succeed yet with the legend. Here is my code:
library(fields)
# Data for demonstration
data_x =
There us no groups argument in your xyplot call, so how is auto.key
supposed to define a legend? (and note that auto.key should be a logical not
a list).
Please re-read the auto.key section in the xyplot man page.
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
Hello all,
I am trying to save a multiband matrix in a tiff file with writeGDAL() from
the rgdal package.
The code below only saves the first band of the matrix StockMat.
Do you have any suggestion to correct my code?
Thanks a lot.
nbRow=nrow(StockMat)
nbCol=ncol(StockMat)
Hi R enthusiasts,
I am dealing with logging data from different sources that contain
data from user activities. The data is all timelined with one column
containing Epoch time and two columns containing data (x and y
coordinates of mouse movements) = three columns for each source. I
have up to 10
I have an external process that is listening to and responding on a
port. This is working fine. If the external service is not running,
though, read.socket does not return and I do not see a way to
interrupt it. Is there a general way to deal with this situation to
have a timeout or of testing
hello,
i'd appreciate help with my glmer.
i have a dependent which is an index (MH.index) ranging from 0-1. this index
can also be considered as a propability. as i have a fixed factor (stage)
and a nested random factor (site) i tried to model with glmer. i read that
it's possible to use a
Dear Charles,
this sounds perfect.
I knew I would find the information on this list.
Thanks to all that answered my question.
Regards,
Dominik Riva
Universitätsspital Basel
Med. Querschnittsfunktionen
Pathologie
Schoenbeinstrasse 40
CH-4031 Basel
Telefon 061 556 53 67
Cordless
I have searched and tried to read before posting but can find nothing to
accomplish change Excel dates in double format to R
Can someone please help
I have a vector of double like this from Excel.
39965.0004549653
and I want to put them in R such that I can display them in any Date and
Time
If one of your variable is discrete,such as 1/0, the logistic regression is
better.Try the glm() function.
--
View this message in context:
http://n4.nabble.com/simple-question-about-contrasts-lm-and-factors-tp1835964p1837091.html
Sent from the R help mailing list archive at Nabble.com.
I mean I want to convert them to whatever is the standard R DateTime class.
--
View this message in context:
http://n4.nabble.com/Excel-date-to-R-format-tp1837208p1837280.html
Sent from the R help mailing list archive at Nabble.com.
__
Dear useRs,
I am glad to announce that that my package mvngGrAd (read: moving
grid) is now available on CRAN.
The package implements moving grid adjustment, which is a spatial method
used in (unreplicated) plant breeding trials to adjust phenotypic values
for environmental effects. The
Actually, Bert, auto.key can be used. See below.
On 2010-04-12 11:41, Bert Gunter wrote:
There us no groups argument in your xyplot call, so how is auto.key
supposed to define a legend? (and note that auto.key should be a logical not
a list).
Please re-read the auto.key section in the xyplot
Noel, thanks a lot. This will help me someday. But I have a question.
When we run Shapiro-Wilk test, the homogenity of variances is a
mandatory condition?
2010/4/12 Yvonnick Noel yvonnick.n...@uhb.fr:
Dear Iurie,
I want to run Shapiro-Wilk test for each variable in my dataset, each
grouped
I am really new with R Graphical user interfacefunctions. I am developing a
software package to calculate pKa (biochemistry) but I want to make it look
aesthetically pleasing and make it user friendly. I have heard that R has
some GUI (Graphical user interface) and you can do some really cool
Dear R-listers
I am writing a manuscript for a scientific journal in clinical medicine.
I have three groups of patients, and I present a 10*3 table of their
characteristics in Table 1. Some of their characteristics, e.g. their
age, are on a continuous scale, others are dichotomous. I am
I have toyed with 'traitr' recently... It is not fancy, but relatively
easy to figure out.
-k
On Mon, Apr 12, 2010 at 4:00 PM, Amitoj S. Chopra amit...@gmail.com wrote:
I am really new with R Graphical user interfacefunctions. I am developing a
software package to calculate pKa (biochemistry)
Dne Po 12. dubna 2010 22:00:21 Amitoj S. Chopra napsal(a):
I am really new with R Graphical user interfacefunctions. I am developing a
software package to calculate pKa (biochemistry) but I want to make it look
aesthetically pleasing and make it user friendly. I have heard that R has
some GUI
Hello,
I am writing a manuscript for a scientific journal in clinical medicine.
I have three groups of patients, and I present a 10*3 table of their
characteristics in Table 1. Some of their characteristics, e.g. their
age, are on a continuous scale, others are dichotomous. I am thinking of
Hi
I like gWidgets, you can choose different api, like rgtk2,
tcltkeasy to use in R, looks pretty(at least for me), try the
examples in the website first, the link below is the answers with
tutorial links, you can also check gWidgets vignette.
I am trying to create a simple window that displays a dataframe. I was able
to put together a function from a couple of examples on the web. For the
most part my function works. The only problem is with resizing the window.
I found a post where someone suggested using
See the relevant article in R News 4/1.
On Mon, Apr 12, 2010 at 10:36 AM, ManInMoon xmoon2...@googlemail.com wrote:
I have searched and tried to read before posting but can find nothing to
accomplish change Excel dates in double format to R
Can someone please help
I have a vector of double
Hi,
On 12 April 2010 22:07, Peter Jepsen p...@dce.au.dk wrote:
3. Are there R packages that can draw tables?
the gplots package has a textplot() function, and the gridExtra
package a tableGrob(),
http://rwiki.sciviews.org/doku.php?id=tips:graphics-grid:table
In theory it should be possible
Gabor Grothendieck ggrothendieck at gmail.com writes:
See the relevant article in R News 4/1.
On Mon, Apr 12, 2010 at 10:36 AM, ManInMoon wrote:
I have a vector of double like this from Excel.
39965.0004549653
and I want to put them in R such that I can display them in any Date
Try this:
sprintf(%.3f, as.numeric(as.POSIXct(39965.004,origin=as.Date(1970-1-1
[1] 39965.004
getOption(digits.secs)
NULL
as.POSIXct(39965.004,origin=as.Date(1970-1-1))
[1] 1970-01-01 06:06:05 EST
sprintf(%.3f, as.numeric(as.POSIXct(39965.004,origin=as.Date(1970-1-1
[1] 39965.004
Is there a setting or easy way I can remove the code that makes R only listen
for web traffic on 127.0.0.1? I've noticed it seems to reject traffic to both
the regular help files (i.e. what starts up when I do ?function), as well as
if I try to connect to pages I've made for the R.rsp package.
Hello all, I'm trying to transform data frames by grouping the rows by the
values in a particular column, ordered by another column, then picking the
first row in each group.
I'd like to convert a data frame like this:
x y z
1 10 20
1 11 19
2 12 18
4 13 17
into one with three rows, like this,
afoo wrote:
Hi,
I am trying to SSH to a remote server through R script. In other words, I
would like to know how I can get a SSH connection to the remote server and
then execute commands on that server with the R script.
So in bash, I would normally type ssh -lusername
Sorry about the double post, but this line:
Sharpie wrote:
ssh u...@webserver.com -e cd ~/.ssh;cat webKey.pub authorized_keys
Should be:
ssh u...@webserver.com cd ~/.ssh;cat webKey.pub authorized_keys
I.e., omit the -e flag.
Apologies,
-Charlie
-
Charlie Sharpsteen
James -
If I understand you correctly:
getone = function(df)df[order(df$x,df$y),][1,]
describes what you want from each data frame corresponding
to a unique value of x. Then, supposing that your data frame
is called df:
sdf = split(df,df$x)
will create a list of data frames for the
Tena koe James
You might try duplicated(), or more to the point !duplicated()
orderedData[!duplicated(orderedData$x),]
HTH
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of James Kebinger
Sent: Tuesday, 13
Hi,
I'm trying to setup some complicated binning with statistics and could
use a little help.
I've found the bin2 function from the ash package, but it doesn't do
everything I need. My intention is to copy some of their code and then
modify as needed.
I have a vector of two columns:
83 matches
Mail list logo