Efficient Optimization Algorithms for Nonlinear Data Analysis

Identification of low-dimensional structures and main sources of variation from multivariate data are fundamental tasks in data analysis. Many methods aimed at these tasks involve solution of an optimization problem. Thus, the objective of this thesis is to develop computationally efficient and theoretically justified methods for solving such problems. Most of the thesis is based on a statistical model, where ridges of the density estimated from the data are considered as relevant features. Finding ridges, that are generalized maxima, necessitates development of advanced optimization methods. An efficient and convergent trust region Newton method for projecting a point onto a ridge of the underlying density is developed for this purpose. The method is utilized in a differential equation-based approach for tracing ridges and computing projection coordinates along them. The density estimation is done nonparametrically by using Gaussian kernels. This allows application of ridge-based methods with only mild assumptions on the underlying structure of the data. The statistical model and the ridge finding methods are adapted to two different applications. The first one is extraction of curvilinear structures from noisy data mixed with background clutter. The second one is a novel nonlinear generalization of principal component analysis (PCA) and its extension to time series data. The methods have a wide range of potential applications, where most of the earlier approaches are inadequate. Examples include identification of faults from seismic data and identification of filaments from cosmological data. Applicability of the nonlinear PCA to climate analysis and reconstruction of periodic patterns from noisy time series data are also demonstrated. Other contributions of the thesis include development of an efficient semidefinite optimization method for embedding graphs into the Euclidean space. The method produces structure-preserving embeddings that maximize interpoint distances. It is primarily developed for dimensionality reduction, but has also potential applications in graph theory and various areas of physics, chemistry and engineering. Asymptotic behaviour of ridges and maxima of Gaussian kernel densities is also investigated when the kernel bandwidth approaches infinity. The results are applied to the nonlinear PCA and to finding significant maxima of such densities, which is a typical problem in visual object tracking.

Data analysis tools and methods for DNA microarray and high-throughput sequencing data

The recent rapid development of biotechnological approaches has enabled the production of large whole genome level biological data sets. In order to handle thesedata sets, reliable and efficient automated tools and methods for data processingand result interpretation are required. Bioinformatics, as the field of studying andprocessing biological data, tries to answer this need by combining methods and approaches across computer science, statistics, mathematics and engineering to studyand process biological data. The need is also increasing for tools that can be used by the biological researchers themselves who may not have a strong statistical or computational background, which requires creating tools and pipelines with intuitive user interfaces, robust analysis workflows and strong emphasis on result reportingand visualization. Within this thesis, several data analysis tools and methods have been developed for analyzing high-throughput biological data sets. These approaches, coveringseveral aspects of high-throughput data analysis, are specifically aimed for gene expression and genotyping data although in principle they are suitable for analyzing other data types as well. Coherent handling of the data across the various data analysis steps is highly important in order to ensure robust and reliable results. Thus,robust data analysis workflows are also described, putting the developed tools andmethods into a wider context. The choice of the correct analysis method may also depend on the properties of the specific data setandthereforeguidelinesforchoosing an optimal method are given. The data analysis tools, methods and workflows developed within this thesis have been applied to several research studies, of which two representative examplesare included in the thesis. The first study focuses on spermatogenesis in murinetestis and the second one examines cell lineage specification in mouse embryonicstem cells.

Optimisation of data analysis method for a fluidized bed combustion test rig

