similar to: 64 Bit R Background Question

I have what R seems to consider a very large dataset, a 12MB text file of lat,long,and height values, 130,000 rows to be exact. Here's what I get: Thomas Colson North Carolina State University Department of Forestry and Environmental Resources (919) 673 8023 tom_colson at ncsu.edu Calendar: www4.ncsu.edu/~tpcolson

>From the following: basin.map <- readAsciiGrid("c:/temp/area.asc", colname="area") I have a SpatialGridDataFrame which has the x and y cordinate of a cell, and the drainage area of that cell. There are many cells with a low drainage area (in my case, 33000 with an area of 37.16) and one cell with the highest drainage area (again, in my case, a drainage area of of

Dear R users, I'm trying to write a small function in order to compute empirical cumulative density function.All seems to work but when I try to plot the function, I always get error messages. This is the function I use mecdf<-function(u,v,z) { u=sort(u) v=sort(v) n=length(u) nb=0 for (i in seq(1,n)) { if (u[i]<z & v[i]<z) { nb<-nb+1 } } nb=nb/n

I have got following error in drawing wavelet fitting. can some one help? > library(faraway) > data(lidar) > newlidar<-lidar[c(1:128),] > library(wavethresh) > wds <- wd(newlidar$logratio) > draw(wds) Error in plot.default(x = x, y = zwr, main = main, sub = sub, xlab = xlab, : formal argument "type" matched by multiple actual arguments [[alternative HTML

Hi All, I am a PhD student in forestry science and I am working with LiDAR data set (huge data set). I am a brand-new in R and geostatistic (SORRY, my background it?s in forestry) but I wish improve my skill for improve myself. I wish to develop a methodology to processing a large data-set of points (typical in LiDAR) but there is a problem with memory. I had created a subsample data-base but

I believe that actually means configure was able to use the BLAS you specified (MKL in this case). You can scan through config.log to make sure. Andy > From: > > Fedora Core 3 > > I installed the Intel MLk, > > and tried : --with-blas="-lmkl -lguide -lpthread" > > and got : External libraries: readline, BLAS(generic) > > thus I'm

Could someone give me a hand with the format of the gsummary function? Basically, I have a large set of xyz coordinates generated by LiDAR data (>37 million points) and I am trying to derive various summary statistics on the z-coordinates by a grid cell. I wrote a function to do this by creating factors from the x- and y- coordinates and then using gsummary. However, I want the function to

Hi I am afraid i am not understanding something very fundamental.... and does not matter how much i am looking into the book "Generalized Additive Models" of S. Wood i still don't understand my result. I am trying to model presence / absence (presence = 1, absence = 0) of a species using some lidar metrics (i have 4 of these). I am using different models and such .... and when i

I'm not sure I'm barking up the right tree here, but would I need to make use of groups to plot two separate datasets within ONE panel in xyplot? The desired end result is a single xy plot of two separate (but similar in values and ranges). Full code follows, xyplot code at bottom #########Determine Frequencies ##########coastal_slope #needs the maptools package to read ESRI grid

OK, I really dread to ask that .... much more that I know some discussion about p-values and if they are relevant for regressions were already on the list. I know to get p-val of regression coefficients - this is not a problem. But unfortunately one editor of a journal where i would like to publish some results insists in giving p-values for the squared deviance i get out from different glm and

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080716/338e44b9/attachment.pl>

Dear R-Help members, I am using Simon Wood`s mgcv package version1.7-22and R version 3.0.0 (2013-04-03) for fitting a GAM-Model to the LIDAR Data contained in the "SemiPar" package. Here is the code for fitting the model and for plotting the result: data("lidar") attach(lidar) ### # mgcv fitting ### gam_fit <- gam(logratio ~ s(range, k = 40, bs = "cr"), gamma

I'm using Samba 3.0.3 on Fedora Core 2 with OpenLDAP 2.1.29 for a backend. I'm getting to typical "The user name could not be found." error upon trying to join a Windows box. I've gone through every digest on lists.samba.org and other sites and nothing has worked yet. Any suggestions: Here's what I've done so far: 1. Installed everything via RPMS: [root@smbtest

After, oh, six months of attempts here and there to read everyone's experiences with Samba/LDAP and inability for a windows 2000/XP machine to join the domain, I finally discovered what was not working properly. In my smb.conf I put: add machine script = /usr/local/sbin/smbldap-useradd -w "%u" As instructed by many How-to's and Idealx. However, I thought to myself, %m

Folks, Mercurial is selectable as an installation in OpenSUSE 10.3 YAST. I have done the following: hg clone http://xenbits.xensource.com/xen-unstable.hg hg clone http://xenbits.xensource.com/linux-2.6.18-xen.hg make dist sh ./install.sh I tried running the following to create the initrd but it failed: Depending on your config, you may need to use ''mkinitrd'' to

Hello, I'm trying to estimate a Cox proportional hazard model with time-varying covariates using coxph. The parameter estimates are fine but there is something wrong with the survival curves I get with survfit (results are not plausible). Let me explain why I think something's wrong. To make sure I'm setting up my data correctly to estimate a model with time-varying covariates, I

I'm using R for PCA and? factor analysis. I want to create biplots of varimax rotated factors that color-code points by their classification. My research is on streams that are urban and rural. So, I want to color code them by this classification. If you just do a biplot from prcomp or princomp, you cannot add this color. So, I have used some code developed by a graduate student in our

I had yam working with ncsu.edu until the things changed when centos seperated the two projects. Now I can't get it to work. This is what I am trying to use but it says "nothing to do". updates = rsync://ftp.linux.ncsu.edu/centos/$release/$repo/$arch/RPMS/ This with gatech on the other hand works but they have a 20 connection limit so I would like to get ncsu.esu working again.

Dear R list, I have the following data: set.seed(1) n <- 5 # number of subjects tt <- 3 # number of repeated observation per subject numco <- 2 # number of covariates x <- matrix(round(rnorm(n*numco),2), ncol=numco) # the actual covariates x > x [,1] [,2] [1,] -0.63 -0.82 [2,] 0.18 0.49 [3,] -0.84 0.74 [4,] 1.60 0.58 [5,] 0.33 -0.31 I need to form a matrix

HI, there, I have a data set with special structure. It is in n*(5*p): n is the number of observations or data points 5*p is the matrix for each data point I'd like to conduct discriminant analysis to this data set. How could I do? And where could I find related references to solve this problem? Thanks a lot! Jieping Zhao PhD student in Bioinformatics, NCSU Lab