similar to: Spatstat - K2 index

As far as I can tell from reading The Fine Documentation (R Language Definition and Intro to R), matrices are supposed to be of homogeneous types. Yet giving matrix() an inhomogeneous list seems to work, although it produces a peculiar object: v = list(1:3,4,5,"a") m = matrix(v,nrow=2) m [,1] [,2] [1,] Integer,3 5 [2,] 4 "a" m[1,] [[1]] [1] 1 2 3 [[2]]

Hello, I have fitted a Log-Gaussian Cox Process on my data but when I try to use "simulate.kppm" of the spatstat package I get this error: Error in rLGCP(model = model, mu = mu, param = param, ..., win = win) : The spatial domain of the pixel image ?mu? does not cover the simulation window ?win? I've used covariates as im, changed the npixel value of the spatstat options to

Hello all, I'm attempting to conduct spatial analysis of trees within a plot. I want to see if the trees are spatially correlated to soil characteristics, say pH, or moisture content. I think one way to do it is with mpl, however, my soil characteristics were not taken at exactly the same locations as my trees and further, the vectors aren't the same length. I'm getting the

Hi Everyone, I am trying to figure out the spatstat package for the first time and am having some trouble. Unfortunately, I can't post my data set but I'll hopefully post enough details for some help. I want to model the intensity of a spatial point process using 2 covariates from my data. After reading through the documentation, I have successfully created 2 "ppp" objects. The

Hi, I started using R about 3 months ago so please excuse my ignorance. I have two datasets. The first consists of 247 fixed position mosquito traps that were serviced weekly for 23 weeks. The second is dengue incidence data indexed to home address for the same time period. I would like to see if there is a space time relationship between mosquito catch rate and dengue incidence. So far I have

Hi, Is there an implementation of the D'agostino K2 omnibus test in an R package. I've found the D test in the fBasics package, but this doesn't give the same answer as my algorithm. I've already coded a version of the K2 and need to check the results. I have a set of spectra with gaussian distributed noise and I want a statistical test to separate those with no signal (so the

Dear all I want to simulate a stochastic jump variance process where N is Bernoulli with intensity lambda0 + lambda1*Vt. lambda0 is constant and lambda1 can be interpreted as a regression coefficient on the current variance level Vt. J is a scaling factor How can I rewrite this avoiding the loop structure which is very time-consuming for long simulations? for (i in 1:N){ ... N <- rbinom(n=1,

Hi R Team!, ? We are opening a new start up soon, codinghoodies.com, we want to make coders feel stylish. ? Out of goodwill I wanted to ask you formally if I can have permission to use the standard R logo on the front of hoodies to sell? I have read your rules but wanted to ask as I feel a moral right to email you asking to show support and respect for the R project. ? If it makes it easier I

I have a problem with file locking on different platforms. System environment: Samba Server 2.2.3.a on Red Hat 7.3 Some Windows 2000 Clients with MS-Office 2000 and OpenOffice 1.0.1 Some Red Hat 8.0 Clients with Open Office 1.0.1 The Linux Clients mounts the Samba Shares with smbmount: smbmount //<Server/shares> $HOME/path -o

My guess: It I clear from the link that they can use the R logo for commercial purposes. The issue is what to do about the "appropriate credit" and "link to the license." How would I do that on a hoodie? Would they need a web address or something? -----Original Message----- From: R-help <r-help-bounces at r-project.org> On Behalf Of John Fox Sent: Tuesday, March 21, 2023

I have the following problem. It is not of earthshaking importance, but still I have spent a considerable amount of time thinking about it. PROBLEM: Is there any way I can have a single textfile that contains both a) data b) programcode The program should act on the data, if the textfile is source()'ed into R. BOUNDARY CONDITION: I want the data written in the textfile in exactly

Hi, I have collected data on trees from 5 forest plots located within the same landscape. Data within the plots are spatially autocorrelated (calculated using Moran's I). I would like to do a ANCOVA type of analysis combining these five plots, but the assumption that there is no autocorrelation in the residuals is obviously violated. Does anyone have any ideas how to incorporate these spatial

Hola a tod en s, estoy realizando unas operaciones con muchos datos con el paquete "vegan", concretamente la función "metaMDS" y me gustaría saber si puedo usar más de un procesador a la vez para acelerar los cálculos. He visto el paquete "multicore" pero no sé usarlo, ¿alguien lo ha usado para algo así? Gracias y saludos, Javier -- http://oikos.inf.um.es

Rainer M Krug writes: > Thorsten Wiegand used in his paper Wiegand T., and K. A. Moloney 2004. > Rings, circles and null-models for point pattern analysis in ecology. > Oikos 104: 209-229 a statistic he called O-Ring statistic which is > similar to Ripley's K, only that it uses rings instead of circles. > > http://www.oesa.ufz.de/towi/towi_programita.html#ring

Dear list, following a previous querry we are still stuck! As pointed out by Erin Hodges the "ts" library includes a PACF function which reports the partial correlation of population density at time t against lagged population density. However, what we are trying to calculate is the partial correlation between rate of population change, Rt=log Nt/Nt-1, against lagged population

Dear all, A question related to the following has been asked on R-help before, but I could not find any answer to it. Input will be much appreciated. I got an unexpected sign of the "slope" parameter associated with a covariate (diam) using zeroinfl(). It led me to compare the estimates given by zeroinfl() and hurdle(): The (significant) negative estimate here is surprising, given

Hello, for my meta-analysis I try to test if two varainces are equal without using the raw scores. I have is the SD's, N's and the Means. I want to test the variances from dependent and independend samples. I assume I can use the var.test procedure for the independent samples, but what about the dependent samples ? Has anyone an idea how to realise this with R ? Thanks in advance

Hi Everyone, I am doing a project based on "Spatially Lagged Predictor Variable Models", I would like to know which package in R would execute this model. Also, I am new to this field of spatial statistics. Any suggestions for a good book on spatial regression analysis would be appreciated. Thanks Again. Cheers Arun -- View this message in context:

Hi I?m trying to fit a nonlinear model to a derivative of the logistic function y = a/(1+exp((b-x)/c)) (this is the parametrization for the SSlogis function with nls) The derivative calculated with D function is: > logis<- expression(a/(1+exp((b-x)/c))) > D(logis, "x") a * (exp((b - x)/c) * (1/c))/(1 + exp((b - x)/c))^2 So I enter this expression in the nls function:

Why do I get different sign of the coefficients of covariates when I run the semi-parametric proportional hazard model (coxph) compared to the parametric proportional hazard model (survreg)? Anyone with experience in extracting information form survreg to make predictions are free to contact me. Cheers, Ullrika [[alternative HTML version deleted]]