Displaying 20 results from an estimated 500 matches similar to: "R question"
2006 Jul 11
2
R newbie: logical subsets
Hello! I'm a newcomer to R hoping to replace some convoluted database
code with an R script. Unfortunately, I haven't been able to figure out
how to implement the following logic.
Essentially, we have a database of transactions that are coded with a
geographic locale and a type. These are being loaded into a data.frame
with named variables city, type, and price. E.g., trans$city
2005 Apr 29
0
Anscombe-Glynn, Bonett-Seier, D'Agostino
Dear useRs,
I was searching CRAN for implementation of kurtosis and skewness tests,
and found that there is some kind of lack on it.
So, I have written three functions:
1. Anscombe-Glynn test for kurtosis
2. Bonett-Seier test based on Geary's kurtosis (which is not widely
known, but I was inspired by original paper describing it, found
coincidentally in Elsevier database)
3.
2006 Sep 08
1
Computing skewness and kurtosis with the moments package
Hi,
I'm a newcomer to R, having previously used SPSS. One problem I have
run into is computing kurtosis. A test dataset is here:
http://www.whinlatter.ukfsn.org/2401.dat
> library(moments)
> data <- read.table("2401.dat", header=T)
> attach(data)
> loglen <- log10(Length)
With SPSS, I get
Skewness -0.320
Kurtosis -1.138
With R:
> skewness(loglen)
[1]
2008 May 29
1
plotting zoo using datetime as xlim
is there a way to use the actual index value for plotting zoo objects
this is the way that the index is set up and a sample range of what I would
like to plot
01/01/06 00:00:00 - 01/01/06 23:45:00
{
library(zoo)
# chron
library(chron)
fmt.chron <- function(x) {
chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x))
}}
x <- structure(c(15.57, 15.5,
2009 Oct 29
1
singular variance-covariance warning in lmer
Dear R Users,
I was hoping for some help with a recurrent error message in lmer. I am trying to model the effect of temperature on metabolic rate in animals (response = int.length) at different temperatures (mean.sst), with repeated measurements on the same individuals (random effect = female). Ideally I would make a random slope and intercept model where the rate can change differently with
2011 Dec 06
0
CentOS-announce Digest, Vol 82, Issue 3
Send CentOS-announce mailing list submissions to
centos-announce at centos.org
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.centos.org/mailman/listinfo/centos-announce
or, via email, send a message with subject or body 'help' to
centos-announce-request at centos.org
You can reach the person managing the list at
centos-announce-owner at centos.org
When
2011 May 06
1
read a netcdf file _Fill_value=-32768
Hello
I am a new user of R .
and I ve problem with R and netcdf .
I succed installation , I could use all examples .
But when I take my netcf it is different .
I want to do statistic on this kind of file .
1)
first calculate mean .
my data is like that
through ncdump -h test.nc
netcdf test {
dimensions:
lat = 301 ;
lon = 401 ;
time = UNLIMITED ; // (80 currently)
2009 Nov 20
1
problem selecting 330 to 30 longitude
Hello,
This is probably very simple but I am trying to select a region from the
OIv2 SST dataset that crosses 0 longitude. The longitude variable is
arranged from 0 to 360 and I get an error message when trying to go from 330
to 30 for example. The data are from a NOAA netCDF file.
*temp = open.ncdf("c:/documents and settings/theoni/desktop/sstoi_v2.nc")
print(temp)
y1 = get.var.ncdf(
2004 Oct 23
0
(PR#7309) misuse of R-bugs (was inappropriate definition
>From the R FAQ
Bug reports on contributed packages should be sent first to the package
maintainer, and only submitted to the R-bugs repository by package
maintainers, mentioning the package in the subject line.
and you are NOT the listed maintainer.
Do learn to read FAQs before causing unnecessary work for other people
(as the R posting guide asks).
Please do as the FAQ asks.
On
2004 Sep 21
2
Ever see a stata import problem like this?
Greetings Everybody:
I generated a 1.2MB dta file based on the general social survey with
Stata8 for linux. The file can be re-opened with Stata, but when I bring
it into R, it says all the values are missing for most of the variables.
This dataset is called "morgen.dta" and I dropped a copy online in case
you are interested
http://www.ku.edu/~pauljohn/R/morgen.dta
looks like this
2006 Nov 26
1
GLM and LM singularities
Hi-
I'm wrestling with some of my data apparently not being called into a GLM or
an LM. I'm looking at factors affecting fish annual catch rates (ie. CPUE)
over 30 years. Two of the factors I'm using are sea surface temperature and
sea surface temperature anomaly. A small sample of my data is below:
CPUE
Year
Vessel_ID
Base_Port
Boat_Lgth
Planing
SST
Anomaly
0.127
2011 Dec 12
2
i can't read large NETCDF file like CRU
i use library(ncdf) to read this file as follow
library(ncdf)
sst.nc = open.ncdf(title)
lonall = get.var.ncdf(sst.nc,'lon')
latall = get.var.ncdf(sst.nc,'lat')
precip = get.var.ncdf(sst.nc,'pre')
close(sst.nc)
if i use this method my pc freeze and not respond until i restart it is
there
--
View this message in context:
2012 Nov 16
2
R-Square in WLS
Hi,
I am fitting a weighted least square regression and trying to compute
SSE,SST and SSReg but I am not getting SST = SSReg + SSE and I dont know
what I am coding wrong. Can you help please?
xnam <-colnames(X) # colnames Design Matrix
fmla1 <- as.formula(paste("Y ~",paste(xnam, collapse=
2009 Nov 25
1
which to trust...princomp() or prcomp() or neither?
According to R help:
princomp() uses eigenvalues of covariance data.
prcomp() uses the SVD method.
yet when I run the (eg., USArrests) data example and compare with my own
"hand-written" versions of PCA I get what looks like the opposite.
Example:
comparing the variances I see:
Using prcomp(USArrests)
-------------------------------------
Standard deviations:
[1] 83.732400 14.212402
2013 Jan 16
1
setting the legend scale in contour plots
Hi,
I'd like to compare SST data for year 2000 with SST for year 2001. I managed
to get filled contour plots showing monthly SST for both years (12 plots for
each year, 24 plots in total). In order to compare year 2000 and year 2001,
however, I'd need to have the same legend scale (same color bar) in both
years plots.
Thus, my question is: is there a way to set the legend scale in
2005 Dec 01
2
about comparison of KURTOSIS in package: moments and fBasics
Hello
I do not know very much about statistics (and English language too :-( ),
then I come in search of a clarification (explanation):
I found two distinct results on KURTOSIS and
I do not know which of them is the correct one.
Any aid will be welcome!
klebyn
################ CODE
rnorm(1000) -> x
library(moments)
kurtosis(x)
skewness(x)
detach("package:moments")
2016 Apr 27
0
error.crosses
Hello all, I have used describeBy to generate the following summary
statistics. I simply need x and y error bars on a plot that has CQN
(xaxis) and Price (yaxis). There should be four total points on the graph
(one for each supplier).
Using "error.crosses(desc$CQN, desc$Price)" does not work.
group: a
vars n mean sd median trimmed mad min max range skew
CQN
2001 Oct 04
0
Summary on random data with zero skew and some kurtosis
Thanks to all who response my problem. Here are my summary :
1. from Dirk Eddelbuettel <edd at debian.org>
We could try a mixture of normals -- ie flip a coin (use a uniform with
some cutoff c where 0 < c < 1 ) to choose between N(0, sigma_1) and N(0,
sigma_2).
2. from Michaell Taylor <michaell.taylor at reis.com>
We could use the gld library to specify the lambdas of
2001 Oct 03
0
Summary : Generate random data from dist. with 0 skewness and some kurtosis
Thanks to all who response my problem. Here are my summary :
1. from Dirk Eddelbuettel <edd at debian.org>
We could try a mixture of normals -- ie flip a coin (use a uniform with
some cutoff c where 0 < c < 1 ) to choose between N(0, sigma_1) and N(0,
sigma_2).
2. from Michaell Taylor <michaell.taylor at reis.com>
We could use the gld library to specify the lambdas of
2008 Jun 24
1
Binding result of a function to a data frame
Hi,
I have the following function:
> kurtosis <-function(x) (mean((x-mean(x))^4))/(sd(x)^4) #x is a vector
and data
> print(mydata)
V1 V2 V3 V4 V5
1 1007_s_at DDR1 2865.1 2901.3 1978.3
2 1053_at RFC2 103.6 81.6 108.0
3