similar to: density estimation

Dear R developpers, I am a user of Splus since many years and I have developped lots of functions to plot graph of data and model results of irregular or regular times series. In Splus regular times series are created using the rts function and irregular time series using the its functions. In both cases, times is given as a numeric vector. There is no problem in plotting a regular and an

Hi I have a netCDF file of volumetric soil water content at four different soil layers and want to convert each soil layer in the netCDF file to a GeoTIFF layer. This code converts the netCDF file to one GeoTIFF layer, i.e. unclear which soil depth. file.nc <- "C:/Soil_Weather_data/Agro/VMC21/VMC21.nc" file.tiff <- "C:/Soil_Weather_data/Agro /VMC21/VMC21" importnetcdf

Dear R-users, I have 1 remark and 1 question on the inclusion of interactions in the gam function from the gam package. I need to fit quantitative predictors in interactions with factors. You can see an example of what I need in fig 9.13 p265 from Hastie and Tibshirani book (1990). It's clearly stated that in ?gam "Interactions with nonparametric smooth terms are not fully

Hi Ahmed, When reading from a ncdf file you can use the 'varname', 'lvar' and 'level' arguments - see the 'Details' section in the docs https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/raster <https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/raster> We can't tell what is in the ncdf file from what you report other

Hello Guys I have problem with internet bank. I have 2 Internet links balancing mode, thus the bank is charging connection down. I tried to force Internet traffic (port 80 and 443) for only a link, however it did not work. How do I make a setting to force the connection to these ports for a specific link. Note: I can not use the file as route_rules have neither the source IP (ltsp) nor of

I am new to OpenSSH (I used SSH before). One thing that I miss from SSH is the "auto-complete". In the sftp from SSH you could complete path using <TAB> just like in Bash, but OpenSSH sftp does not do that. Is there any plan to add that feature ? If not, could it be add ? ;-) Thanks... -- Federico Petronio petronio at agro.uba.ar Linux User #129974 --- Unix

Dear all, Has anyone written R code for partial Mantel Tests- as described for instance in Legtendre & Legendre (1998) ? In other words, in a community ecology analysis, I would like to calculate the correlation between two dissimilarity matrices, controlling for a third distance matrix representing geographical distances between sites. Thanks! Christophe Bouget Biodiversité et gestion des

Hi all, i'm new to the list and here's my first question : I'm trying to install several printer drivers on my samba server with imprints. I did a package with the NT/2000 and W9x drivers and the mkprintpkg command. Then I used the command install_printer_drivers.pl to install the files in the appropriate directories. So now i'd like to check the drivers available and the

Dear R users, Up to now, I only used precompiled packages. As I am working on vegetation ecology, I would be interested in using a package not stored by CRAN called labdsv_0.9-1.tar.gz and developped by Dave Roberts at the National Center for Ecological Analysis and Synthesis (unfortunately I did not find any corresponding .zip file). Although I tried to follow the proposed guidelines to

Dear Isabelle, for most of my work I have switched from S-Plus to R, because it is easier to make scripts with due to the better documentation and some other more reasonable implementations. And there was also a cost issue with the S-Plus Server, which is not for free at all. In R I use the its package for time series, which is quite sufficient for my needs. The only problem I encountered was

Dear all, To modelize the abundance of fish (4 classes) with a set of environmental variables, I used the polr and predict.polr functions. I would like to know how to bring the cumulated probabilities back to a discrete ordinal scale. For the moment I used the predict.polr function with the argument "class". Is there an other way? polrf <- polrf <- polr_mod(formula =

Dear R-helpers, this mail is intended to mention a rather trange result and generate potential useful comments on it. I am not aware of another posts on this issue ( RSiteSearch("quasipoisson lmer version dispersion")). MUsing the exemple in the reference of the lmer function (in lme4 library) and turning it into a quasi-poisson or quasi-binomial analysis, we get different results,

Hi all, I have a very large binary vector, I wish to calculate the number of 1's over sliding windows. this is my very slow function slide<-function(seq,window){ n<-length(seq)-window tot<-c() tot[1]<-sum(seq[1:window]) for (i in 2:n) { tot[i]<- tot[i-1]-seq[i-1]+seq[i] } return(tot) } this works well for for reasonably sized vectors. Does

Hi all, Has anyone got a function to perform clustering by Flexible UPGMA (Un-weighted Pair Group Using Arithmetic Averaging) . It does not seem to exist in any R package. Regards, Ton -- Laboratoire d'Hydro?cologie Quantitative 3 bis, Quai Chauveau - CP 220 69336 LYON cedex 09 t?l. : +33 (0)4.72.20.89 20 fax : +33 (0)4.78.47.78.75 mail : ton.snelder at cemagref.fr

What you may need is to compile this package from source. It is explained in one of the R manuals (I can't remember which one though). You can also look at my "R Guide for Windows Users" at http://www.stat.auckland.ac.nz/~kwan022/rinfo.php (Section 5). I'm not sure the exact contents in the labdsv_0.9-1.tar.gz , but in general you can first unpack it with, e.g.: tar zxvf

Dear all, I have the following data frame: > raw.count Var1 Freq 1 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 707 2 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC 14 3 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAT 3 But why when printint it using 'cat', it doesn't print the desired string "AAA" ? > cat(raw.count$Var1,

Hello, >I am hoping for some advice as to how I might create a subset of a >matrix. The matrix is 176 x 3530. The rows are individual records >and the columns words. I want to create a new matrix that only >consists of words which occur in at least 5 records. For example, >if column 7 is "charges" and this only appears in 4 records/rows >this variable would not

Hello every one Working under Windows Xp in French, and R version 2.4.0 (and : > localeToCharset() : "ISO8859-1" ) I am trying to source a file with basic commands. source("R files/InputFile Analysis (Calculs et mises à l'échelle - gestion des pb d'unités).r") and it works fine as long as I do not use command like these ones :

Hello All, I have had considerable bad luck with attempting the following with for loops. Here is the problem: # Suppose we have a data.frame with the following data, which can be considered a type of edgelist (for those with networks backgrounds): # # V1 V2 # 1 A # 1 A # 1 B # 2 A # 3 C # 3 A # 3 C # 3 B # # I want the output of the function to produce a matrix, such that #each factor of

Hi all, Is any way to estimate the parameters of a curve, not manualy, from many subsets of my dataset [[alternative HTML version deleted]]