Perhaps what you want to explore is the object returned from the summary method
rather than the fit object?
coef(summary(kinfit))
---
Jeff NewmillerThe . . Go Live...
DCN:Bas
I don't really know what you want, but if you have many columns with the same
names I am wondering why this is so. Do you really want to merge, which puts
all of the non-key columns side-by-side in one data frame? If so, why don't you
start by renaming the columns so they will make sense in the
Hi, R experts
I am currently using lmekin() function in coxme package to fit a
mixed effect model for family based genetic data. How can I extract the p
value from a lmekin object? When I print the object in R console, I can
see the p value and Z value are just over there. But I can not ex
Hello all,
Evidently my previous message met some filter due to subject line. I am
re-sending my message. I apologize if this was sent out twice.
Based on "Ripley & Thompson, Analyst, 1987
", I am trying to do a regression of my data which assumes a linear
relationship between measurement
Hi,
I couldn't find any attached data. Could you dput() the data?
A.K.
- Original Message -
From: HJ YAN
To: r-help@r-project.org
Cc:
Sent: Wednesday, September 5, 2012 7:57 PM
Subject: [R] Summarizing data containing data/time information (as factor)
Dear R user
I want to create a
On Sep 5, 2012, at 8:51 PM, Sam Steingold wrote:
> I have a list of data frames:
>
>> str(data)
> List of 4
> $ :'data.frame': 700773 obs. of 3 variables:
> ..$ V1: chr [1:700773] "200130446465779" "200070050127778" "200030633708779"
> "200010587002779" ...
> ..$ V2: int [1:700773] 0 0
I have a list of data frames:
> str(data)
List of 4
$ :'data.frame': 700773 obs. of 3 variables:
..$ V1: chr [1:700773] "200130446465779" "200070050127778" "200030633708779"
"200010587002779" ...
..$ V2: int [1:700773] 0 0 0 0 0 0 0 0 0 0 ...
..$ V3: num [1:700773] 1 1 1 1 1 ...
$
On Sep 5, 2012, at 4:57 PM, HJ YAN wrote:
> Dear R user
>
> I want to create a table (as below) to summarize the attached data
> (Test.csv, which can be read into R by using 'read.csv(Test.csv, header=F)'
> ),
Unfortunately you did not read theinformation about posting attachments
carefully en
There are several packages to do this, with probably rgdal as the most
general for reading - it will read many image formats as a
SpatialGridDataFrame (class from package sp). This class behaves like
a data.frame in many ways and can be coerced to one with
as.data.frame. Non-spatial image formats w
On Sep 5, 2012, at 1:58 PM, mahout user wrote:
> Hello dear,
>
>
> I am new to R, Have developed the model for prediction. I dont know
> exactly about the followed terms
> residual standard error
> degrees of freedom,
> multiple R-squared,
> adjusted R-squared
> F-statistics
> p-values
>
> T
Dear R user
I want to create a table (as below) to summarize the attached data
(Test.csv, which can be read into R by using 'read.csv(Test.csv, header=F)'
), to indicate the day that there are any data available, e.g.value=1 if
there are any data available for that day, otherwise value=0.
Hello all,
Based on "Ripley & Thompson, Analyst, 1987 ", I am trying to do a regression of
my data which assumes a linear relationship between measurements by two
modalities of the same physiological parameter. The complication is that my
errors are heterogeneous, i.e. not only both X & Y var
Hi,
Not sure how your data looks like. With the sample data below, the code works.
Try this:
set.seed(1)
dat1<-data.frame(MIGRATION=sample(c(0,1),100,replace=TRUE),distance=sample(40:80,100,replace=TRUE))
RR.rebuild<-glm(MIGRATION~distance,data=dat1,subset=!(1:100 %in%
c(56,23,20,9,19)),family=b
Hello dear,
I am new to R, Have developed the model for prediction. I dont know
exactly about the followed terms
residual standard error
degrees of freedom,
multiple R-squared,
adjusted R-squared
F-statistics
p-values
Thanks in advance.
[[alternative HTML version deleted]]
I'm interested in using R to perform some statistical analysis of bitmaps.
When I search, I see a lot of information about outputting bitmaps, but I'm not
finding much about loading a bitmap into a data frame so that it can be
analyzed. Would someone have a quick pointer to help me out?
On Sep 5, 2012, at 2:05 PM, Andrew Yee wrote:
> The answer to this may be obvious, but I was wondering in the rms
> package and the survfit(), how you can plot the censored time points
> as ticks.
survfit() is not in rms.
>
> Take for example,
>
> library(survival)
> library(rms)
>
> foo <-
You will probably have to contact the maintainer of package sqldf (G.
Grothendieck ) as it appears that the package does
not support missing values. I tried stripping out the NAs so that the fields
were empty, but the results are the same. Fields specified as numeric are
converted to character. Ins
On Wed, Sep 5, 2012 at 5:31 PM, David Reiner wrote:
> I'm trying to use sqdf's function read.csv.sql to read CSV files in which the
> missing values are represented by NA's.
> Plain old read.csv works fine on these files, but they are rather large and
> I'd like to filter using sql-like statemen
I see the attachment didn't go through.
here is the text:
2012-07-01,15:50:00,NA,0,NA,0,NA
2012-07-01,15:51:00,NA,0,NA,0,NA
2012-07-01,15:52:00,NA,0,NA,0,NA
2012-07-01,15:53:00,NA,0,NA,0,NA
2012-07-01,15:54:00,NA,0,NA,0,NA
2012-07-01,15:55:00,NA,0,NA,0,NA
2012-07-01,15:56:00,NA,0,NA,0,NA
2012-07-0
I'm trying to use sqdf's function read.csv.sql to read CSV files in which the
missing values are represented by NA's.
Plain old read.csv works fine on these files, but they are rather large and I'd
like to filter using sql-like statements.
However, even if I specify field.types correctly and nrow
The answer to this may be obvious, but I was wondering in the rms
package and the survfit(), how you can plot the censored time points
as ticks.
Take for example,
library(survival)
library(rms)
foo <- data.frame(Time=c(1,2,3,4,5,6,10), Status=c(1,1,0,0,1,1,1))
answer <- survfit(Surv(foo$Time, fo
Hallo folks,
I know I should not ask the same question again. But I have a problem I cannot
solve and without the solution I am stuck and lost, unable to get along with my
work!
Someone suggested I should try the code below in order to eliminate the
outliers from my data. I did as I was told
Hi,
Does anyone know how to upload a table to green plumb & have it be
distributed?
I know how to upload using sqlSave(chann, d, tablename="castaneg.wh_d",
rownames=F, colnames=T)
but how can I make my table be distributed randomly on the server.
In SAS you can use the option "distribute_on=rand
I can't replicate your problem. I created a data set with "Male" and
"Female" since that is what you indicate, but your commands use "M" and "F"
which is different. When I use "Male" and "Female" the recoding is just as
expected, but you don't even need to do this. You probably already have a
facto
Hi Conrad,
On Wed, Sep 5, 2012 at 3:14 PM, Conradsb wrote:
> I currently have a data set in which gender is inputed as "Male" and "Female"
> , and I'm trying to convert this into "1" and "0".
This is usually not necessary, and makes things more confusing. "Male"
and "Female" is clear and self-ex
This is very off-topic here.
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Pla
pdb wrote:
>> as.POSIXlt('2004-10-31 02:00:00')
> [1] "2004-10-31"
>> as.POSIXlt('2004-10-31 03:00:00')
> [1] "2004-10-31 03:00:00"
>
> How do I tell as.POSIXlt() to ignore daylight savings and just convert to
> a time as is? I've read about the 'isdst' but it is still unclear what to
> do.
Ther
I currently have a data set in which gender is inputed as "Male" and "Female"
, and I'm trying to convert this into "1" and "0".
I found a website which reccomended using two commands:
data$scode[data$sex=="M"] <- "1"
data$scode[data$sex=="F"] <- "2"
to convert to numbers, and:
data$scode <- fa
April Lindeman wrote on 09/04/2012 01:43:06 PM:
>
> I am trying to compare Vbert growth curves from several years of
> fish data. I am following the code provided by: http://
> www.ncfaculty.net/dogle/fishR/gnrlex/VonBertalanffy/VonBertalanffy.pdf
> . Specifically the section on VBGM Comparisons
Hello all,
Can someone please tell me how I can simulate a key press (particularly
arrow keys and space bar) on a Mac?
Any help is much appreciated.
Cheers
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.et
Hello everybody,
I have a problem with the model.table of anova.
I have data (datos2) with 4 cluster (V7), calculated with daisy and hclust.
str(datos2)
'data.frame': 56 obs. of 7 variables:
$ Estacion: Factor w/ 56 levels "Abradelo","AltoDoRodicio",..: 1 2 3 4 5 6
7 8 9 10 ...
$ Invierno:
Maybe I should have added further details to my previous message.
So, after running
dd <- dredge (M.gam.auto)
avgmod <- model.avg(dd, subset = delta < 2)
I ask for the confidence intervals using
confint (avgmod)
What you are suggesting would result in the visualization of the initial
model,
On Sep 5, 2012, at 8:06 AM, Nicolai Schneider wrote:
> Hello everyone,
>
> a problem with the plot.zoo function. In the parameters of the function,
> cex.lab is ignored. I tried to reduce the size of the yaxis labels by at
> least 50%.
>
> --
>
> Example:
>
> sample <- as.zoo(
Hi,
If i understand what you mean, you can do
Plot(M.gam.auto)
You Will get the centered predictions for smooths, you can also extract
predictions for each term using (predict(M.gam.auto, type="terms")).
The results you got are thé coefficients for each of you smooths pieces (as
k=4, you get 3)
Hi,
It is not clear to me what a model-averaged (sensu MuMIn package) GAM is...
How is the smoother averaged and how can I see it ?
My model is:
M.gam.auto <- gam (Leksite ~ s(Altmean,k=4) + s(Class_unfav_prop, k=4) +
+ s(Visib2000, k=4) + s(Ddisturbance, k=4) + s(Steppe_prop,
k=
a) This sounds like homework. This is not a homework support forum.
b) If it is not homework, you should take one or more classes on statistics.
Your questions are more about theory than R and this is not a statistics theory
mailing list.
c) You ask questions about the use of your data, but you
These questions are off topic for this list. Try a statistical list
like stats.stackexchange.com.
Probably better yet, as your statistical skills sound like they are
somewhat limited, consult a local statistician for help.
-- Bert
On Wed, Sep 5, 2012 at 7:54 AM, Vignesh Prajapati wrote:
>
> Hel
A horribly pedestrian alternative to plotrix:
barplot(x,beside=TRUE, ylim=c(0,90), etc.)
par('usr')# just getting plot limits
[1] 0.44 15.56 0.00 90.00
par(new=T)
lines(c(0.44,15.56),c(20,20))
lines(c(0.44,15.56),c(40,40))
etc.
[[alternative HTML version deleted]]
___
On 09/05/2012 02:41 AM, Hermann Norpois wrote:
To whom it may concern.
As I would like to analyse some array data I was keen on downloading the
lumi package that depends obviously on hdrcde that is not available for r
2.12.1. I did not find instructions to solve or circumvent this problem.
Insta
Hello folks,
I am on learning phase of R. I have developed Regression Model over six
predictor variables. while development, i found my all data are not very
linear. So, may because of this the prediction of my model is not exact.
Here is the summary of model :
Call:
lm(formula = y ~ x_1 +
Hi,
I think the simplest way to do this is using a small piece of php code.
See
http://stats.stackexchange.com/questions/4279/how-can-i-integrate-r-with-php
There are several webpages about this topic, but all of them are very
ancient. I have the plan to write a text about this
Greetz
Frans
www.
Hello everyone,
a problem with the plot.zoo function. In the parameters of the function,
cex.lab is ignored. I tried to reduce the size of the yaxis labels by at
least 50%.
--
Example:
sample <- as.zoo(EuStockMarkets)
par(las=1)
plot.zoo(sample, plot.type="multiple", main="Tim
You might prefer to read list via the web or a news reader or an RSS
feed, in which case http://dir.gmane.org/gmane.comp.lang.r.general might
be helpful.
Keith J
On 05/09/2012 13:47, David Winsemius wrote:
On Sep 5, 2012, at 3:00 AM, Marcus Tullius wrote:
Hello there,
is there a way I can
Hi R users
I'm looking for a way to do a contour plot of the attached 2D matrix. Using
contour() gives me some ugly contour lines that I would like to smooth. Any
ideas?
oliverjack38
http://r.789695.n4.nabble.com/file/n4642291/param1-4.txt param1-4.txt
--
View this message in context:
http
On 05.09.2012 11:41, Hermann Norpois wrote:
To whom it may concern.
As I would like to analyse some array data I was keen on downloading the
lumi package that depends obviously on hdrcde that is not available for r
2.12.1. I did not find instructions to solve or circumvent this problem.
Instal
On Fri, Aug 31, 2012 at 12:15 PM, David L Carlson wrote:
> Using a data.frame x with columns bins and counts:
>
> x <- structure(list(bins = c(3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5,
> 11.5, 12.5, 13.5, 14.5, 15.5), counts = c(1, 1, 2, 3, 6, 18,
> 19, 23, 8, 10, 6, 2, 1)), .Names = c("bi
On Wed, 5 Sep 2012, Nordlund, Dan (DSHS/RDA) wrote:
The line is drawn at 0.01 because there is a censored record where quant=0.01,
i.e.
sb.t[which(sb.t$quant==0.01 & sb.t$ceneq1==TRUE),]
site sampdateera param quant ceneq1 floor ceiling
34169 MW-5 1995-01-10 Pre-mining SbTot
On Sep 5, 2012, at 3:00 AM, Marcus Tullius wrote:
> Hello there,
>
> is there a way I can erase my name from the mailing list and still get mails
> concerning the topics
If you change your subscription options to daily digest ( at the page where you
signed up), you can just read the top of th
Hi
The key information is how landuse.raster and fire.did.not.occurr objects looks
like.
e.g. for matrix you can do
> set.seed(333)
> mat<-matrix(rnorm(20),5,4)
> mat[mat>0]<-NA
> fat<-matrix(sample(1:5, 20, replace=T),5,4)
> lat<-100
> mat[fat==1]<-lat
> mat
[,1][,2] [
Since noone but you knows what you are interested in - unlikely.
You could try to set up your mailclient with rules, deleting mails coming from
the list that do not contain certain keywords.
On 05.09.2012, at 12:00, Marcus Tullius wrote:
> Hello there,
>
> is there a way I can erase my name fr
Hello,
Em 05-09-2012 07:26, David Winsemius escreveu:
On Sep 4, 2012, at 4:39 PM, s.s.m. fauzi wrote:
Hi,
I have two table matrix, and I would like to compare the different between
two matrix.
For example:
Matrix 1:
A B C
A 0 1 0
B 0 0 1
C 0 0 0
Matrix 2:
A B C
A 0 1 0
B 0 0 0
agent dunham hotmail.com> writes:
> I've a tree which included at first 23 variables. Then I've pruned this
> tree, and there are only 8 variables involved.
Are you sure the others are not being considered as surrogate variables?
> Is there a way to introduce only this 8 values?
No. Why? See
Hello,
Jeff already answered to your question,
as.POSIXlt('2004-10-31 02:00:00', format = "%Y-%m-%d %H:%M:%S")
And in my system and timezone (GMT) I get
as.POSIXlt('2004-10-31 02:00:00')
[1] "2004-10-31 02:00:00"
So Jeff's point (b) should also be followed.
sessionInfo()
R version 2.15.1 (201
Hello there,
is there a way I can erase my name from the mailing list and still get mails
concerning the topics
I am interested in? I cannot handle all those emails coming in constantly. I
do not have time for that.
Greetings,
Francisco
[[alternative HTML version deleted]]
__
hi,
I am said to create a small form in html and after i click the submit
button it should execute R function and move to next page of html
can I know how can it be done???
can we link R and HTML???
which packages will help me in doing this???
are there any examples
To whom it may concern.
As I would like to analyse some array data I was keen on downloading the
lumi package that depends obviously on hdrcde that is not available for r
2.12.1. I did not find instructions to solve or circumvent this problem.
Installing hdrcde by hand did not work either. It was
On 09/05/2012 05:40 PM, Marcus Tullius wrote:
Hallo there,
greetings from Germany.
I have a simple question for you.
I have run a binary logistic model, but there are lots of outliers distorting
the real results.
I have tried to get rid of the outliers using the following commands:
Ank,
Welcome to R.
I suspect that what you are missing is:
require(betategarch)
If you don't know it, you may be interested in
the blog post:
http://www.portfolioprobe.com/2012/07/06/a-practical-introduction-to-garch-modeling/
Also questions about garch are generally best
sent to R-sig-finan
Hallo there,
greetings from Germany.
I have a simple question for you.
I have run a binary logistic model, but there are lots of outliers distorting
the real results.
I have tried to get rid of the outliers using the following commands:
remove = -c(56, 303, 365, 391, 512, 746, 859, 940,
I'll rephrase the question...
If you try...
as.POSIXlt('2004-10-31 02:00:00')
you get
[1] "2004-10-31"
What do I need to do to make it return
[1] "2004-10-31 02:00:00"
--
View this message in context:
http://r.789695.n4.nabble.com/POSIXlt-and-daylight-savings-time-tp4642253p4642272.ht
60 matches
Mail list logo