Nowadays the used fuel variety in power boilers is widening and new boiler constructions and running models have to be developed. This research and development is done in small pilot plants where more faster analyse about the boiler mass and heat balance is needed to be able to find and do the right decisions already during the test run. The barrier on determining boiler balance during test runs is the long process of chemical analyses of collected input and outputmatter samples. The present work is concentrating on finding a way to determinethe boiler balance without chemical analyses and optimise the test rig to get the best possible accuracy for heat and mass balance of the boiler. The purpose of this work was to create an automatic boiler balance calculation method for 4 MW CFB/BFB pilot boiler of Kvaerner Pulping Oy located in Messukylä in Tampere. The calculation was created in the data management computer of pilot plants automation system. The calculation is made in Microsoft Excel environment, which gives a good base and functions for handling large databases and calculations without any delicate programming. The automation system in pilot plant was reconstructed und updated by Metso Automation Oy during year 2001 and the new system MetsoDNA has good data management properties, which is necessary for big calculations as boiler balance calculation. Two possible methods for calculating boiler balance during test run were found. Either the fuel flow is determined, which is usedto calculate the boiler's mass balance, or the unburned carbon loss is estimated and the mass balance of the boiler is calculated on the basis of boiler's heat balance. Both of the methods have their own weaknesses, so they were constructed parallel in the calculation and the decision of the used method was left to user. User also needs to define the used fuels and some solid mass flowsthat aren't measured automatically by the automation system. With sensitivity analysis was found that the most essential values for accurate boiler balance determination are flue gas oxygen content, the boiler's measured heat output and lower heating value of the fuel. The theoretical part of this work concentrates in the error management of these measurements and analyses and on measurement accuracy and boiler balance calculation in theory. The empirical part of this work concentrates on the creation of the balance calculation for the boiler in issue and on describing the work environment.

Spurring innovation through idea contests : data analysis and process improvements

This research concerns the Urban Living Idea Contest conducted by Creator Space™ of BASF SE during its 150th anniversary in 2015. The main objectives of the thesis are to provide a comprehensive analysis of the Urban Living Idea Contest (ULIC) and propose a number of improvement suggestions for future years. More than 4,000 data points were collected and analyzed to investigate the functionality of different elements of the contest. Furthermore, a set of improvement suggestions were proposed to BASF SE. Novelty of this thesis lies in the data collection and the original analysis of the contest, which identified its critical elements, as well as the areas that could be improved. The author of this research was a member of the organizing team and involved in the decision making process from the beginning until the end of the ULIC.

Tilastollisia perusmenetelmiä yleistajuisesti

Kirjallisuusarvostelu

Configuring Health Care Services - Service Control Improvement

The costs of health care are increasing, and at the same time, population is aging. This leads health care organizations to focus more on home based care services. This thesis focuses on the home care organization of the South Karelian District of Social and Health Services (Eksote), which was established in 2010; how its operation is organized and managed, and which problem types are faced in the daily operation of home care. This thesis examines home care services through an extensive interview study, process mapping and statistical data analysis. To be able to understand the nature of services and special environment theory models, such as service management and performance measurement, service processes and service design are introduced. This study is conducted from an external researcher‟s point of view and should be used as a discussion opener. The outcome of this thesis is an upper level development path for Eksote home care. The organization should evaluate and build a service offering, then productize home care services and modularize the products and identify similarities. Service processes should be mapped to generate efficiency for repeating tasks. Units should be reasonably sized and geographically located to facilitate management and operation. All this can be done by recognizing the different types of service products: runners repeaters and strangers. Furthermore, the organization should not hide behind medical issues and should understand the legislative, medical and operational frameworks in health care.

Analysis of patterns in mining time series using Qlucore

Raw measurement data does not always immediately convey useful information, but applying mathematical statistical analysis tools into measurement data can improve the situation. Data analysis can offer benefits like acquiring meaningful insight from the dataset, basing critical decisions on the findings, and ruling out human bias through proper statistical treatment. In this thesis we analyze data from an industrial mineral processing plant with the aim of studying the possibility of forecasting the quality of the final product, given by one variable, with a model based on the other variables. For the study mathematical tools like Qlucore Omics Explorer (QOE) and Sparse Bayesian regression (SB) are used. Later on, linear regression is used to build a model based on a subset of variables that seem to have most significant weights in the SB model. The results obtained from QOE show that the variable representing the desired final product does not correlate with other variables. For SB and linear regression, the results show that both SB and linear regression models built on 1-day averaged data seriously underestimate the variance of true data, whereas the two models built on 1-month averaged data are reliable and able to explain a larger proportion of variability in the available data, making them suitable for prediction purposes. However, it is concluded that no single model can fit well the whole available dataset and therefore, it is proposed for future work to make piecewise non linear regression models if the same available dataset is used, or the plant to provide another dataset that should be collected in a more systematic fashion than the present data for further analysis.

Nursing students and older people nursing. Towards a future career

The purpose of this two-phased study is to examine the interest of nursing students in choosing a career in older people nursing. First, the scoping phase explores the different premises for choosing older people nursing as a career. Second, the evaluation phase investigates the outcomes of the developed educational intervention involving older people as promoters of choosing a career in older people nursing, factors related to these outcomes, and experiences with educational intervention. The ultimate goal is to encourage more nursing students to choose older people nursing as their career. The scoping phase applies an exploratory design and centres around a descriptive, cross-sectional survey, documentary research and a scoping literature review. The information sources for this phase include 183 nursing students, 101 newspaper articles and 66 research articles. The evaluation phase applies a quasi-experimental design and a pre-post-test design with a non-equivalent comparison group and a post-intervention survey. The information sources for this phase include 87 nursing students and 43 older people. In both phases, statistical and narrative methods are applied in the data analysis. Nursing students neutrally regarded the idea of a career in older people nursing. The most consistent factors related to the nursing students’ career plans in older people nursing were found to be nursing work experience and various educational preparations in the field. Nursing students in the intervention group (n=40) were more interested in older people nursing and had more positive attitudes towards older people than did students in the comparison group (n=36). However, in both groups, the interest that students had at the baseline was associated with the interest at the one-month follow-up. There were no significant differences between the groups in terms of the students’ knowledge levels about ageing. The nursing students and older people alike highly appreciated participating in the educational intervention. It seems possible to positively impact nursing students and their choices to pursue careers in older people nursing, at least in the short-term. The involvement of older people as promoters of this career choice provides one encouraging alternative for impacting students’ career choices, but additional research is needed.

Monitoring of crystallization processes by using infrared spectroscopy and multivariate methods

Väitöstutkimuksessa on tarkasteltuinfrapunaspektroskopian ja monimuuttujaisten aineistonkäsittelymenetelmien soveltamista kiteytysprosessin monitoroinnissa ja kidemäisen tuotteen analysoinnissa. Parhaillaan kiteytysprosessitutkimuksessa maailmanlaajuisesti tutkitaan intensiivisesti erilaisten mittausmenetelmien soveltamista kiteytysprosessin ilmiöidenjatkuvaan mittaamiseen niin nestefaasista kuin syntyvistä kiteistäkin. Lisäksi tuotteen karakterisointi on välttämätöntä tuotteen laadun varmistamiseksi. Erityisesti lääkeaineiden valmistuksessa kiinnostusta tämäntyyppiseen tutkimukseen edistää Yhdysvaltain elintarvike- ja lääkeaineviraston (FDA) prosessianalyyttisiintekniikoihin (PAT) liittyvä ohjeistus, jossa määritellään laajasti vaatimukset lääkeaineiden valmistuksessa ja tuotteen karakterisoinnissa tarvittaville mittauksille turvallisten valmistusprosessien takaamiseksi. Jäähdytyskiteytyson erityisesti lääketeollisuudessa paljon käytetty erotusmenetelmä kiinteän raakatuotteen puhdistuksessa. Menetelmässä puhdistettava kiinteä raaka-aine liuotetaan sopivaan liuottimeen suhteellisen korkeassa lämpötilassa. Puhdistettavan aineen liukoisuus käytettävään liuottimeen laskee lämpötilan laskiessa, joten systeemiä jäähdytettäessä liuenneen aineen konsentraatio prosessissa ylittää liukoisuuskonsentraation. Tällaiseen ylikylläiseen systeemiin pyrkii muodostumaan uusia kiteitä tai olemassa olevat kiteet kasvavat. Ylikylläisyys on yksi tärkeimmistä kidetuotteen laatuun vaikuttavista tekijöistä. Jäähdytyskiteytyksessä syntyvän tuotteen ominaisuuksiin voidaan vaikuttaa mm. liuottimen valinnalla, jäähdytyprofiililla ja sekoituksella. Lisäksi kiteytysprosessin käynnistymisvaihe eli ensimmäisten kiteiden muodostumishetki vaikuttaa tuotteen ominaisuuksiin. Kidemäisen tuotteen laatu määritellään kiteiden keskimääräisen koon, koko- ja muotojakaumansekä puhtauden perusteella. Lääketeollisuudessa on usein vaatimuksena, että tuote edustaa tiettyä polymorfimuotoa, mikä tarkoittaa molekyylien kykyä järjestäytyä kidehilassa usealla eri tavalla. Edellä mainitut ominaisuudet vaikuttavat tuotteen jatkokäsiteltävyyteen, kuten mm. suodattuvuuteen, jauhautuvuuteen ja tabletoitavuuteen. Lisäksi polymorfiamuodolla on vaikutusta moniin tuotteen käytettävyysominaisuuksiin, kuten esim. lääkeaineen liukenemisnopeuteen elimistössä. Väitöstyössä on tutkittu sulfatiatsolin jäähdytyskiteytystä käyttäen useita eri liuotinseoksia ja jäähdytysprofiileja sekä tarkasteltu näiden tekijöiden vaikutustatuotteen laatuominaisuuksiin. Infrapunaspektroskopia on laajalti kemian alan tutkimuksissa sovellettava menetelmä. Siinä mitataan tutkittavan näytteenmolekyylien värähtelyjen aiheuttamia spektrimuutoksia IR alueella. Tutkimuksessa prosessinaikaiset mittaukset toteutettiin in-situ reaktoriin sijoitettavalla uppoanturilla käyttäen vaimennettuun kokonaisheijastukseen (ATR) perustuvaa Fourier muunnettua infrapuna (FTIR) spektroskopiaa. Jauhemaiset näytteet mitattiin off-line diffuusioheijastukseen (DRIFT) perustuvalla FTIR spektroskopialla. Monimuuttujamenetelmillä (kemometria) voidaan useita satoja, jopa tuhansia muuttujia käsittävä spektridata jalostaa kvalitatiiviseksi (laadulliseksi) tai kvantitatiiviseksi (määrälliseksi) prosessia kuvaavaksi informaatioksi. Väitöstyössä tarkasteltiin laajasti erilaisten monimuuttujamenetelmien soveltamista mahdollisimman monipuolisen prosessia kuvaavan informaation saamiseksi mitatusta spektriaineistosta. Väitöstyön tuloksena on ehdotettu kalibrointirutiini liuenneen aineen konsentraation ja edelleen ylikylläisyystason mittaamiseksi kiteytysprosessin aikana. Kalibrointirutiinin kehittämiseen kuuluivat aineiston hyvyyden tarkastelumenetelmät, aineiston esikäsittelymenetelmät, varsinainen kalibrointimallinnus sekä mallin validointi. Näin saadaan reaaliaikaista informaatiota kiteytysprosessin ajavasta voimasta, mikä edelleen parantaa kyseisen prosessin tuntemusta ja hallittavuutta. Ylikylläisyystason vaikutuksia syntyvän kidetuotteen laatuun seurattiin usein kiteytyskokein. Työssä on esitetty myös monimuuttujaiseen tilastolliseen prosessinseurantaan perustuva menetelmä, jolla voidaan ennustaa spontaania primääristä ytimenmuodostumishetkeä mitatusta spektriaineistosta sekä mahdollisesti päätellä ydintymisessä syntyvä polymorfimuoto. Ehdotettua menetelmää hyödyntäen voidaan paitsi ennakoida kideytimien muodostumista myös havaita mahdolliset häiriötilanteet kiteytysprosessin alkuhetkillä. Syntyvää polymorfimuotoa ennustamalla voidaan havaita ei-toivotun polymorfin ydintyminen,ja mahdollisesti muuttaa kiteytyksen ohjausta halutun polymorfimuodon saavuttamiseksi. Monimuuttujamenetelmiä sovellettiin myös kiteytyspanosten välisen vaihtelun määrittämiseen mitatusta spektriaineistosta. Tämäntyyppisestä analyysistä saatua informaatiota voidaan hyödyntää kiteytysprosessien suunnittelussa ja optimoinnissa. Väitöstyössä testattiin IR spektroskopian ja erilaisten monimuuttujamenetelmien soveltuvuutta kidetuotteen polymorfikoostumuksen nopeaan määritykseen. Jauhemaisten näytteiden luokittelu eri polymorfeja sisältäviin näytteisiin voitiin tehdä käyttäen tarkoitukseen soveltuvia monimuuttujaisia luokittelumenetelmiä. Tämä tarjoaa nopean menetelmän jauhemaisen näytteen polymorfikoostumuksen karkeaan arviointiin, eli siihen mitä yksittäistä polymorfia kyseinen näyte pääasiassa sisältää. Varsinainen kvantitatiivinen analyysi, eli sen selvittäminen paljonko esim. painoprosentteina näyte sisältää eri polymorfeja, vaatii kaikki polymorfit kattavan fysikaalisen kalibrointisarjan, mikä voi olla puhtaiden polymorfien huonon saatavuuden takia hankalaa.

