search for: pik

...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 Proff...

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

...Model output climate of GFDL-ESM2M r1i1p1 Interpolated to 0.5 degree and bias corrected using observations from 1960 - 1999 for EU WATCH project CDO: Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo) product_id: input model_id: gfdl-esm2m institute_id: PIK experiment_id: historical ensemble_id: r1i1p1 time_frequency: daily creator: isimip at pik-potsdam.de description: GFDL-ESM2M bias corrected impact model input prepared for ISIMIP2 Now I want to extract a subset from this dataset using pair of lon and lat po...

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 Jones Branch Dr. Suite 600 McLean, VA 22102 703.371.7301 (m) 703.637.9421 (o)

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.

...ion for factor weighted mean fff<-function(x,fact,w) { ws<-tapply(w,fact,sum) newx<-x*w tapply(newx,fact,sum)/ws } which can handle particular case but does exist some more general solution how to use FUN(X1,X2) in aggregation procedures (tapply, aggregate, by) directly? Thank you Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body",...

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 is that mysql2 hangs whenever I invok...

...FDL-ESM2M r1i1p1 Interpolated > to 0.5 degree and bias corrected using observations from 1960 - 1999 > for EU WATCH project > CDO: Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo) > product_id: input > model_id: gfdl-esm2m > institute_id: PIK > experiment_id: historical > ensemble_id: r1i1p1 > time_frequency: daily > creator: isimip at pik-potsdam.de > description: GFDL-ESM2M bias corrected impact model input > prepared for ISIMIP2 > > Now I want to extract a subset from this d...

...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 help me what is wrong with my code? By the way idx<-dropout(df$y) df$y[idx]<-NA works OK Thanks a lot beforehand Best regards. Petr Pikal Precheza a.s., Nab?.Dr.E.Bene?e 24, 750 62 P?erov tel: +420581 252 257 ; 724 008 364 petr.pikal at precheza.cz; p.pik at volny.cz fax +420581 252 561

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

...rm is distribution. And we want, for sure, the distribution to be uniform, represent one class only as the best. Where deviance come from at all? I look at a formula and see that the only difference to entropy is use of *number* of each class points, instead of *probability* as a multiplier of log(Pik). So, it looks like the deviance and entropy differ by factor 1/N (or 2/N), where N is total number of cases. Then WHY to say "deviance"? Any historical reason? Or most likely I do not understand something very basic. Please, help. Thanks, Alexander Skomorokhov, -------------- next part...

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 well and now I don't know what is a problem? Any ideas how to solve this problem...

...amples 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 message Error in data.frame((1:nrow(data))[srswor(size, nrow(data)) == 1], rep(size/nrow(data), : arguments imply differing number of rows: 3, 122 My stratanames and size vectors have the same length, 41 (the number of strata), and 122 is the number of s...

...any clue other then I(....) but I can not resolve how to use it in that case. Is there any possibility how to handle different data types in one Excel sheet column directly by using ODBC connection and keep all values in such column as they are when moving data from Excel to R? Thank you. Petr Pikal Precheza a.s. Nabr.Dr.E.Benese 24 tel: 00420 641 25 2257 +420 (0)724 008 364 petr.pikal at precheza.cz p.pik at volny.cz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send...

...* 3.69683+286.6*(1/sklon)-0.03100345*otac*sklon) ifelse(test<65, (otac/sklon * 141.76), (otac * 3.69683+286.6*(1/sklon)- 0.03100345*otac*sklon)) } > fff(seq(20,50,5),70) [1] 40.50286 50.62857 60.75429 70.88000 65.15783 72.79077 80.42371 #correct Can anybody help? Thank you in advance Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body",...

...out previous reformating to ISO like structure and adding a seconds to it. I tried direct transformation but it does not work as I expected (or as I understood from help pages :-). > strptime(timevec,"%m.%d.%y %H:%M") [1] NA NA NA NA NA Thank you for any hint. Best regards. Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz

...; "2003-06-12 09:03:00" "2003-06-12 09:06:00" works. So it is definitely my lack of understanding how to implement a vector of cut points (constructed by seq()), to break parameter of cut() function. Please can you give me any hints? R 1.7.0, Windows NT Thanks a lot Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz

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

...bit more general solution fff<-function(p) { p1<-p[1] p2<-p[2] v<-dnorm(x,30,p1)+dnorm(x,60,p2) s<-sum(v) sum((y-v/s)^2) } fff2<-function(p) { p1<-p[1] p2<-p[2] v<-dnorm(x,30,p1)+dnorm(x,60,p2) s<-sum(v) v/s } Thank you very much for any help. Best regards. Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz