search for: piks

...on't know how serious my trouble is. After running a script: ** *t**<-c(14598417794,649693)* *data**=data.frame(read.spss("C:\\Ginters\\Kalibracija\\cal_data.sav"))* *Xs=**as.matrix(data[,1:2]) * *koef**=data.frame(read.spss("C:\\Ginters\\Kalibracija\\f.sav")) * *piks=**as.matrix(koef[,1**])* *g=regressionestimator(Xs,piks,t)* I get: *Error: cannot allocate vector of size 1614604 Kb In addition: Warning messages: 1: Reached total allocation of 255Mb: see help(memory.size) 2: Reached total allocation of 255Mb: see help(memory.size) * My OS is Win 2000 Proffe...

Hallo all just a small question I did not find an answer in help pages. Is it possible to use identify() after plotting with plot.POSIX to label points and/or to find out some points? Thanks a lot. Best regards Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz

I have a series of nertCDF files containing global data for a particular variable, e.g. tmin/tmax/precipiation/windspeed/relative humuidity/radiation etc. I get the following information when using *nc_open* function in R: datafile: https://www.dropbox.com/s/xpo7zklcmtm3g5r/gfdl_preci.nc?dl=0 File gfdl_preci.nc (NC_FORMAT_NETCDF4_CLASSIC): 1 variables (excluding dimension variables):

I would like to try a likelihood ratio test in place of waldtest. Ideally I'd like to provide two glm models, the second a submodel of the first, in the style of lrt (http://www.pik-potsdam.de/~hrust/tools/farismahelp/lrt.html). [lrt takes farimsa objects] Does anyone know of such a likelihood ratio test? Chris Elsaesser, PhD Principal Scientist, Machine Learning SPADAC Inc. 7921

Dear R experts, I want to plot a line chart with another secondary axis placed right to the standard secondary axis which one can access with the axis command, so that the data lines are seen in the same plot. Is there any way to do this in R? Many thanks, Kirsten.

Dear R users I searched some sources but i did not find an answer.Please give me some hint to following problem. I would like to compute a summary statistic for some vector for different factor levels. I know I can use tapply or aggregate but I do not know if there is a way how to use function with several (two) variable input (like weighted.mean). I wrote a simple a function for factor

Previously we''ve used Rails2.xx with no problems and we''re trying to upgrade to ruby 1.9.2p180 and Rails 3.0.7 with pik to keep ruby versions separate. After some initial difficulties with MySql2 gem installation, following <a href=http://rorguide.blogspot.com/2011/03/installing- mysql2-gem-on-ruby-192-and.html>these tips</a> got mysql2 v0.2.7 installed. My problem

Two questions: 1. Is the order of the dimensions shown what is shown if you look at str(ncin) - I mean shown at the end where it describes the variable and its dimensions? 2. Is you problem how to subset the netcdf file, how to write to the .csv file, or both? -Roy > On Aug 28, 2017, at 2:21 PM, Eeusha Nafi <eshad002 at gmail.com> wrote: > > I have a series of nertCDF files

Dear all I have a function in which I would like to change some values to NA according to some condition. dropout<-function(y, nahr=FALSE,...) { <some stuff for computing an index> if (nahr) y[index]<<-NA invisible(index) } in case y is a vector all works OK but if it is a part of data frame by calling dropout(df$y) or dropout(df[,number]) no change is done. Please can you

Hello. My target architecture supports natively 16 bit integers (i16). Whenever I write in C programs using only short types, clang compiles the program to LLVM and converts the i16 data to i32 to perform arithmetic operations and then truncates the results to i16. Then, the InstructionCombining (INSTCOMBINE or IC) pass removes these conversions back and forth from i16, except for

Hello, The question looks like simple. It's probably even stupid. But I spent several hours searching Internet, downloaded tons of papers, where deviance is mentioned and... And haven't found an answer. Well, it is clear for me the using of entropy when I split some node of a classification tree. The sense is clear, because entropy is an old good measure of how uniform is distribution.

Hi R users. I have a problem with function strata in sampling packages. > st0 = strata(dom, stratanames="stratas", size=sample.size, method="systematic",pik, FALSE) Error in sort.list(y) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list? In previous version of R 2.9.1 and previous version of package sampling this code worked

Dear all, I have to select 122 stratified random samples from a population of >3900 cells. I have 41 strata and I have to draw a different number of samples from them(between 2 and 8). I have tried to apply the funtion strata following the instruction in the manual: strata(dataframe, stratanames=NULL, size, method=c("srswor"), pik,description=TRUE) but I get the error

Dear R users I use R (1.3.0) Windows version usually with connection to Excel sheets via RODBC Version: 0.8-2. Everything is OK until I use mixed data types (numeric and character) in one column. In that case I will get a particular column of data frame with NA's. Example vzorek lpas apas NA 51.240 20.26 NA 51.340 20.42 NA 49.770 22.05 RB 12 49.840 19.30 RB 12 49.880

Dear all I try to persuade if statement in my function to work as I want (but I am not very successful:( My question is why my function with if statement is evaluated correctly in case of atomic variables and incorrectly in case of vector variables. Here is an example: #function with if statement fff <- function(otac,sklon) { test <- (otac *

Dear all I have a character vector of dates something like: timevec<-c("15.5.2003 00:00", "15.5.2003 00:01", "15.5.2003 00:02", "15.5.2003 00:03","15.5.2003 00:04") and I would like to transform it to some more convenient date class. Is there a way how to do it directly without previous reformating to ISO like structure and adding a

Dear all I would like to aggregate a data frame in a three minute interval. I have a time vector: > str(cas.param) `POSIXct', format: chr [1:181] "2003-06-12 09:00:00" "2003-06-12 09:01:00" "2003-06-12 09:02:00" and I want to turn it into a factor by cut. It works if I want to break it in let say hours > str(cut(cas.param,breaks="hour"))

Greetings, Does anyone know how to get an id number in the little header above each individual plot within a trellis plot? The default seems to be to print the word id and add a line indicating on a linear scale where the current id sits. Thanks in advance for any help you can send, Suzette ================================= Suzette Blanchard, Ph.D. Research Scientist Frontier Science

Hi, I want to connect to a share on a Samba server over the Internet. I have the Samba server installed and working OK (for NT clients at least) but when I connect from my NT workstation over the Internet I get 'Invalid Handle' dialogs when I try to map a network drive. I think it should be OK but one major difference is the remote NT workstation is not a member of the domain in which

Dear all Please, is there any function or package for dealing with multimodal distributions? I try to fit multimodal distribution or more precisely to find out mixture of normal distributions which can lead to my actual data. I use optim to find (in that case) two parameters but what I want is to let the function find out arbitrary number of normal distributions underlaying my actual data