Perusterveydenhuollon ja erikoissairaanhoidon käytön seuranta ja ennakointi kuntanäkökulmasta

Follow-up of utilisation and prediction of primary health care and hospital care from the municipality point of view. Planning, follow-up, and evaluation of primary health care within municipality entail comprehensive information about factors that influence health. In addition to populationbased research, various statistical data and registries serve as sources of information. The present study examined utilisation of primary health care and hospital care with the existing databases, registries, and categorization of Diagnosis Related Groups (DRGs) from the municipality (purchaser) point of view. Research involving the cases of Paimio, Sauvo, and Turku as examples of municipalities pointed out that, even in the small municipalities, it is possible to assess and predict health services to be offered to the inhabitants by following databases and registries. Health-related databases and registries include a plenty of possible uses that have not adequately been employed at the level of municipality. Descriptive futures research and community analysis formed the framework of the study. Descriptive futures research may be used to establish predictions based on past developmental traditions, and quantitative time trend analyses may be employed to make estimations about future events. Community analysis will assist in making conclusions about population- based health care needs, in assessing the functionality or effectiveness of the health care system, and in appropriately targeting limited resources. The aim of the present study was to describe the health service profile so that the arrangements and planning of health services as well as the contract negotiations of hospital care become easier within municipalities. Another aim was to assess the application of Hilmo (registry for posting hospital care periods), Aitta and Sotka (statistical databases) for the purposes of resource planning in the procurement of hospital care. A third aim was to evaluate how the system of the DRGs adapts in the prediction of retaining health services within short (1-year), intermediate (5-year) and long range (10-15-year) intervals. The findings indicated that the follow-up of primary health care utilisation combined with follow-up of hospital care utilisation allows municipalities to plan and predict health services when databases are applied. Information about the past contacts with the databases has indicated that the health care culture and incidence of disease change rather slowly in the area of investigation. For the purposes of health care research, it is recommended that methods of application used in making predictions about health care utilisation need to be further developed

Puunjalostuskonsernin voiteluhuollon analysoiminen ja kehittäminen

Diplomityön tarkoituksena on analysoida ja kehittää erään puunjalostuskonsernin voiteluhuoltoa. Tutkimuksen kirjallisessa osassa on selvitetty voitelun perusteet, voiteluaineet sekä voiteluhuoltoon liittyviä yleisiä käsitteitä ja termejä. Työssä on selvitetty ohjeita voiteluaineiden valintaan ja voiteluöljyjen kunnonvalvontamenetelmiin sekä selvitetty voiteluaineiden puhtauden merkitystä. Lisäksi on esitetty katselmus eri voitelumekanismista. Näiden sovellusesimerkkinä on laskettu paperikoneen kuivaussylinterin laakerin voitelukalvon ominaispaksuus. Osan lopussa on tutkittu voiteluaineiden varastoinnin ja käytettyjen voiteluaineiden jatkokäsittelyn merkitystä ympäristöasiat huomioiden. Soveltavassa osassa on edellistä hyödyntäen analysoitu puunjalostuskonsernin yksitoista eri paikkakuntaa voiteluhuollon toimivuuden kannalta. Tämä tilastollinen tieto on kerätty sekä kyselykaavakkeen avulla että haastattelemalla jokaisen paikkakunnan voiteluhuollon asiantuntijoita. Lopputuloksena on saatu voiteluhuollon kehityskohteet ja SWOT-analyysit paikkakunnittain sekä yhteenveto voiteluhuollon eri toimintojen parhaimmista käytännöistä. Lisäksi analyysin perusteella laadittiin konsernin voiteluhuollon kehityssuunnitelma.

Advanced chemometric methods: Applicability on industrial data

Recent years have produced great advances in the instrumentation technology. The amount of available data has been increasing due to the simplicity, speed and accuracy of current spectroscopic instruments. Most of these data are, however, meaningless without a proper analysis. This has been one of the reasons for the overgrowing success of multivariate handling of such data. Industrial data is commonly not designed data; in other words, there is no exact experimental design, but rather the data have been collected as a routine procedure during an industrial process. This makes certain demands on the multivariate modeling, as the selection of samples and variables can have an enormous effect. Common approaches in the modeling of industrial data are PCA (principal component analysis) and PLS (projection to latent structures or partial least squares) but there are also other methods that should be considered. The more advanced methods include multi block modeling and nonlinear modeling. In this thesis it is shown that the results of data analysis vary according to the modeling approach used, thus making the selection of the modeling approach dependent on the purpose of the model. If the model is intended to provide accurate predictions, the approach should be different than in the case where the purpose of modeling is mostly to obtain information about the variables and the process. For industrial applicability it is essential that the methods are robust and sufficiently simple to apply. In this way the methods and the results can be compared and an approach selected that is suitable for the intended purpose. Differences in data analysis methods are compared with data from different fields of industry in this thesis. In the first two papers, the multi block method is considered for data originating from the oil and fertilizer industries. The results are compared to those from PLS and priority PLS. The third paper considers applicability of multivariate models to process control for a reactive crystallization process. In the fourth paper, nonlinear modeling is examined with a data set from the oil industry. The response has a nonlinear relation to the descriptor matrix, and the results are compared between linear modeling, polynomial PLS and nonlinear modeling using nonlinear score vectors.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

The Power of Data - Advocating for Library by Showing Evidence of Impact

Markku Laitinen's keynote presentation in the QQML conference in Limerick, Ireland the 23rd of April, 2012.

Immigration Strategy for Northern Ostrobothnia 2015

A strategy process was completed in the ESF project “Promotion of Work-related Immigration”, which was implemented at Centre for Economic Development, Transport and the Environment for North Ostrobothnia, and an immigration strategy was drawn up for Northern Ostrobothnia on the basis of the process. Information was collected about the situation in Northern Ostrobothnia from the point of view of immigration and the future availability of labour. The intention was to use the information as background material for the strategy. Employers’ need for support in recruiting foreign labour was investigated with a broad inquiry, to which 1000 respondents replied. The strategy process was carried out together with an outside consultant (Net Effect Oy) by arranging three workshops and a seminar where the workshop results were summarised. A large number of companies, authorities, municipalities, associations, project actors and immigrants engaged in immigration issues participated in the workshops. The draft strategy is based on their experiences about immigration and on statistical data, background inquiries and surveys. To ensure the accuracy of the draft strategy, comments were requested from several parties and received from 64 organisations. The core of the immigration strategy consists of an initial analysis, values, a vision and priorities. The strategy is composed of three priorities. The key aim of the priority Internationalisation and Supporting Diversity is to support diversity in schools, workplaces and people’s everyday lives e.g. through attitude development and by promoting internationalisation in companies and education institutions. The aim of the priority Supporting Entrepreneurship and Recruiting Foreign Labour is to promote entrepreneurship among immigrants and the recruitment of foreign labour and to develop the forecasting of educational needs. The priority Developing Integration Services, Regional Cooperation and Networks, in turn, seeks to develop the service structure and policies of immigrant integration and to increase cooperation and exchange of information between regional actors engaged in integration issues. The aim is to use the strategy as a guideline document for immigration issues in Northern Ostrobothnia. The strategy is used to coordinate the existing organisations and operations dealing with immigration issues. In addition, it contains a future-oriented focus and underlines the management of new immigration projects and operations. The main party responsible for the implementation of the strategy is the Immigration Committee. In addition, responsible parties have been assigned to each measure. The implementation of the immigration strategy will be monitored annually on the basis of indicators.